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FOREWORD 


This document is submitted in accordance with the requirements of NASA 
Contract NAS1-14946, Modification Number 3, Study of the Prediction of 
Cruise Noise and Laminar Flow Control Noise Criteria for Subsonic Air 
Transports . 

The final technical reports of this program comprise three volumes. The 
first volume, CR-159104, describes the technical selection and development 
of cruise noise prediction and LFC noise criteria procedures. The second 
volume, CR-159105, is a Methods Manual for applying the cruise noise pre- 
diction procedures to actual aircraft designs. The final volume, this 
report, consists of a User's Guide for a computer program developed from 
the noise prediction procedures described in CR-159105. 
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USERS GUIDE FOR NEAR-FIELD CRUISE NOISE PREDICTION PROGRAM 


SUMMARY 

A computer program for the prediction of near-field aircraft noise under 
cruise conditions has been developed. The total noise for free-field, loss- 
less conditions at selected observer locations is obtained by summing the 
contributions from up to nine acoustic sources. These noise sources, 
selected by the user, include the: fan/compressor, turbine, core (com- 

bustion), jet, shock, and airframe (trailing-edges and turbulent boundary 
layers). The effects of acoustic suppression materials (e.g. engine inlet 
treatment) may also be included in the noise prediction. The program is 
available for use on the NASA/Langley Research Center CDC computer. This 
document provides detailed instructions for using this program, a program 
listing, and a sample case with output. A few comparisons of the programs 
predictions with measured data are given and some possible reasons for their 
lack of agreement presented. 


INTRODUCTION 

The overall objective of NASA Contract NAS1-14946 (A Study of the Prediction 
of Cruise Noise and Laminar Flow Control Noise Criteria for Subsonic Air 
Transports) was to develop procedures for identifying exterior surfaces on 
laminar flow control (LFC) aircraft which are subject to noise levels which 
would adversely affect the operation of LFC. Interest in the development 
of this technology was precipitated by the desire to develop LFC as one 
method for improving the fuel efficiency of commercial air transports. Two 
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of the primary tasks required to accomplish this objective were (1) to 
develop from existing technology general procedures for the prediction of 
aircraft noise incident upon aircraft surfaces during subsonic cruising 
flight, and (2) to summarize and explicitly define all the prediction 
methods in a Cruise Noise Prediction Methods Manual. The results of this 
work are documented in References 1 and 2. This Users Guide represents 
work done under an extension to the basic contract to develop a Cruise 
Noise Prediction Computer Program from the Methods Manual procedures of 
Reference 2. It includes a sample of the comparisons made between noise 
levels measured on the NASA/Dryden JetStar aircraft and noise levels pre- 
dicted by the program. Where possible, the prediction method was to have 
been modified and improved as a result of these comparisons. However, the 
program has not been modified as a result of these comparisons since those 
factors responsible for lack of agreement between measured and predicted 
noise levels have not been identified. This is discussed further with 
the presentation of the comparisons. 


PREDICTION METHODS 

The noise prediction methods employed in this program were selected or 
developed from currently available technology. They are believed to be 
the best available methods for application to aircraft noise prediction 
in the near field during cruising flight. It should be pointed out, how- 
ever, that the term "near field" as applied here actually refers to the 
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close-in acoustic far field . Discussions of the aircraft noise sources 
considered and the selection and development of the prediction procedures 
used herein are contained in Reference 1. The details of the resulting 
computational procedures are given in Reference 2 in the form of a "Methods 
Manual." As noted in both References, no methods were found in the avail- 
able literature which were directly applicable to near-field noise predic- 
tion in the cruise regime. In the case of prediction methods for propul- 
sion system (or LFC suction system) noise sources, those methods available 
have typically been derived for the acoustic far field at sea level con- 
ditions with no forward speed or for the relatively slow speeds associated 
with takeoff or landing. This is especially true of the turbomachinery 
noise sources: fan/compressor, turbine, and core. For these sources, in 

particular, it is believed that the available methods will produce satis- 
factory predictions for the close-in far field since these sources are 
expected to behave essentially as acoustic point sources. Consequently, 
the basic far-field prediction procedures for the turbomachinery sources' 
have been applied here with modifications as outlined below which attempt 
to adapt the methods to the cruise environment. 


The acoustic near-field is generally considered to be the region within 
a couple of wavelengths of the noise source or the region where motion 
in the medium is strongly influenced by the source. The close-in far- 
field is used here to describe the region which lies just beyond the near- 
field and extending into the far field. This is the primary region of 
interest for cruise noise predictions. In this region, unlike the far- 
field, the physical dimensions or physical distribution of the noise source 
is important and must be considered. These effects are discussed in some 
detail in Reference 1. 
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The other propulsion system noise sources are the jet, jet shock, and shock 
screech. In the case of jet mixing noise, the procedure used is a near- 
field predictor but is strictly only applicable to a single circular nozzle 
under static conditions. It has been modified as suggested in the litera- 
ture (see Reference 1) for application to coaxial jets and to account for 
forward-speed effects on noise source strength. The procedure employed 
for the prediction of jet shock-associated broadband noise is another far- 
field prediction method. This source is a distributed source, neverthe- 
less the far-field procedure utilized is considered the best available 
for the present case. It has shown reasonable agreement with measured 
near-field data (see Reference 2). The method used to predict shock 
screech (discrete tone) noise is one developed partially from a near-field 
data base. See Reference 1. Finally, for the propulsion noise sources, 
the directivity indices of the basic prediction procedures have been 
extrapolated where necessary to cover the full range from 0° to 180° 
from the engine inlet. 

The available airframe noise source prediction procedures were less well 
defined than those of the propulsion system sources. Consequently, for 
the airframe sources considered, it has been possible to develop the 
evolving prediction technologies into methods somewhat more applicable to 
the near-field case. Conversely, the prediction procedures are less 
general than the propulsion source predictions. The trailing-edge noise 
prediction procedure is a far-field predictor modified in the literature 
for predictions in the close-in far field. In this procedure, the observer 
location is restricted to a plane perpendicular to the trailing-edge at its mid- 
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point. For this reason, most of the trailing-edge noise inputs are a 
function of observer location. Unfortunately, further development of this 
method to handle other observer locations was beyond the scope of the 
contracted effort. The other airframe source included in the present pro- 
cedure is turbulent boundary layer noise. The prediction procedure was 
developed from methodologies available in the literature as described in 
Reference 1. This widely distributed noise source must be subdivided 
into smaller elements to adapt its prediction procedure to the close-in 
far field. As a result, considerable input may be required for this 
source prediction in the form of observer coordinates for each element as 
a function of observer location. 

The prediction methods outlined above are applied to the high altitude, 
high forward speed case of the program presented here by accounting for 
the following effects: 


Forward- 

-speed propagation effects 


(a) 

source emission angle 1 

forward-speed 

(b) 

sound propagation path length j 

t 

transformation 

(c) 

convective amplification 


(d) 

dynamic amplification 



o Relative ve locity/ forward-speed effects on noise source strength 
o High altitude effects (acoustic impedance and atmospheric 
attenuation) 
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Forward-speed propagation effects are discussed in detail in Reference 1. 
Nonetheless, it is noted here that a Doppler frequency shift correction 
is not applied or required since the source and receiver move together. 

A convective amplification correction is applied to all sources except 
jet mixing noise, which has this effect built into its prediction model. 

A correction for dynamic amplification is included for jet mixing noise, 
shock broadband noise, and turbulent boundary layer noise. 

Forward speed and/or relative velocity effects on source strength are in- 
cluded in the selected jet-mixing noise prediction procedure. Otherwise, 
a forward speed correction to sound power is not applied to the other pro- 
pulsion sources. In the case of turbomachinery noise, forward speed is 
not expected to affect source strength since the noise generating mechan- 
ism occurs internal to the engine. Available data for shock broadband 
noise indicate no relative velocity effect on sound power for this source. 
Some recent data show a possible reduction of shock screech noise level 
with forward speed. However no correction is applied here since pre- 
dictions developed from the measured data appear inconsistent. For the 
airframe noise sources, forward speed is, of course, a key factor in the 
basic prediction equations and no separate correction is required. Finally, 
it is noted that forward speed does impact all engine cycle parameters 
including those employed to predict propulsion system noise and thus its 
affect on sound power (and spectral content) is accounted for in this 
way . 
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The effects of high altitude have been accounted for in the form of a 
correction to the local acoustic impedance. A correction for atmospheric 
attentuation is not included in the program since its effect is estimated 
to be small. Neglecting this radiation effect will result in slightly 
conservative noise estimates. The prediction procedures are only appli- 
cable to free-f ield conditions. Surface effects such as reflection, or 
refraction or scattering effects such as in boundary layers, jet exhaust 
wakes, or airfoil wakes are not included in the procedures or methods of 
this program. Shielding effects are not included in the program although 
a method for their estimation was included in References 1 and 2. The 
reason for their omission is discussed with the presentation of the data 
comparisons . 

Finally, the effects of acoustic sound suppression materials may be included 
in the noise analysis. Acoustic treatment can be accounted for in the form 
of duct liners and/or acoustic splitter rings in the engine inlet, or 
treated duct liners in a fan exhaust duct or the primary exhaust duct. 
Attenuation effects can thus be included in the fan/compressor, turbine, 
and/or core noise predictions. The method used is empirical and based on 
extensive test data. The result is a procedure which avoids details of 
the acoustic liner design. The liners, where used, are assumed to be 
acoustically optimized for maximum absorption at a specified design fre- 
quency. The effect of high Mach number inlets on forward-radiated fan 
noise and compressor noise can also be included. 
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SOURCE SELECTION GUIDELINES 


The sources to be included in the noise prediction is a decision which 
rests primarily on the judgment of the user. Any number or combination of 
noise sources is permissible. Some guidelines and considerations for this 
source selection follow: 

Fan/ compressor : This source will normally be included if propulsion system 

noise is desired. To predict the compressor noise of a turbojet engine or 
the forward-radiated fan noise of a turbofan engine the user should specify 
SRC(l) = 1. The "compressor" noise prediction option (SRC(3)) is included 
with some reservations to provide an estimate of the noise radiated by the 
compressor of a turbofan engine. The prediction procedure used is identical 
to that of a first-stage fan with inlet guide vanes. Since the fan noise 
procedure does not account for any effect of sound propagation through 
upstream fan stages, the compressor noise predictions obtained here should 
be considered very preliminary. This stipulation is not expected to 
seriously limit program usefulness, however, since fan noise generally 
dominates the compressor noise of turbofan engines. If a compressor noise 
source is selected by the user, note that the combination tone noise 
component may be deleted from the prediction. This option was provided as 
an investigative tool since the literature indicate that combination tones 
may not propagate through an upstream rotor or stator vane array. Never- 
theless, the method will compute combination tones for first-stage fans 
with or without inlet guide vanes. The fan noise prediction is not recom- 
mended for turbofan engines having more than two fan stages. 
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The selected method is based on a somewhat limited range of fan types. It 


is most applicable to single stage fans of relatively high bypass ratio and 
multiple (say three or more) low pressure turbine stages. This is consistent 
with the turbofan engines expected for the 1985-1990 time frame. Further- 
more, the fan noise prediction method assumes aerodynamically clean, relatively 
short-duct nacelles with hard walls. 

It should also be noted that this prediction procedure has been developed 
for cruising flight, and engine inflow turbulence is assumed to be low. 
Significant inflow turbulence, as might exist during ground static engine 
operation, can be expected to result in fan/compressor noise levels some- 
what higher than will be predicted here. The source prediction method can 
be consulted for a suitable correction factor if desired. 

Turbine noise : As with fan/compressor noise, this source would normally 

be included in predictions where engine noise is of interest. The method 
applied here is not recommended for engines whose low pressure turbines 
contain less than three stages. It may be useful to note that the turbine's 
blade passage frequency can normally be expected to fall beyond the frequency 
range of the predictions (11220 Hz). As a result, the turbine discrete tone 
and peak broadband level will not be included in the predicted spectra. 

Core (combustion noise) : The core noise source is recommended for all cases 

where engine noise may be important. The procedure is applicable to turbojets, 
turbofans, and turboprops. It includes transmission losses through the 
turbine stages. 
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Jet mixing noise: This source should normally be included where propulsion 


system noise is of interest. The program will handle single-flow circular 
exhaust or dual-flow coannular exhaust nozzles. Mixer nozzles are treated 
as single-flow exhaust where the flow parameters are those of the mixed 
exhaust flow. Some care is required when selecting observer locations. 

The program will terminate if the observer location is more than 30 nozzle 
diameters from the nozzle exit plane or from the engine centerline. The 
observer must also be located outside of the exhaust flow region defined 
by a 7.5° angle from the nozzle exit lip. Finally, the observer must be 
located forward of the engine inlet, if the computed source location (a 
function of frequency) relative to the engine centerline is greater than 
the actual observer distance from the engine centerline. 

Shock broadband noise : The procedure is applicable to the shock broadband 

noise associated with the supersonically underexpanded jet exhaust of 
conical nozzles. With appropriate modifications to the prediction parameters, 
it is assumed to also apply to the annular jet exhaust of dual-flow nozzles. 
The method is not applicable to source emission angles of more than 150° 
relative to the engine inlet axis. If this source is specified for a jet 
with a subsonic flow exhaust, the shock noise component will be set to zero. 

Shock discrete tone ( screech) noise : This source, which is believed to 

result from a feedback mechanism between shocks in a jet exhaust and the 
nozzle lip, may be neglected for nozzle pressure ratios less than 2. If 
requested by the user at lower pressure ratios, the program will set the 
predicted levels to zero. The decision of whether or not to include this 
source will most likely be the most difficult. A most important aspect 
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of "screech" noise is that it does not always appear to be present and no 
criteria are available to determine its presence (or absence). It appears 
that this source, if present, may be eliminated by a minor nozzle redesign. 

To be on the conservative side the user may desire to include this source 
in the predictions. If this is done, the noise levels predicted for 
screech should be considered as upper bounds for the true values. In 
any event, noise levels are predicted for only the fundamental and second 
harmonic discrete frequencies. For separate-flow nozzles, the procedure 
is applicable to either or both of the exhaust flows. 

Trailing-edge noise : This source may be included if airframe noise is of 

interest to the user. The method is presumed to be applicable to "turbulent" 
or partially laminar-flow-controlled surfaces as long as the boundary layer 
at the trailing edge consists of fully developed turbulent flow. In its 
present form, the procedure used is somewhat cumbersome when applied to 
normal prediction situations. There is a restriction that the observer 
location must be located in a plane perpendicular to the trailing-edge of 
interest at its mid-point. As a result, the input for trailing-edge noise 
predictions is required at each observer location. This may frequently 
result in extensive input for this source. If the peak frequency is com- 
puted to fall beyond 11220 Hz, this source will be set to zero. 

Turbulent boundary layer noise ; The implementation of this method will 
normally involve approximating the surface of interest (say a portion of the 
fuselage) by one or more flat surfaces subdivided into several elemental 
areas. The center of each element is treated as a point noise source. The 


11 



total noise is estimated by summing, on an energy basis, the noise contri- 
bution from each of these sub-areas. The noise contribution of an element 
whose peak frequency is greater than 11220 Hz will be set to zero. If the 
observer location lies on an element, that element's contribution will also 
be set to zero. In the close-in far field, the total radiated noise levels 
computed for a given surface will be dependent upon the number (and orien- 
tation) of these elemental areas. Unfortunately, no guidance is presently 
available as to a selection criteria which might relate the number and orien- 
tation of the representative elemental areas to the accuracy of the prediction 
procedure. Note, however, that for near-field estimates the more subdivi- 
sions made the more accurate will be the procedure (maximum dimension of the 
elemental area small compared to source-to-receiver distance), although fewer 
"elements" may give a conservative (high) estimate. For far-field estimates, 
one representative segment may be sufficient. 

INPUT DESCRIPTION 

The input data and format required by the Cruise Noise Prediction Program 
(LFCNO) are described in this section. A "bar chart" is used to indicate 
each line of input required. Not all lines are required; the actual input 
required will depend on the sources considered or the type of input selected 
for a given source. The "bar chart" begins in column one (1) for all input. 
The input can be specified in any consistent set of dimensional units. 

Format codes are given as Fa. or la. The "F" indicates a floating point 
value (real number with decimal such as 215.6 or 28.). "I" indicates an 
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integer value (no decimal). The "a" indicates the field or number of 
columns within which the value must appear. It is important to note that 
integer input must appear in the right-most columns within its specified 
field. No spaces are allowed between the fields specified for a line of 
input. The input specification is broken down into four categories: (1) 

general input, (2) configuration inputs, (3) flight condition inputs, and 
(4) observer location inputs. See the input flowchart on the following 
page. The input in all categories except (1) is further broken down by 
noise source. The input specified under these source breakdowns is to be 
included if and only if that source has been specified as a source for 
the prediction. Many other lines of input, as indicated, contain conditional 
input. These lines of input must be omitted if the coded condition, a pre- 
vious input, is not satisfied. 

The control cards for running this program (LFCNO) must include a card to 
preset computer storage allocated for data arrays to zero (LDSET,PRESET= 
ZERO.). The line-by-line input description follows the input flowchart. 
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INPUT FLOWCHART 
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TITLE 


1st 72 Columns 


TITLE - Any title desired by user (2 lines). 


GENERAL INPUT 


SRC ( 1 ) , 

1=1, 

12 

12 

12 

12 

12 


12 


SRC ( I ) , 1=1, 12 


Noise source code; determines noise sources to be 
included in noise prediction. 

SRC ( I ) =00 Do not include noise source I in the 
prediction. 

=01 Include noise source I in noise prediction. 


Noise Source Code: 


I 


= 1 
= 2 
= 3 
= 4 
= 5 
= 6 
= 7 
= 8 
= 9 
= 10 
= 11 
= 12 


fan 


no 1 se 
noi se 
(radiated 

noi se 


from inlet) 


Forward-radiated 
Aft-radiated fan 
Compressor noise" 

Turbine noise 
Core (combustion) 

Jet mixing noise 

Primary-jet shock-associated broadband noise 
Secondary-jet shock-associated broadband noise 
Primary-jet shock discrete tone (screech) noise 
Secondary-jet shock discrete tone (screech) noise 
Trai 1 ing-edge noise 
Turbulent boundary layer noise 


Airframe noise 


Notes : 

(1) Set SRC(3) =02, if a combination tone noise 
component is to be computed for compressor noise. 

(2) "Primary" jet exhaust means single-flow exhaust 
including mixed-flow exhaust nozzles, or the 
primary (or core) flow of a dual-flow exhaust. 
"Secondary jet" means the annular exhaust of a 
dual-flow exhaust, e.g. the fan exhaust of a 
short-duct turbofan. 


See remarks concerning compressor noise prediction on page 8 , 
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TCI 


TC2 TC3 TC4 TC5 


12 

12 



12 

TCI, .. 

. .. TC5 

- 

Acoustic trea 


it code. 

Set =00 if no acoustic treatment suppression effects are 
to be included in the source noise prediction. 


TCI = 01 , 02, or 03 “ For engine inlet acoustic treatment 

effects on fan and/or compressor noise. 

= 01 Engine inlet acoustic treatment 
only 

= 02 High Mach number inlet for 
acoustic suppression 
= 03 Both inlet treatment and high 
Mach number inlet 

TC2 = 01 - For aft-radiated fan noise treatment 

effects. 


TC3 = 01 , 02, or 03 - For compressor noise treatment effects 

= 01 Compressor inlet treatment 
effects only 

= 02 Engine inlet treatment effects 
only 

= 03 Both compressor inlet and 
engine inlet treatment 
effects to be applied 

TC4 = 01 - For turbine noise suppression* 

effects. 

TC5=01 - For core noise suppression* effects 


(Primary exhaust duct treatment) 


RH00 


FI 0. 


CO 


FI 0. 


RH00 

CO 


Sea level, standard day atmospheric density 

Sea level, standard day atmospheric speed of sound 
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P01 


P02 P03 



P01 , P02, P03 - Output codes =00 if the output as specified below is not 

desi red. 


P01 =01 Output 1/3 octave-band sound pressure levels 
(SPL's) and overall sound pressure levels 
(OASPL's) for the total noise of each source 
and the total noise. 

P01 =02 Same as P01 =01, except includes the contri- 
bution (same information as above) of each 
engine at each observer location. 

P02=01 Print out the fan noise component contributions 
of each engine at each observer location. 

P03 = 01 Print out the compressor noise component con- 
tributions of each engine at each observer 
location . 


ENGINE AND AIRCRAFT CONFIGURATION AND DESIGN 
POINT INPUTS PLUS REFERENCE PARAMETERS 


NC 

13 


NC 


Engine exhaust nozzle configuration code. 

=001 Single-flow exhaust 
=002 Dual -flow exhaust 


DX 1 
FI 0. 


DX1 - Overall or maximum engine nacelle length, i.e. 

distance from primary nozzle exit to the inlet. 
See Figure 1 . 
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Input the following line, If NC=002, 


DX2 
FI 0. 


DX2 


Distance between primary and secondary nozzle exit 
planes. See Figure 1. 


FAN INPUTS [SRC (l ) OR SRC(2)=01] 


NS IGV 



NS 


Number of fan stages (maximum of 2) 


IGV - Code for fan first-stage inlet guide vanes 

= 00 No inlet guide vanes 
= 01 Inlet guide vanes present 


Input following two lines for each fan stage, 1=1, NS: 


MTRD ( I ) 
FI 0. 


NRB(I) 

NSV(I) 

RSSF(I) 

13 

13 

FI 0. 


Repeat this input for the 
second fan stage if NS = 2. 


MTRD ( I ) 

NRB(I) 

NSV(I) 

RSSF(I) 


Rotor tip relative Mach number at engine design point (ratio 
of the flow velocity relative to the fan tip to the local speed 
of sound at the inlet to the fan stage). 

Number of rotor blades. 

Number of stator vanes. 

Rotor-stator spacing (percent) see Figure 2. 
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Input the following line, if SRC (2) = 01 and TC2 = 01 , 
XLFA HDFA FDFA 



XLFA - Effective length of fan exhaust duct acoustic treatment 

HDFA - Effective duct height between opposite acoustic liner faces 

FDFA - Acoustic liner design frequency (peak attenuation frequency) (Hz) 


Input the following line for fan and/or compressor noise if TCI =01 or 03; 
and i f SRC (1 ) = 01 , or i f SRC (3) > 0 and TC3 > 1 , 


XLFF HDFF FDFF NSFF 



XLFF - Effective length of engine inlet acoustic treatment 
HDFF - Effective duct height 'between opposite acoustic liner faces 
FDFF - Acoustic liner design frequency (peak attenuation frequency) (Hz) 
NSFF - Number of acoustic splitter rings in the inlet 


COMPRESSOR NOISE INPUTS [SRC(3)=01 OR 02] 


The following input applies to the first compressor stage only. 


MTRDC 



MTRDC - Rotor tip relative Mach number at engine design point. 
(See definition of MTRD(l) above.) 


NRBC NSVC RSSC 



NRBC - Number of rotor blades 
NSVC - Number of stator vanes 

RSSC - Rotor-stator spacing (percent), see Figure 2. 
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Input the following line if TC3=01 or 03, 
XLC HDC FDC 



XLC - Effective length of compressor Inlet acoustic treatment 
HDC - Effective duct height between opposite acoustic liner faces 
FDC - Acoustic liner design frequency (Hz) 


Input the following reference values if either fan or compressor noise is 
to be computed , i.e. S RC ( 1 ) =01,, or SRC (2) =01, or SRC (3) = 01 or 02. 


DTREFF MFREF RREFF 



DTREFF - Reference fan stage temperature rise, 44.4 C° (80 F°) 

MFREF - Reference mass flow rate, 142.9 kg/s (315 lbm/s, 9.79 slug/s) 

RREFF - Reference source-to-observer distance, 1 m (3 - 281 ft, 39-37 in) 


TURBINE NOISE INPUTS [SRC(4)=01] 



NTB - Number of rotor blades in last turbine stage. 
DTT - Blade tip diameter of last turbine stage 
AP - Primary or core nozzle exit area 


UREFT AREFT 



UREFT - Reference tip speed, 340.3 m/s (1116.4 f/s) 

AREFT - Reference nozzle area, 0.0929 m 2 (1.0 f 2 , 144 in 2 ) 
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Input the following line if TC4=01, 


FDT 
FI 0. 


FDT - Turbine acoustic treatment design frequency (peak attenuation 
frequency) (Hz) 


CORE (COMBUSTION) NOISE INPUTS [SRC(5)=01] 


DTTD WREF DTREFC 


FI 0. 


FI 0. 


FI 0. 


DTTD 

WREF 

DTREFC 


Total temperature drop across the turbine (all stages) at the 
engine design point. 

Reference combustor airflow, 15-9 kg/s (35 lbm/s, 1.09 slug/s) 
Reference combustor temperature rise, 777-8 C° (1400 F°) 


Input the following line if TC5 = 01 , 

FDCR 
FI 0. 

FDCR - Core (combustion) noise acoustic treatment design frequency (Hz) 


Input the following two lines (as indicated) if either turbine [SRC(4)=01] 
or core noise [SRC (5) =01] i s to be computed: 


RREFTC 
FI 0. 


RREFTC - Reference source-to-observer distance, 70.4 m (230.9 ft) 


21 







If TCk =01 , or TC5 =01, 


XLTP HDTP 



XLTP - Effective length of acoustic treatment in core exhaust nozzle for 
turbine and/or core noise suppression 

HDTP - Effective duct height between opposite acoustic liner faces 


JET MIXING NOISE INPUTS [SRC (6) =01] 



D1 - Primary exhaust nozzle exit diameter 

TREF - Reference exhaust total temperature, 555. 6°K (1000°R) 


NPLUG 



NPLUG - Primary nozzle code 

= 00 No plug in primary nozzle 
= 01 If primary exhaust is a plug nozzle 


Input the following line, if NPLUG = 01 


HAP 



HAP - Annulus height of primary plug nozzle 
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Input the following line, if NC = 002 



D2 - Diameter of secondary jet nozzle exit 

AR - Ratio of secondary jet exit area to primary jet exit area 

RW - Ratio of secondary jet mass-flow rate to primary jet mass-flow rate 


SHOCK NOISE INPUT 


Input the following line, if SRC(7) =01 or SRC(9) =01, 



DEI - Primary nozzle exit equivalent diameter 


Input the following line, if SRC (8) =01 or SRC(10) =01, 



DE2 - Secondary nozzle exit equivalent diameter 
HA - Secondary nozzle exit annulus height 


Input the following line, if SRC(9) =01, 



A1 - Primary nozzle exit area 
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Input the following line, if SRC (10) = 01 , 



A2 - Secondary nozzle exit area 


Input the following line, if SRC(9) =01 or SRC(lO) =01, 
AREFSS 



AREFSS - Reference nozzle exit area, 2.60 (10) 3 m 2 (0.028 f 2 , 4.03 in 2 ) 


Input the following 1 ine, if SRC (6) =01 or SRC(7) =01 or SRC(8) =01 , 



RG - Gas constant, 286.95 N. m/kg. °K (53-34 1 bf. f/1 bm °R, 1716.16 1 bf . f/slug. °R) 
GC - Dimensional constant, 1.00 kg.m/N.s 2 (32. 174 1 b m . f/1 bf . s 2 , 

1.00 s lug. f/1 bf . s 2 ) 

TRAI L I NG-EDGE NOISE INPUTS [ SRC (11) = 01 ] 

DELC VREF 



DELC - Code for turbulent boundary layer thickness at trailing edge (<$ te ) 
= 00 Estimate 6 te from equation for the turbulent boundary layer 
thickness on a flat plate 

= 01 6 te to be input for each observer location 
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TURBULENT BOUNDARY LAYER NOISE [SRC (12) =01] 


NJ DELSC 



NJ - Number of elemental areas into which surface being considered is 
subdivided for computational purposes, maximum value =20. 

DELSC - Code for turbulent boundary layer displacement thickness (S*) 

= 00 Estimate 6* from equation for turbulent boundary layer 
thickness on a flat plate 

= 01 5* to be input for each elemental area and for each 

observer location 


Note: The aircraft surface area being considered for turbulent 

boundary layer noise generation will normally be subdivided 
into elements. Each element is considered as a flat surface 
and a point source. This should be considered in selecting 
the number and orientation of the surface elements. 


AT(J), J 

= 1, NJ 


FI 0. 

FI 0. 



FI 0. 

5 values per line 


AT(J) - Surface area of the surface element, see Figure 3. 


UREFBL AREFBL RREFBL 


FI 0. 


FI 0. 


FI 0. 


UREFBL 

AREFBL 

RREFBL 


Reference freestream velocity, 236 m/s (774.3 f/s, 458.7 knots, 
527-9 MPH) 

Reference element area, 0.0929 m 2 (1.0 f 2 , 144 in 2 ) 

Reference source-to-observer distance, 1.0 m (3.281 ft, 39.37 in) 
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Input the following line, if-DELSC=00, 


LT (J) , J = 1 , NJ 


FI 0. 

FI 0. 



FI 0. 


5 values per line 


LT(J) - Distance from point at which boundary layer growth begins to center 
of J th surface element. See Figure 3. 


INPUT FOR DESIRED FLIGHT CONDITIONS 

NFLTC 

13 


NFLTC - Number of flight conditions 


* * 

NOTE : Repeat all remaining input for eaah flight condition. 

* ’ * 


ALT 


MA 


PSC 


FI 0. 


FI 0. 


FI 0. 


ALT - Flight altitude* 

MA - Flight Mach number 
PSC - Power setting code* 

"For output identification purposes only. 


RHOA 


CA 


TA 


FI 0. 


F10. 


FI 0. 


RHOA - Ambient air density 

CA - Ambient speed of sound 

TA - Ambient temperature (absolute) 
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FAN NOISE INPUTS [SRC(l)=01 or SRC (2) -01] 


MFF FRPH 



MFF - Mass flow rate through the fan 

FRPM - Fan rotational speed (RPM) 


Input the following line for each fan stage (where 1=1, NS) 


DTSG(I) MTR(I) MT( I ) 

Repeat this line for the 
second stage if NS = 2. 

DTSG(l) - Total temperature rise across the I 11 * 1 fan stage 
MTR(l) - Rotor tip relative Mach number for the 1^ fan state 
[see definition of MTRD(l) above] MTR(l) <MTRD(l) 

MT(l) - Rotor tip Mach number for the |th fan stage (i.e. ratio of 
tip speed to local speed of sound) 


Input the following line, if SRC(2) =01 and if TC2 = 01 , 


MDFA CDFA 



MDFA - Mean flow Mach number through the treated portion of the fan 
exhaust duct 

CDFA - Local speed of sound in the treated portion of the fan exhaust duct 
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Input the following 1 ine, if TCI =01 or 03; and SRC ( 1 ) = 01 , or 
SRC (3) >0 and TC3 > 1 , 


MDFF CDFF 


FI 0. 


FI 0. 


MDFF - Mean flow Mach number through the treated portion of the engine inlet 
CDFF - Local speed of sound in the treated portion of the engine inlet 


Input the following 1 i ne, i f TCI > 1 ; and SRC ( 1 ) = 01, or SRC (3) >0 
and TC3 > 1 , 


Ml 


FI 0. 


Ml - Engine inlet throat Mach number 


COMPRESSOR INPUTS [ SRC (3) = 01 OR 02] 


The following input applies to the first compressor stage: 


MFC CRPM 


FI 0. 


FI 0. 


MFC - Mass-flow rate through the compressor 
CRPM - Compressor rotation speed (RPM) 


DTC 


MTRC 


MTC 


FI 0. 


FI 0. 


FI 0. 


DTC - Total temperature rise across the first compressor stage 
MTRC - Rotor tip relative Mach number [see definition of MTRD(l) above] 
MTRC < MTRDC 

MTC - Rotor. tip Mach number [see definition of MT ( I ) above] 
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Input the following line, if TC3 = 01 or 03, 
MDC CDC 



MDC - Mean flow Mach number through the treated section of the compressor 
inlet 

CDC - local speed of sound in the treated section of the compressor inlet 


TURBINE INPUTS [SRC(4)=01] 



PR - Low-pressure turbine (LPT) total-to-static pressure ratio (i.e. ratio 
of total pressure at LPT inlet to the turbine exit static pressure. 
TRPM - Low-pressure turbine rotational speed (RPM) 


CORE INPUTS [SRC (5) =01] 

W3 RH03 T3 Tfr 

FI 0. FI 0. FI 0. FI 0. 

W3 - Combustor air mass-flow rate 
RH03 - Combustor inlet air density 
T3 - Combustor inlet total temperature 

T4 - Combustor exit total temperature 


Input the following line, if TC4=01 or TC5-01, 


MDTP CDTP 



MDTP - Mean flow Mach number through the acoustically treated portion of 
the engine core exhaust duct 

CDTP - local speed of sound associated with MDTP above 
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JET MIXING AND SHOCK NOISE INPUTS 


Input the following line, if SRC (6) =01, or SRC (7) =01, or SRC (9) =01, 


VI 

FI 0. 


VI - Primary jet exhaust velocity (fully expanded flow) 


Input the following line, if SRC (6) =01 and NC = 002, or SRC (8 ) =01, 
or SRC(10) =01 , 


V 2 


FI 0. 


V2 - Secondary jet exhaust velocity (fully expanded flow) 


Input the following line, if SRC(6)=01, or SRC(7)=01, 


T1 

FI 0. 


T1 - Primary exhaust total temperature (absolute) 


Input the following line, if SRC(6)=01 and NC=002, or SRC(8)=01, 

T2 
FI 0. 

T2 - Secondary exhaust total temperature (absolute) 
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Input the following line, if SRC ( 9) =01, 


NPR1 
FI 0. 


NPR1 - Primary exhaust nozzle pressure ratio. 


Input the following line, if SRC(lO) =01, 

NPR2 
FI 0. 

NPR2 - Secondary exhaust nozzle pressure ratio. 


TURBULENT BOUNDARY LAYER NOISE INPUTS [SRC (12) =01] 


Input the following line, if DELSC = 01 , 


DELS (j) , J = 1 , NJ 

FI 0. 

FI 0. 


FI 0. 



5 values per line 


DELS(J) - Turbulent boundary layer displacement thickness at the 
element of the surface to be considered. 


TEN OR TBLN INPUTS 

Input the following line, if SRC (11) =01 and DELC=00, or 
SRC ( 1 2) =01 and DELSC = 00, 

MUA 
FI 0. 

MUA - Ambient viscosity 
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INPUTS FOR OBSERVER LOCATIONS 


NLOC 

13 


NLOC - Number of observer locations to be considered for this flight case 

* * 

NOTE: Repeat all remaining input for each observer location for 

the specified flight case . 

* * 

LOCID 

13 


LOCID - Observer location identification number 


TRAI LING-EDGE NOISE INPUTS [ SRC (11 ) = 01 ] 


XLTE 


XTE 


YTE 


FI 0. 


FI 0. 


FI 0. 


XLTE - Length of trailing edge considered at this observer location 
XTE - Distance from center of trailing edge to observer location measured 

in the surface flow direction (usually along the wing chord) 
(positive forward). See Figure 4. 

YTE - Observer location measured perpendicular to the extended wing 
chord line. 


Input the following line, if DELC = 00, 


CBAR 
FI 0. 


CBAR - Wing chord or flow-surface length used to compute turbulent boundary 
layer thickness at the surface trailing edge. 
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Input the following line, if DELC=01, 


DELTA 
FI 0. 


DELTA - Turbulent boundary layer thickness at trailing edge. 


TURBULENT BOUNDARY LAYER NOISE INPUTS [SRC (12) =01] 


XT(J) YT (J) ZT(J) 


FI 0. 


FI 0. 


FI 0. 


J = 1 , NJ 


Repeat NJ times for each element J =1 } NJ 

XT(J), YT(J), ZT(J) - Coordinates of the observer location relative to the 

center of the Jth element of the surface considered 
for turbulent boundary layer noise generation. See 
Figure 3. 


NENG 

13 


NENG - Number of engines considered for this observer location 


Repeat the following two lines for each engine for the last specified 
observer location: 


NOENG(J) 

13 


NOENG(J) - Number identification of the J 1 -* 1 engine considered, where J=1, 
NENG. 


33 








X,Y - Coordinates of observer location relative to engine NOENG(J). 

The origin is at the center of the primary nozzle exit plane. 
X is measured positive forward along the engine centerline. 

Y is measured perpendicular to the engine centerline. 
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FAN DUCT 



Figure 1. Description of geometry input for fan or compressor noise, 
observer location for propulsion noise sources, and 
turbomachinery source locations. 
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RSSF = (c^/c^) x 100 

(a) GEOMETRY FOR 1 ST -STAGE FANS WITHOUT INLET GUIDE VANES. 


STATOR 



RSSF or RSSC = MINIMUM OF: {c^Cy) x 100 OR (c 4 /c 3 ) x 100 

(b) GEOMETRY FOR iST-sTAGE FANS WITH INLET GUIDE VANES (IGV), 
2ND-STAGE FANS, OR 1 ST -STAGE OF COMPRESSORS. 


Figure 2. Fan or compressor rotor-stator geometry. 
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NOTE: TOTAL AREA FOR PREDICTION = ^ AT(J) 

J=1 


Figure 3. Geometry for turbulent boundary layer noise. 
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OBSERVER 



Note: The observer location must lie in a plane perpendicular 

to the trailing edge at its midpoint. 


Figure 4. Source/observer geometry for trailing-edge noise. 
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PROGRAM OUTPUT 


The Cruise Noise Prediction Program (LFCNO) predicts one-third octave band 
sound pressure levels (from 50 to 10,000 Hz) and overall sound pressure 
levels (OASPL's) for up to nine aircraft noise sources and the total noise 
at selected locations on an aircraft surface for each selected flight 
condition. The predicted levels are for free-f ield , lossless conditions. 
Optional output includes the fan or compressor component noise contribu- 
tions at each observer location for each engine. Another option allows 
the user to suppress the output for individual engines and receive only 
the total and total component noise at each observer location for each 
flight condition. 

In most cases, where input data or computed noise prediction parameters 
fall outside allowable or desirable ranges, the program will printout a 
message indicating the exceedance. The computer run may or may not termi- 
nate depending on the problem. Adherence to the limitations and input 
requirements specified in the input description and to the input flowchart 
will avert most problems. Table I indicates those occurrences which will 
result in a diagnostic or warning message. 



TABLE I 


PROGRAM DIAGNOSTICS AND WARNING MESSAGES 


PROGRAM ROUTINE 

PROBLEM 

RUN 

TERMINATED 

MAIN 

Airplane Mach number, M^ ^ 1.0 

Yes 


Blade passage frequency outside 
range of 44.7 to 22390 Hz 

Yes 

FAN 

Rotor tip relative Mach number 
exceeds design value 

Yes 

TURB 

(Turbine) 

Computed peak 1/3 O.B. sound 
pressure level exceeds computed 
overall sound pressure level. 
(This check is a holdover 
from initial program checkout; it 
would indicate a basic problem 
with the prediction equations) 

Yes 


Observer location is outside allowable 

Yes 


range of 30 equivalent nozzle diam- 
eters from either the nozzle exit 



plane or centerline. 


JETMX 



(Jet Mixing 
Noise) 

Observer located within the iet 

J 0 

exhaust region defined by a 7.5 

Yes 


angle from the exhaust centerline 
at the nozzle lip. 



The observer must be located forward 

Yes 


of the engine inlet if the observer 
distance from the engine centerline 
is less than the computed source 
distance from the engine centerline 
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TABLE I - continued 


PROGRAM ROUTINE 

PROBLEM 

RUN 

TERMINATED 

EXTFRQ 

(called by Jet 
Mixing Noise 

The center frequency of one or 
more of the three octave bands 
of the basic procedure fall out- 
side the Stroubal number range 
of the prediction method. 

Yes 

Routine) 

Extrapolation of the frequency range 
to cover 50 to 10,000 Hz exceeds the 
frequency range of the normalized or 
"reduced" (THOMSON) spectrum shapes 
used in the extrapolation. 

No 

SHOCK 

(Shock 

broadband) 

One or more frequencies (50 to 10,000 Hz) 
fall outside the Stroubal number range 
of the prediction procedure. 

No 

XFORM 
( forward- 
speed trans- 
formation 
routine) 

Observer location same as source 
location 

Yes 
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PROGRAM LISTINGS 


42 



MAIN 


1 PROGRAM MAIN ( INPUT r OUTPUT t T APE5= INPUT r T APE6=0UTPUT ) 


2 C • 

3 C DRIVER PROGRAM FOR NEAR-FIELD CRUISE NOISE PREDICTION 

4 C 

5 C LATEST REVISION ON 12-12-79 

6 C 

7 REAL MFF r MFC r MFREF r DTSG ( 2 ) »MTR ( 2 ) t MTRD (2 ) r MT (2 ) t MTRC t MTRDC r MTC t 

8 2 RSSF (2) f MDFF f MDFA r NSFF » MDC > MI > MDTP > TITLE ( 24 ) t 

9 3 SFFBB ( 24 r 4 ) * SFFDT (24 >4) »SFFCTN(24r4) > SFF ( 24 * 4 ) * SFABB ( 24 » 4 ) r 

10 4 SFADT ( 24 1 4 ) rSFA(24>4) t SFAN (24 r 4 ) r OSFFBB ( 4 ) t OSFFDT ( 4 ) t OSCTN ( 4 ) t 

11 5 OASFF (4 ) r OSFABB (4 ) j OSFADT (4 ) * OASFA (4 ) 1 0 ASF AN (4 ) » 

12 6 SCBB (24 r 4 ) t SCDT ( 24 1 4 ) t SCCTN ( 24 r 4 ) fSCMX(24r4) >0SCBB(4) rOSCDT(4 ) t 

13 7 OSCCTN ( 4 ) f OASCMX ( 4 ) > ST ( 24 1 4 ) f OASTUR ( 4 ) r SCR (24 f 4 ) r OASCOR ( 4 ) 

14 8 r SUP (24r 4 ) f S JS (24 » 4 ) > SJMX (24 f 4 ) r OSUP ( 4 ) t OSJS (4 ) r OAS JMX (4) 

15 REAL MAf F (24 ) r NPR1 > NPR2 » MUA » LT (20 ) r AT (20 ) » XT (20 ) t 

16 2 YT (20 ) » ZT (20 ) r DELS (20 ) » STEN (24 ) 

17 REAL SPSBB (24 r 4 ) > SSSBB (24 »4 ) > SBB1 ( 24 ) r SBB2 ( 24 ) t OSPSBB ( 4 ) r 

18 2 OSSSBB ( 4 ) r SSP (24 * 4 ) > SSS ( 24 » 4 ) r OSSSP ( 4 ) p OSSSS ( 4 ) 

19 3 p STOT ( 24 p 4 ) p OASTOT ( 4 ) 

20 REAL SFFT (24) fSFAT(24) pSFANT (24 ) pSCMXT(24) pSTT ( 24) pSCRT(24) p 

21 2 S JPT ( 24 ) p S JST ( 24 ) r S JMXT (24 ) p SPSBBT ( 24 ) t SSSBBT ( 24 ) p 

22 3 SSPT (24 ) p SSST ( 24 ) p STOTAL (24 ) 

23 INTEGER SRC ( 12 ) p ERR p SCI p SC2 t NRB ( 2 ) p NSV ( 2 ) r POl p 

24 2 P02 1 P03 p DELC p DELSC pTC1pTC2pTC3pTC4pTC5pN0ENG(4) 

25 COMMON MApRHOApCApRHOOp CO pCIMPDp DX1 p DX2 p ERR 

26 COMMON/FAN/ IGV p SI ( 24 } p S2 (24 ) p S3 (24 ) p S 4 (24 ) p S5 ( 24 ) p DTREFF p MFREF p 

27 2 RREFF p XLFF p HDFF p MDFF p CDFFp FDFF pNSFFpMIpXLFAp HDFA pMDFApCDFApFDFAp 

28 3 XLC p HDC p MDC p CDC p FDC 

29 COMMON /TC/STC (24) p AREFT p UREFT p WREF p DTREFC p XLTP p HDTP p MDTP p CDTP 

30 COMMON /JMX/SJ1 (24) pSJ2(24 ) pTREFpRGpGCpHAPpRWpD2pV2pT2pAR 

31 COMMON /TBL/UREFBL p AREFBL p RREFBL p STBL (24 ) 

32 C 

33 C PERMANENT DATA 

34 C 

35 DATA F/50. p 63. p80. p 100. p 125. rl60, r200. >250. p315. r400 , p500. p 

36 2 630. p 800. >1000. >1250. p 1600 . p 2000. >2500. p3150 . p 4000. p 

37 3 5000. >6300. >8000. p 10000,/ 

38 C 

39 C 

40 C READ GENERAL INPUT 

41 C 

42 C 

TITLE 

( SRC ( I ) p 1=1 p 12 ) 

TC1pTC2pTC3p TC 4p TC5 
RHOOpCO 
POl p P02 p P03 
NC 
DX1 


43 

READ(5p104) 

44 

READ ( 5 p 100 ) 

45 

READ ( 5 p 100 ) 

46 

READ ( 5 p 103 ) 

47 

READ ( 5 p 100 ) 

48 

READ(5p 105) 

49 

READ(5f 103) 

50 

DX2=0,0 
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MAIN 


51 IF (NC.EQ.2) READ ( 5 r 103 ) DX2 

5? C 

53 ERR=0 

54 C 

55 C INPUT FORMAT STATEMENTS 

56 C 

57 100 FORMAT ( 1212 ) 

58 101 FORMAT ( 5F10 . 0 ) 

59 102 FORMAT (3F10«0 ) 

60 103 FORMAT (8F10.0 ) 

61 104 FORMAT ( 12A6 ) 

62 105 FORMAT (2I3»F10«0) 

63 106 FORMAT ( 13 f 5F10 *2 ) 

64 C 

65 C READ FAN INPUTS 

66 C 

67 IF(SRC(1)«EQ«0* AND. SRC ( 2) • EQ • 0 ) GO TO 6 

68 READ ( 5 » 100 ) NS>IGV 

69 DO 4 1=1 t NS 

70 READ ( 5 r 103 ) MTRD(I) 

71 READ ( 5 r 105 ) NRB ( I ) t NSV ( I ) t RSSF ( I ) 

72 4 CONTINUE 

73 IF(SRC(2) .EQ.1.AND.TC2.EQ.1) READ(5>103) XLFArHDFArFDFA 

74 6 CONTINUE 

75 IF(SRC(1) .EQ.O.OR.TC1.EQ.O) GO TO 7 

76 GO TO 8 

77 7 IF(SRC(3) . EQ. 0 . OR . TC3. LE. 1 ) GO TO 9 

78 8 IF(TC1.EQ.1.0R.TC1.EQ.3) READ(5>103} XLFF , HDFF t FDFF , NSFF 

79 9 CONTINUE 

80 C 

81 C 

82 C READ COMPRESSOR INPUTS 

83 C 

84 C 

85 IF ( SRC ( 3 ) « EQ . 0 ) GO TO 10 

86 READ ( 5 r 103 ) MTRDC 

87 READ ( 5 1 105 ) NRBC t NSVC t RSSC 

88 IF(TC3.EQ.1,0R.TC3.EQ.3) READ(5rl03) XLC t HDC t FDC 

89 10 CONTINUE 

90 C 

91 C REF VALUES FOR FAN OR COMPRESSOR 

92 C 

93 IF(SRC(1) .EQ.1.0R.SRC(2) . EQ. 1 . OR . SRC ( 3 ) .GE.l) READ(5rl03) DTREFF t 

94 2 MFREF t RREFF 

95 C 

96 C READ TURBINE INPUTS 

97 C 

98 IF ( SRC ( 4 ) . EQ • 0 ) GO TO 12 

99 READ ( 5 » 106 ) NTBr DTT t AP 

100 READ ( 5 f 103 ) UREFTfAREFT 
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MAIN 


101 

102 

103 C 

104 C 

105 C 

106 

107 

108 

109 

110 C 

111 C 

112 C 

113 

114 

115 C 

116 C 

117 C 

118 C 

119 

120 
121 
122 

123 

124 

125 C 

126 C 

127 C 

128 

129 

130 

131 

132 C 

133 

134 C 

135 C 

136 C 

137 

138 C 

139 C 

140 C 

141 

142 C 

143 C 

144 C 

145 

146 

147 

148 

149 

150 


IF ( TC4 . EQ. 1 ) READ { 5 f 103 ) FDT 
12 CONTINUE 

READ CORE INPUTS 

IF ( SRC ( 5 ) • EQ . 0 ) GO TO 14 
READ ( 5 » 103 ) DTTD» WREF »DTREFC 
IF ( TC5 • EQ. 1 ) READ (5 » 103 ) FDCR 
14 CONTINUE 

READ REFERENCE DISTANCE FOR TURBINE AND CORE 

IF ( SRC ( 4 ) .EQ.1.0R.SRC{5) .EQ.l) READ(5>103) RREFTC 
IF(TC4. EQ.l. OR ,TC5. EQ.l) READ(5rl03) XLTP r HDTP 


READ JET MIXING NOISE INPUTS 

IF ( SRC (6 ) . EQ • 0 ) GO TO 18 
READ ( 5 r 103 ) D1 1 TREF 
READ ( 5 f 100 ) NPLUG 
IF (NPLUG.EQ.l ) READ (5 » 103 ) HAP 
IF (NC.EQ.2) READ ( 5 » 103 ) D2* AR >RW 
18 CONTINUE 

READ SHOCK NOISE INPUTS 

IF(SRC(7). EQ.l. OR, SRC(9). EQ.l) READ(5*103) DEI 
IF(SRC(8). EQ.l. OR. SRC(IO). EQ.l) READ(5d03) DE2 t HA 
IF(SRC(9) .EQ.l) READ ( 5 » 103 ) A1 
IF(SRC(10) .EQ.l) READ ( 5 » 103 ) A2 

IF(SRC(9) .EQ.l. OR. SRC(IO) .EQ.l) READ(5rl03) AREFSS 

COMMON INPUTS FOR JET MIXING AND SHOCK NOISE 

IF (SRC (6) . EQ • 1 .OR , SRC ( 7 ) . EQ . 1 . OR . SRC ( 8 ) .EQ.l) READ (5 , 103) RG> GC 

READ TR AILING-EDGE NOISE INPUTS 

IF(SRC(11) .EQ.l) READ (5 » 106 ) DELC t VREF 

READ TURB. BOUNDARY LAYER NOISE INPUTS 

IF ( SRC (12 ) . EQ . 0 ) GO TO 22 
READ (5» 100 ) NJ» DELSC 
READ (5» 101 ) ( AT( J) , J=1,NJ) 

READ ( 5 » 103 ) UREFBL » AREFBL »RREFBL 

IF ( DELSC , EQ . 0) READ(5»101) ( LT ( J ) » J=1 t N J ) 

22 CONTINUE 
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MAIN 


151 C 

152 C FLIGHT CONDITION INPUTS = FUNC { ALT , MA , POWER SETTING) 

153 C 

154 READ ( 5 » 105 ) NFLTC 

155 NFC=1 

156 23 READ (5 » 103 ) ALT,MA,PSC 

157 READ (5 » 103 ) RHOA»CA»TA 

158 IF(MA.GE.l.O) GO TO 400 

159 C 

160 C FAN INPUTS 

161 C 

162 IF ( SRC ( 1 ) . EQ. 0 • AND# SRC (2 ) . EQ • 0 ) GO TO 26 

163 READ (5 r 103 ) MFF t FRPM 

164 DO 24 1=1 * NS 

165 READ { 5 1 103 ) DTSG ( I ) »MTR < I > r MT ( I ) 

166 24 CONTINUE 

167 IF(SRC<2) .EQ.l. AND. TC2. EQ.l) READ(5rl03) MDFA t CDFA 

168 26 CONTINUE 

169 IF(SRC(1) .EQ.O.OR.TC1.EQ.O) GO TO 28 

170 GO TO 30 

171 28 IF (SRC(3) . EQ.0.0R.TC3.LE.1) GO TO 32 

172 30 IF(TC1.EQ.1.0R.TC1.EQ.3) READ(5»103) MDFFrCDFF 

173 IF (TC1.EQ.2.0R.TC1.EQ.3) READ(5»103) MI 

174 32 CONTINUE 

175 C 

176 C COMPRESSOR INPUTS 

177 C 

178 IF ( SRC ( 3 ) • EQ . 0 ) GO TO 34 

179 READ ( 5 » l63 ) MFC » CRPM 

180 READ ( 5 » 103 ) DTC»MTRC»MTC 

181 IF(TC3.EQ.1.0R.TC3.EQ.3) READ(5»103) MDC^CDC 

182 34 CONTINUE 

183 C 

184 C TURBINE INPUTS 

185 C 

186 IF (SRC(4) .EQ.l) READ(5»103) PR t TRPM 

187 C 

188 C CORE INPUTS 

189 C 

190 IF(SRC(5) .EQ.l) READ(5»103) W3 » RH03 r T3 » T4 

191 C 

192 C TURBINE OR CORE INPUT 

193 C 

194 IF (TC4. EQ.l. OR. TC5. EQ.l ) READ(5»103) MDTP r CDTP 

195 C 

196 C JETMX AND SHOCK INPUTS 

197 C 

198 IF(SRC(6). EQ.l. OR. SRC<7). EQ.l. OR. SRC(9). EQ.l) READ (5 » 103 ) VI 

199 IF(SRC(6),EQ.l. AND. NC . EQ.2.0R .SRC ( 8 ) . EQ . 1 . OR . SRC ( 1 0 ) .EQ.l) 

200 2 READ ( 5 r 103 ) V 2 
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MAIN 


201 IF (SRC (6 ) • EQ. 1 .OR, SRC (7 ) «EQ« 1 ) READ(5*103) T1 

202 IF(SRC(6) .EQ.1.AND.NC.EQ.2.0R.SRC (8) .EQ.l) READ(5*103) T2 

203 C 

204 C SHOCK INPUTS 

205 C 

206 IF(SRC(9) .EQ.l) READ(5*103) NPR1 

207 IF(SRC(10) .EQ.l) READ(5*103) NPR2 

208 C 

209 C TBLN INPUT 

210 C 

211 IF (DELSC.EQ, 1 ) READ(5*101) ( DELS ( J ) t J=1 * NU ) 

212 C 

213 C READ COMMON INPUTS FOR TEN AND TBLN 

214 C 

215 IF (SRC (11) .EQ.l. AND . DELC . EQ . 0 . OR . SRC (12) . EQ . 1 . AND . DELSC . EQ . 0 ) 

216 2 READ (5*103) MUA 

217 CIMPD=10 . *ALOG10 (RHOA*CA/ (RHOO*CO ) ) 

218 C 

219 C OBSERVER INPUTS 

220 C 

221 READ (5*105) NLOC 

222 NO=l 

223 40 READ (5*105) LOCID 

224 IF ( SRC ( 11 ) . EQ . 0 ) GO TO 41 

225 READ (5*103) XLTE * XTE * YTE 

226 DELTA=0.0 

227 IF (DELC.EQ«0 ) READ(5*103) CBAR 

228 IF ( DELC .EQ.l) READ<5*103) DELTA 

229 41 CONTINUE 

230 IF ( SRC ( 12 ) . EQ . 0 ) GO TO 42 

231 READ (5*102) ( XT ( J ) * YT ( J ) * ZT ( J ) t J=1 * N J ) 

232 42 CONTINUE 

233 READ ( 5 » 105 ) NENG 

234 J=1 

235 45 READ (5*105) NOENG(J) 

236 READ (5 » 103 ) Xt Y 

237 XXrX 

238 YY=Y 

239 C 

240 C 

241 C 

242 C COMPUTE ENGINE COMPONENT NOISE LEVELS 

243 C 

244 C 

245 C 

246 C INITIALIZE COMPUTED DATA STORAGE ARRAYS 

247 C 

248 DO 50 1=1*24 

249 SFFBB ( I * J ) =0 , 0 

250 SFFDT ( I * J ) =0 « 0 
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251 SFFCTNdr J)=0.0 

252 SFF ( I r J ) =0 . 0 

253 SFABBdr J)=0.0 

254 SFADTdf J)=0.0 

255 SFA(I,J)=0.0 

256 SFANdr J)=0.0 

257 50 CONTINUE 

258 OSFFBB ( J) =0 • 0 

259 05FFDT ( J) =0 • 0 

260 OSCTN(J)=0.0 

261 OASFF ( J ) =0 • 0 

262 OSFABB ( J) =0 • 0 

263 OSFADK J)=0.0 

264 0A5FA ( J ) =0 • 0 

265 OASFAN ( J ) =0 • 0 

266 C 

267 IF (SRC ( 1 ) »EQ» 0 • AND.SRC (2) .EQ.O ) 60 TO 60 

268 C 

269 C COMPUTE DESIRED FAN NOISE COMPONENTS 

270 C 

271 SC1=SRC(1) 

272 SC2=SRC(2) 

273 SUM1=0.0 

274 SUM2=0.0 

275 SUM3=0.0 

276 SUM4=0.0 

277 SUM5=0.0 

278 C 

279 C 

280 DO 56 K=1»NS 

281 DELT=DTSG(K) 

282 XMTR=MTR(K) 

283 XMTRD=MTRD(K) 

284 XMT=MT(K) 

285 XNR=NRB(K) 

286 XNS=NSV(K) 

287 RSS=RSSF(K) 

288 NCODE=K 

289 CALL FAN(XXrYYrMFFfFRPMfSClrSC2f NCODEf TClf TC2rTC3f DELTrXMTRr 

290 2 XMTRD t XMT t XNR t XNS > RSS » BPFF r FREFF ) 

291 IF(ERR.EQ.l) GO TO 1000 

292 IF(K.EQ.2) GO TO 51 

293 BPFF1=BPFF 

294 FREFF1=FREFF 

295 GO TO 52 

296 51 BPFF2=BPFF 

297 FREFF2=FREFF 

298 52 CONTINUE 

299 C 

300 DO 55 1=1 > 24 
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301 

302 
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304 

305 C 

306 

307 

308 

309 

310 

311 

312 

313 


314 

315 

316 C 

317 

318 

319 

320 C 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 C 

332 

333 

334 C 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 


IF(SCl.EQ.O) GO TO 53 

IF (SI ( I ) .GT.O.O) SUM1=SUM1+10.**(S1 ( I )/10. ) 

IF(S2(I) .GT.O.O) SUM2=SUM2+10 . **(S2(I)/10.) 

IF ( S3 ( I) .GT.O.O) SUM3=SUM3+10.**(S3( I )/10. ) 

IF (SFFBB (IrJ) .GT.O.O. OR. SKI) .GT.O.O) 

2 SFFBB( If J)=10.*ALOG10(10.**(SFFBB( Ir J)/10. )+10.**(Sl ( I )/10. ) ) 

IF (SFFDT (IrJ). GT.O.O, OR ,S2 (I) .GT.O.O) 

2 SFFDT ( If J)=10.*ALOG10 ( 10 .** (SFFDT ( If J)/10. )+10.**(S2( D/10, ) ) 

IF (SFFCTN (IrJ). GT.O.O. OR. S3 (I). GT.O.O) 

2 SFFCTN (I r J)=10.*ALOG10 ( 10. **( SFFCTN ( Ir J)/10. ) +10 . ** ( S3 (I ) /10 . ) ) 
IF(K.NE.NS) GO TO 53 

IF (SFFBBt IrJ) . GT.O.O. OR. SFFDT ( Ir J) . GT.O.O .OR. SFFCTN (I r J) .GT.O.O) 

2 SFFdr J)=10.*ALOG10(10.**(SFFBBdr J)/10. ) + 

3 10.** (SFFDT (IrJ)/10. )+10. ** ( SFFCTN ( I r J ) /10 , ) ) 

53 IF (SC2.EQ.0 ) GO TO 54 

IF (S4 (I ) .GT.O.O) SUM4=SUM4+10.**(S4(I)/10. ) 

IF(S5( I) .GT.O.O) SUM5=SUM5+10.**(S5(I)/10. ) 

IF(SFABBdrJ) .GT.O.O. OR. S4 (I) .GT.O.O) 

2 SFABB dr J)=10.*ALOG10 ( 10 . ** ( SFABB ( If J)/10. )+10.**(S4( D/10, ) ) 

IF ( SFADT (IrJ). GT.O.O. OR .55 (I). GT.O.O) 

2 SFADT (Ir J ) =10 . *ALOG10 ( 10.** (SFADT (IrJ) /10. ) +10.** (S5( D/10. ) ) 

IF (K.NE.NS) GO TO 54 

IF (SFABB (IrJ).GT.O, 0 .OR .SFADT (IrJ) .GT.O.O) SFAdr J) = 

2 10.*ALOG10(10.**(SFABB(Ir J)/10. ) +10 . ** ( SFADT (I r J)/10. ) ) 

54 IF (K.NE.NS) GO TO 55 

IF (SFF (IrJ) .GT.O.O.OR.SFAdr J) .GT.O.O) SFANdr J) = 

2 10.*ALOG10(10.**(SFF(Ir J)/10. )+10.**(SFA(Ir J)/10. ) ) 

55 CONTINUE 

56 CONTINUE 

IF(SCl.EQ.O) GO TO 57 

IF (SUM1. GT.O.O) OSFFBB ( J ) =10 • *ALOG10 ( SUM1 ) 

IF (SUM2. GT.O.O) OSFFDT ( J ) =10 . *AI_OG10 ( SUM2 ) 

IF (SUM3. GT.O.O) OSCTN ( J ) =10 . *ALOG10 ( SUM3 ) 

SUMX=SUM1+SUM2+SUM3 

IF (SUMX. GT.O.O) OASFF ( J ) =10 , *ALOG10 ( SUMX ) 

I F ( SC2 , EQ , 0 ) 0 ASF AN ( J ) =OASFF ( J ) 

57 IF ( SC2 • EQ , 0 ) GO TO 58 

IF (SUM4. GT.O.O) OSFABB ( J) =10 • *ALOG10 (SUM4 ) 

IF ( SUM5 .GT.O.O) OSFADT ( J ) =10 . *ALOG10 ( SUM5 ) 

SUMY=SUM4+SUM5 

IF (SUMY. GT.O.O) OASFA ( J ) =10 . *ALOG10 ( SUMY ) 

IF(SCl.EQ.O) OASFAN ( J ) =OASFA ( J ) 

SUMXY=SUMX+SUMY 

IF ( SCI. EQ.l. AND. SUMXY. GT.O.O) OASFAN ( J ) =1 0 . *AL0G1 0 ( SUMXY ) 

58 CONTINUE 
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353 C 

354 C 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 C 

366 C 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 C 

394 C 

395 C 

396 

397 

398 

399 

400 C 


60 CONTINUE 

COMPUTE COMPRESSOR NOISE 

DO 61 1=1,24 
SCBBdr J)=0.0 
SCDT(Ir J)=0.0 
SCCTNU, J)=0.0 
SCMX( I, J)=0.0 

61 CONTINUE 
OSCBB( J)=0.0 
OSCDT( J)=0.0 
OSCCTN ( J ) =0 , 0 
OASCMX ( J ) =0 « 0 


IF (SRC (3) .EQ.O ) GO TO 70 

XNRC=NRBC 

XNSC=NSVC 

NC0DE=2+SRC(3) 

CALL FAN ( XX , YY , MFC , CRPM »0,0,NCODE,TC1,TC2,TC3,DTC, MTRC , MTRDC , 

2 MTC , XNRC » XNSC » RSSC , BPFC , FREFC ) 

IF(ERR.EQ.l) GO TO 1000 

SUM1=0.0 

SUM2=0.0 

SUM3=0.0 

DO 62 1=1,24 

IF (SI (I) .GT.0,0) SUM1=SUM1+10.**(S1(I)/10. ) 

IF ( S2 ( I ) .GT.0.0) SUM2=SUM2+10 . ** ( S2 ( I ) /10 • ) 

IF { S3 ( I ) «GT.0«0) SUM3=SUM3+10 • ** ( S3 ( I ) /10 • ) 

SCBB ( I » J ) =S1 ( I ) 

SCOT ( I » J) =S2 ( I ) 

' SCCTN ( I r J ) =S3 ( I ) 

IF(S1(I). GT. 0,0. OR. S2(I). GT.0.0. OR. S3(I). GT.0.0) SCMX(I»J)= 

2 10.*ALOG10 (10.** (SI ( I )/10. )+10.**(S2 ( I )/10. ) +10 . ** ( S3 ( I ) /10 . ) ) 
62 CONTINUE 

IF (SUM1. GT.0.0) OSCBB( J ) =1 0 . *AL0G10 (SUM1) 

IF (SUM2. GT.0.0) OSCDT ( J ) =10 . *AL0G1 0 ( SUM2 ) 

IF (SUM3. GT.0.0) OSCCTN ( J ) =10 . *ALOGlO ( SUM3 ) 

SUMC=SUM1+SUM2+SUM3 

IF (SUMC. GT.0.0) OASCMX ( J ) =10 . *ALOG10 (SUMC) 

70 CONTINUE 

COMPUTE TURBINE NOISE 

DO 71 1=1,24 
ST ( I » J ) =0 • 0 

71 CONTINUE 
OASTUR ( J ) =0 . 0 
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401 IF { SRC ( 4 ) . EQ . 0 ) GO TO 80 

402 XNB=NTB 

403 CALL TURB(F,PRr TRPM»DTT»AP»XNBrXX,YYfRPEFTC>TC4>FDT) 

404 IF(ERR.EQ.l) GO TO 1000 

405 SUM=0.0 

406 DO 72 1=1 » 24 

407 STd»J)=STCd) 

408 IF(STC(I) .GT.0.0) SUM=SUM+10 .** (STC ( I ) /10 . ) 

409 72 CONTINUE 

410 IF(SUM, GT.0.0) OASTUR ( J ) =10 ,*ALOG10 ( SUM ) 

411 80 CONTINUE 

412 C 

413 C COMPUTE CORE NOISE 

414 C 

415 DO 81 1=1 f 24 

416 SCR ( I r J ) =0 , 0 

417 81 CONTINUE 

418 OASCOR ( J ) =0 . 0 

419 C 

420 IF ( SRC ( 5 ) • EQ • 0 ) GO TO 90 

421 CALL C0RE(FrW3»RH03»T3»T4>DTTD»XXrYY»NC»RREFTC>TC5rFDCR) 

422 IF ( ERR . EQ . 1 ) GO TO 1000 

423 SUM=0 . 0 

424 DO 82 1=1 r24 

425 SCR ( I r J ) =STC ( I ) 

426 IF(STC(I) .GT.0.0) SUM=SUM+10 . ** ( STC ( I ) /10 . ) 

427 82 CONTINUE 

428 IF(SUM, GT.0.0) OASCOR (J) =10. *AL0G10 (SUM) 

429 90 CONTINUE 

430 C 

431 C COMPUTE JET MIXING NOISE 

432 C 

433 DO 91 1=1 >24 

434 SJP(I»J)=0.0 

435 SJS(I»J)=0.0 

436 S JMX ( I » J ) =0 * 0 

437 91 CONTINUE 

438 OSJP(J)=0.0 

439 OSJS(J)=0.0 

440 OASJMX ( J ) =0 • 0 

441 C 

442 IF ( SRC ( 6 ) . EQ . 0 ) GO TO 200 

443 XJ=-XX 

444 Y J=YY 

445 NCOAX=NC— 1 

446 CALL JETMX (F»TA>D1 » VI »Tlf NPLUGrNCOAXr XJr YJ) 

4 U 7 IF(ERR.EQ.l) GO TO 1000 

448 SUM1=0 • 0 

449 SUM2=0 • 0 

450 DO 93 1=1 » 24 
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456 

457 
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459 

460 

461 

462 

463 C 

464 C 

465 C 

466 

467 
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469 
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471 

472 

473 

474 

475 

476 

477 
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486 

487 

488 

489 
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493 C 

494 C 

495 C 

496 

497 

498 

499 

500 


IF(SJ1(I) .GT.0.0) SUM1=SUM1+10.**(SJ1(I)/10.) 

IF ( SJ2 ( I ) • GT. 0 « 0 ) SUM2=SUM2+10.**(SJ2(I)/10. ) 
SJP(Ir J)=SJ1(I> 

SJS ( I > J ) =S J2 ( I ) 

IF(SJ1(I) .GT.0.0. 0R.SJ2<I) .GT.0.0) SJMX(IrJ)= 

2 10.*ALOG10 (10.**(SJ1 ( I )/10. ) +10 . ** ( S J2 ( I ) /10 . ) ) 
93 CONTINUE 

IF (SUM1. GT.0.0) OSJPt J ) =10 . *ALOG10 ( SUM1 ) 

IF (SUM2. GT.0.0) 0SJS(J)=1Q.*AL0G10( SUM2 ) 
SUMT=SUM1+SUM2 

IF (5UMT. GT.0.0) 0A5 JMX ( J ) =10 . *ALOG10 ( 5UMT ) 

200 CONTINUE 

COMPUTE SHOCK BROADBAND NOISE 

DO 250 1=1 f 24 
SPSBB ( I * J ) =0 . 0 
SSSBB ( I f U) =0 • 0 
250 CONTINUE 

OSPSBB ( J ) =0 . 0 
OSSSBB ( J ) =0 • 0 
IF ( SRC ( 7 ) . EQ . 0 ) GO TO 255 

CALL SHOCK(F,VlrTlfDElfDElr7fRG,GCfXXfYY,SBBl) 

IF (ERR.EQ.l ) GO TO 1000 

SUM=0.0 

DO 253 1=1 »24 

IF(SBB1(I) .GT.0.0) SUM=SUM+10 • ** (SBB1(I)/10, ) 
SPSBB ( I r J ) =SBB1 ( I ) 

253 CONTINUE 

IF ( SUM , GT . 0 , 0 ) OSPSBB ( J ) =10 . *ALOG10 (SUM) 

255 CONTINUE 

IF ( SRC ( 8 ) .EQ . 0 ) GO TO 260 

CALL SHOCK (Ff V2 r T2 fDE 2* HA»8r RGrGCf XXrYYrSB62) 

IF ( ERR , EQ . 1 ) GO TO 1000 
SUM=0 . 0 
DO 258 1=1 r 24 

IF ( SBB2 ( I ) .GT.0.0) SUM=SUM+10 . ** ( SBB2 ( I ) /1 0 . ) 
SSSBB( I , J)=SBB2( I ) 

258 CONTINUE 

IF (SUM. GT.0.0) OSSSBB ( J ) =1 0 . *AL0G1 0 (SUM) 

260 CONTINUE 

COMPUTE SHOCK SCREECH NOISE 

DO 262 1=1 r 24 
SSP(I>J)=0.0 
262 CONTINUE 

OSSSP( J)=0.0 

IF ( SRC ( 9 ) • EQ , 0 ) GO TO 265 
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529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 C 

541 

542 

543 C 

544 C 

545 C 

546 
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IF ( NPR1 . LT , 2 • 0 ) GO TO 265 

CALL SCRCH<V1>DE1>DE1>NPR1>A1> AREFSSr XX >YY> 9 > SS1>SS2 >NB1>NB2) 
IF(ERR.EQ.l) GO TO 1000 
DO 264 1=1 >24 
IF(I.EQ.NBl) SSP(I>J)=SS1 
IF { I , EQ.NB2 ) SSP ( I > J ) =SS2 

264 CONTINUE 

IF(SS1.GT.O.O.OR.SS2.GT.O.O) OSSSP(J)= 

2 10.*AL0610 (10. **(551/10. ) +10.** (SS2/10. ) ) 

265 DO 267 1=1 >24 
SSS(I»J)=0.0 

267 CONTINUE 
OSSSS ( J)=0.0 

IF ( SRC ( 10 ) . EQ. 0 ) GO TO 270 
IF( NPR2.LT, 2.0) GO TO 270 

CALL SCRCH(V2>HA> DE2 > NPR2 >A2>AREFSS>XX>YY>10>SS1>SS2>NB1> NB2 ) 

IF ( ERR , EQ. 1 ) GO TO 1000 
DO 268 1=1 >24 
IF(I.EQ.NBl) S5S ( I > J ) =SS1 
IF ( I .EQ.NB2) SSS ( I > J ) =SS2 

268 CONTINUE 

IF(SS1,GT,0. 0 .0R.SS2. GT .0,0) OSSSS (J)= 

2 10.*ALOG10(10.**(SS1/10. ) +10 . ** ( SS2/10 . ) ) 

270 CONTINUE 

DO 598 1=1 >24 
STOT ( I > J ) =0 • 0 
598 CONTINUE 

OASTOT ( J ) =0 . 0 
SUMT=0 , 0 
DO 600 1=1 >24 

S=10.**(SFAN(I>J)/10,)+10.**(SCMX(I>J)/10, ) +10 . ** ( ST ( I > J ) /I 0 . ) + 

2 10.**(SCR(I> J)/10. )+10.**{SJMX(I> J)/10, ) +10 . ** ( SPSBB ( I > J ) /10 . ) + 

3 10 , ** (SSSBB(I>J)/10,)+10.**(SSP(I>J)/10, ) +10 . ** { SSS ( I > J ) /10 . ) 
IF(S.GT.O.O) STOT ( I > J ) =10 , *ALOG10 ( S ) 

SUMT=SUMT+S 
600 CONTINUE 

IF (SUMT.GT.0.0) OASTOT ( J ) =10 , *ALOG10 ( SUMT ) 

J=J+1 

IF ( J.LE.NENG) GO TO 45 

COMPUTE AIRFRAME NOISE 

DO 272 1=1 >24 
STEM { I ) =0 . 0 
STBL ( I ) =0 , 0 
272 CONTINUE 
OASTEN=0 , 0 
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551 OASTBL=0.0 

552 C 

553 IF ( SRC ( 11 ) .EQ.O) GO TO 280 

554 CALL TEN ( F t XLTE » DELTA f CBAR r MUA r VREF > XTE » YTE * STEN ) 

555 IF (ERR.EQ.l ) GO TO 1000 

556 SUM=0.0 

557 DO 274 I=l»24 

558 IF(STENd) .GT.O.O) SUM=SUM+10 . ** ( STEN ( I ) /10 . ) 

559 274 CONTINUE 

560 IF(SUM, GT.O.O) OASTEN=10.*ALOG10 (SUM) 

561 280 CONTINUE 

562 C 

563 IF(SRC(12) .EQ.O) GO TO 290 

564 CALL TBLNCXTrYTrZTr AT r LT r DELS f N J , MUA r DELSC ) 

565 IF (ERR.EQ.l ) GO TO 1000 

566 SUM=0.0 

567 DO 282 1=1 r24 

568 IF(STBL(I) .GT.O.O) SUM=SUM+10 . ** ( STBL (I ) /10 . ) 

569 282 CONTINUE 

570 IF<SUM. GT.O.O') OASTBL=10 . *ALOG10 ( SUM ) 

571 290 CONTINUE 

572 C 

573 IF(POl.EQ.O) GO TO 210 

574 IFCP01.EQ.1) GO TO 203 

575 DO 202 J=1 t NENG 

576 WRITE ( 6 1 305 ) TITLE 

577 WRITE (6r300 ) ALT t MA r P5C t LOC ID t NOENG ( J ) 

578 WRITE (6»301 ) 

579 WRITE ( 6 > 302 ) ( F (I ) , SFF (I t J ) , SFA (I , J ) , SPAN (I , J ) , SCMX (I , J ) t 

580 2 ST(I»J) rSCR(I» J) rSJP(IrJ) »5JS(I> J) »SJMXd» J) »SPSBBd»J) » 

581 3 SSSBB(I»J) »SSP(I»J) fSSS(IrJ) rSTOTdi J) »I=1»24) 

582 WR ITE ( 6 f 304 ) OASFF ( J ) ,OASFA(J) 1 0A5FAN ( J ) >OASCMX(J) rOASTUR(J ) t 

583 2 OASCOR(J) rOSJP(J) t OS JS (J) »OASJMX(J) »05PSBB(J ) t 

584 3 0S55BB ( J ) 1 0S5SP ( J ) r 05SSS ( J ) » 0A5T0T ( J ) 

585 C 

586 300 FORMAT (T30 r25HFLIGHT CONDITION: ALT =»F7.0 f 

587 2 7H MA =rF6.3>20H POWER SET. CODE = t F5.2 » // t TUO , 

580 3 23H0BSERVER LOCATION NO, = t 14 r // » TU5 t 10HENG INE N0.»I4»//) 

589 301 FORMAT (T10r51Hl/3 OCT. BAND SPL*S AND OASPL’S DB RE 0.00002 N/M2 » 

590 2 //»T3f4HFREQrTllr 17HF AN FAN FAN » T54 , 16HJET MIXING NOISE > 

591 3 5X r 25HSH0CK BB SHOCK SCREECH ,/ f T3 » 4H ( HZ )» Til , 10HFWD AFT t 

592 4 50H TOTAL CMX TURB CORE PRI SEC TOT AL » 3X r 

593 5 39HPRI SEC PRI SEC TOTAL >/ ) 

594 302 FORMAT (1XfF 6.0» 13F7 . 1 f 7X r F7 . 1 ) 

595 304 FORMAT ( / » 2X p * OASPL * rl3F7.1*7XrF7.1) 

596 305 FORMAT (lHl»//t8X»12A6r/»8X» 12A6 1 ///) 

597 202 CONTINUE 

598 C 

599 203 CONTINUE 

600 DO 518 1=1 » 24 


54 



601 

5T0T AL ( I ) =0 • 0 

602 

SFFT ( I ) =0 « 0 

603 

SFAT ( I ) =0 • 0 

604 

SFANT { I ) =0 • 0 

605 

SCMXT ( I ) =0 • 0 

606 

STT(I)=0.0 

607 

SCRT ( I ) =0 • 0 

608 

S JPT ( I ) =0 . 0 

609 

S J5T ( I ) =0 • 0 

610 

S JMXT ( I ) =0 « 0 

611 

SPSB8T ( I ) =0 . 0 

612 

SSSBBT ( I ) =0 • 0 

613 

SSPT ( I ) =0 • 0 

614 

5S5T ( I ) =0 • 0 

615 

518 CONTINUE 

616 

OSUM1=0.0 

617 

OSUM2=0.0 

618 

0SUMF=0.0 

619 

OSUM3=0.0 

620 

OSUM4=0.0 

621 

OSUM5=0.0 

622 

OSUM6P=0.0 

623 

OSUM6S=0.0 

624 

OSUM6=0.0 

625 

OSUM7=0.0 

626 

OSUM8=0.0 

627 

OSUM9=0.0 

628 

OSUM10=0.0 

629 

OSUMI=0.0 

630 

OAS1=0.0 

631 

OAS2=0.0 

632 

OASF=0.0 

633 

OAS3=0.0 

634 

OAS4=0.0 

635 

OA55=0.0 

636 

OAS6P=0.0 

637 

OAS6S=0.0 

638 

OAS6=0.0 

639 

OA57=0.0 

640 

OAS8=0.0 

641 

OAS9=0,0 

642 

OA510=0.0 

643 

OAST =0.0 

644 

DO 520 1=1 f 24 

645 

SUM1=0 * 0 

646 

SUM2=0.0 

647 

SUM3=0 • 0 

648 

SUM4=0 • 0 

649 

SUM5=0.0 

650 

SUM6P=0 • 0 
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SUM6S=0.0 

SUM7=0.0 

SUM8=0.0 

SUM9=0.0 

SUM10=0.0 

DO 510 J=1 r NENG 

IF(SFFdrJ) .GT.O.O) SUM1=SUM1+10 . **(SFF(IfJ)/10.) 

IF(SFA(If J) .GT.O.O) SUM2=SUM2+10.**(SFA dr JJ/10. ) 

IF(SCMXdrJ) .GT.O.O) SUM3=SUM3+10.**(SCMX(Ir J)/10. ) 

IF (ST (I r J ) .GT.O.O) SUM4=SUM4+10 . ** (ST(IrJ)/10.) 

IF (SCR (I r J) .GT.O.O) SUM5=SUM5+10 • ** (SCR(IrJ)/10.) 

IF <5 JR (If J) .GT.O.O) SUM6P=SUM6P+10.**(SJP(If U)/10. ) 

IF(SJSdfJ) .GT.O.O) SUM6S=SUM6S+10.**(SJSdr J)/10. ) 

IF ( SP5BB (I r J ) .GT.O.O) SUM7=SUM7+10 . ** ( SPSBB (I t J ) /10 . ) 
IF(SSSBBdrJ) .GT.O.O) SUM8=SUM8+10 . ** ( SSSBB (I r J ) /10 . ) 
IF(SSPdfJ) .GT.O.O) SUM9=SUM9+10.**(SSP(If J)/10. ) 

IF (SSS ( I f J) .GT.O.O) SUM10=SUM10+10.**(SSS(If J)/10. ) 

510 CONTINUE 

SUMI=SUM1+SUM2+SUM3+SUM4+SUM5+SUM6P+SUM6S+SUM7+SUM8+SUM9+SUM10+ 
2 10 • ** ( STEN (I ) /10 . ) +10 • **( STBL ( I ) /10 . ) 

IF (SUMI. GT.O.O) STOTALd )=10.*ALOG10 (SUM I ) 

IF (SUM1. GT.O.O) SFFTd ) =10 . *ALOG10 ( SUMI ) 

IF (SUM2. GT.O.O) SFAT ( I ) =10 . *ALOG10 ( SUM2 ) 

SUMF=SUM1+SUM2 

IF (5UMF. GT.O.O) SFANT (I ) =10 . *ALOG10 ( SUMF ) 

IF (SUM3. GT.O.O) SCMXT ( I ) =10 .*ALOG10 (SUM3 ) 

I F ( SUM4 .GT.O.O) 5TT (I ) =10 . *AL0G1 0 ( SUM4 ) 

I F ( SUM5 .GT.O.O) SCRT (I ) =1 0 . *ALOG10 ( SUM5 ) 

IF (5UM6P. GT.O.O) 5 JPT (I ) =10 . *ALOG10 ( SUM6P ) 

IF (5UM6S. GT.O.O) S JST (I ) =10 . *ALOG10 ( 5UM6S ) 

SUM6=SUM6P+SUM6S 

IF (SUM6. GT.O.O) S JMXT (I ) =10 . *ALOG10 ( 5UM6 ) 

IF (SUM7. GT.O.O) SP5BBT ( I ) =10 . *ALOG10 ( SUM7 ) 

IF (5UM8. GT.O.O) SSSBBT ( I ) =10 . *ALOG10 ( SUM8 ) 

IF (SUM9. GT.O.O) SSPT (I)=10.*ALOG10( SUM9 ) 

IF (SUM10. GT.O.O) SSSTd )=10.*ALOG10 (SUM10) 

OSUMI=OSUMI+SUMI 

0SUM1=05UM1+SUM1 

0SUM2=05UM2+SUM2 

OSUMF=OSUMF+SUMF 

0SUM3=05UM3+5UM3 

0SUM4=05UM4+5UM4 

05UM5=0SUM5+SUM5 

0SUM6P=0SUM6P+SUM6P 

0SUM6S=0SUM6S+SUM6S 

05UM6=0SUM6+SUM6 

05UM7=0SUM7+5UM7 

05UM8=0SUM8+SUM8 
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701 0SUM9=0SUM9+SUM9 

702 05UM10=OSUM10+SUM10 

703 520 CONTINUE 

704 C 

705 

706 

707 

708 

709 

710 

711 

712 

713 

714 

715 

716 

717 

718 

719 C 

720 

721 

722 

723 

724 

725 

726 

727 

728 C 

729 605 FORMAT (T30*25HFLIGHT CONDITION: ALT =*F7.0*7H MA = * 

730 2 F6.3*20H POWER SET. CODE = * F5.2 ► // * T40 * 

731 3 23HOBSERVER LOCATION NO. =*I4*//*T25* 

732 4 31HT0TAL ENGINE AND AIRFRAME NOISE*//) 

733 601 FORMAT (T10*51Hl/3 OCT. BAND SPL’S AND OASPL'S DB RE 0.00002 N/M2 * 

734 2 //*T3*4HFREQ*T11*17HFAN FAN FAN * T54 * 16HJET MIXING NOISE* 

735 3 5X * 25HSH0CK BB SHOCK SCREECH * 3X * 8HAIRFRAME */* T3 * 

736 4 4H(HZ) *T11* 10HFWD AFT* 

737 5 50H TOTAL CMX TURB CORE PR I SEC TOTAL* 

738 6 3X * 24HPR I SEC PR I SEC*4X* 

739 7 19HTEN TBLN -TOTAL*/) 

740 602 FORMAT(1X*F6.0*15F7.1*F8.1) 

741 604 F0RMAT(/*2X* 5HOASPL * 15F7 • 1 * F8 , 1 ) 

742 210 CONTINUE 

743 IF(SRC(l).EQ.O. AND. SRC ( 2 ) . EQ , 0 • OR . P02 • ME • 1 ) GO TO 220 

744 DO 206 J=1*NENG 

745 WR I TE ( 6 * 305 ) TITLE 

746 WRITE (6*300 ) ALT * MA * PSC * LOC ID * NOENG ( J ) 

747 WRITE ( 6 * 306 ) 

748 WRITE (6*308) ( F ( I ) * SFFBB ( I * J ) * SFFDT ( I * J ) * SFFCTN ( I * J ) * SFF ( I * J ) * 

749 2 SFABB(I*J) *SFADT(I*J) *SFA(I*J) *SFAN(I*J) *1=1*24) 

750 WRITE (6*309 ) OSFFBB ( J ) * OSFFDT ( J ) * OSCTN ( J ) * OASFF ( J ) * 


IF ( OSUM1 . GT . 0 . 0 ) OAS1=10 . *ALOG10 ( OSUM1 ) 

IF ( 0SUM2 . GT. 0 . 0 ) OAS2=10 . *ALOG10 ( 0SUM2 ) 

IF ( OSUMF , GT. 0 . 0 ) 0ASF=10 . *ALOG10 < OSUMF ) 

IF ( 0SUM3 • GT • 0 • 0 ) OAS3=10 . *ALOG10 ( 0SUM3 ) 

I F ( OSUM4 . GT . 0 . 0 ) OAS4=10 . *ALOG10 ( 0SUM4 ) 

IF { 0SUM5 . GT , 0 . 0 ) OAS5=10 . *AL0G10 ( 0SUM5 ) 

IF ( 0SUM6P . GT. 0 • 0 ) OAS6P=10 . *ALOG10 ( 0SUM6P ) 

IF ( 0SUM6S • GT . 0 . 0 ) OAS6S=10,*ALOG10 (0SUM6S) 

I F ( 0SUM6 . GT , 0 . 0 ) OAS6=10 .*ALOG10 (0SUM6 ) 

IF ( 0SUM7 • GT, 0 , 0 ) OAS7=10 . *ALOG10 ( 0SUM7 ) 

IF ( 0SUM8, GT ,0.0) OAS8=10 . *ALOG10 ( 0SUM8 ) 

IF ( 0SUM9 , GT, 0 , 0 ) OAS9=10 . *ALOG10 ( 0SUM9 ) 

IF ( OSUMIO , GT. 0 • 0 ) OAS10=10.*ALOG10 (OSUMIO ) 

I F ( OSUM I . GT . 0 . 0 ) OAST=10 , *ALOG10 ( OSUM I ) 

WR I TE ( 6 * 305 ) TITLE 

WR I TE ( 6 * 605 ) ALT * MA * PSC * LOC ID 

WRITE (6*601) 

WRITE (6* 602) ( F ( I ) * SFFT ( I ) * SFAT ( I ) * SFANT ( I ) * SCMXT ( I ) * STT ( I ) * 

2 SCRT ( I ) * S JPT ( I ) * S JST ( I ) * SJMXT ( I ) * SPSBBT ( I ) * SSSB8T ( I ) * 

3 SSPT ( I ) * SSST ( I ) * STEN ( I ) * STBL ( I ) * STOTAL ( I ) *1=1*24) 

WRITE (6*604) OAS1 * OAS2 * OASF * 0AS3 * 0AS4 * 0AS5 * 0AS6P * 0AS6S * 

2 0AS6 * 0AS7 * 0AS8 * 0AS9 * OASIO * OASTEN * OASTBL * OAST 
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751 2 OSFABB ( J ) r OSFADT ( J ) # OASFA ( J ) » OASFAN ( J ) 

752 C 

753 306 FORMAT ( T24 » 20HF AN NOISE COMPONENTS # 

754 2 // # T12 # 20H1/3 OCT. BAND SPL’S » 

755 2 32HAND OASPL’S - DB RE 0.00002 N/M2»//»T40# 

756 3 19H** EXHAUST DUCT **#/#T8# 

757 4 51H*** INLET RADIATED NOISE *** * RADIATED NOISE *#/» 

758 5 T63 r 5HTOTAL#/# T4 #25HFREQ BROAD DSCRT CMBNTN»T40 1 

759 6 12HBR0AD DSCRT , T64 t 3HFAN t /# T4 r 

760 7 55H ( HZ ) -BAND TONES TONES TOTAL -BAND TONES TOTAL r 

761 8 9H NOISE#/) 

762 308 F0RMAT(1X#F6.0»4F7.1#2X#3F7.1#2X#F7.1) 

763 309 FORMAT ( / # 2X # 5HOASPL ► 4F7 .1#2X#3F7.1#2X#F7.1) 

764 C 

765 206 CONTINUE 

766 220 CONTINUE 

767 C 

768 IF (SRC(3) .EQ.O.OR, P03 . NE. 1 ) GO TO 230 

769 DO 228 J=1#NENG 

770 WRITE ( 6 #305 ) TITLE 

771 WRITE ( 6 # 300 ) ALT # MA # PSC » LOC ID # NOENG ( J ) 

772 WRITE ( 6 # 222 ) 

773 WR ITE ( 6 # 224 ) ( F ( I ) »SCBB(I»J) #SCDT(I#J) # SCCTN ( I , J ) , SCMX ( I , J ) , 

774 2 1 = 1 » 24 ) 

775 WRITE (6 r226 ) OSCBB { J ) , OSCDT ( J ) , OSCCTN ( J ) , OASCMX ( J ) 

776 C 

777 222 FORMAT (T13r27HC0MPRESS0R NOISE COMPONENTS t//t Til t 

778 2 31H1/3 OCT. BAND SPL’S AND OASPL » S t / t T17 » 

779 3 18HDB RE 0.00002 N/M2 » // » Til 1 22HFPEQ BROAD DSCRT t 

780 4 9H CMBNTN t 

781 5 / » T1 1 r 22H ( HZ ) -BAND TONES 1 18H TONES TOTAL r / ) 

782 224 FORMAT (8XrF6.0r4F9.1 ) 

783 226 FORMAT (/»9Xf5H0ASPL»4F9.1) 

784 228 CONTINUE 

785 230 CONTINUE 

786 C 

787 C 

788 NO=NO+l 

789 IF (NO.LE.NLOC ) GO TO 40 

790 NFC=NFC+1 

791 IF (NFC.LE. NFLTC ) GO TO 23 

792 GO TO 1000 

793 C 

794 400 WRITE ( 6 ? 402 ) 

795 402 FORMAT ( // 1 5X t 

796 2 59HAIRPLANE MACH NUMBER MUST BE LESS THAN 1.0 - RUN TERMINATED# 

797 3 /) 

798 1000 STOP 

799 END 
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1 SUBROUTINE FAN < X p Y p Wp N p SCI p SC2 p NCODE p TCI p TC2 p TC3 p DELT p MTR p MTRD p 

2 2 MT pXNRpXNSpRSSpBPFpFREF) 

3 DIMENSION ANG<19) p BBDIRI (19) p BBDIRD (19 ) p DTDIR I ( 19 ) p DTDIRD ( 19 ) p 

4 2 CTNDIR ( 19 ) f F ( 27 ) t SCTN12 (24 ) p SCTN14 (24 ) p SCTN18 ( 24 ) p DSPL ( 24 ) 

5 REAL N p MTR p MTRD pMT pMApLBBARpLBBPKpLBBpLTBARpLTPKpLT p 

6 2 LCB12pLCBi4pLCB18pLC12pLC14rLC18pLC12PKp LC14PKp LC18PK 

7 3 p MDFF p NSFF fMIpMDFAp MDC 

8 INTEGER SC p SCI p SC2 p ERR p BN p TCI p TC2 t TC3 

9 COMMON /FAN/IGV p SPLIBB (24 ) pSPLIDT ( 24 ) r SICTN ( 24 ) p SPLDBB (24 ) p 

10 2 SPLDDT (24 ) p DTREF p WREFp RREF p XLFF p HDFF r MDFF p CDFF p FDFFp NSFF p MI p 

11 3 XLFAp HDFAp MDFAp CDFAp FDFAp XLC p HDC pMDC p CDC p FDC 

12 COMMON MA r RHOA t CAp RHOO r CO p CIMPDp DX1 » DX2 1 ERR 

13 C 

14 C ***** DIRECTIVITY DATA ***** 

15 C 

16 DATA ANG/O.OpIO. p20. p30. p40. p50. p60. p 70. p80 . p90. pIOO. pIIO. p 

17 2 120. p130. p140. p150. p160. p170. p180./ 

18 DATA BBDIRI/-2.2p-1.0p0. pO. pO. p-2.0 p-4.5p-7.5p-11 . Op-15. Op 

19 2 —19 • 5 p —25 . 0 p— 30 • 6 p —36 • 3 p**42 • 1 p —47.6 p "53 « 3 p —58 • 8 p —64 « 6/ 

20 DATA BBDIRD/-41.7 p -37 .4 p -33 , 1 p -28 .8 p -24 , 3 p -20 , 1 p-15.8p-11,5p 

21 2 -8. Op- 5. Op-2. 7 p-1. 2p-0. 3p 0. p-2. Op-6 . Op-IO. Op- 15. Op-20.0/ 

22 DATA DTDIRI/-2.9p-1.5pO. pO. pO. p -1 .2 p -3 . 5 p -6 . 8 p -10 . 5 p -14. 5 p 

23 2 -19.0 p-23.3 p-27.8p-32.4p- 36.9 p-41 .5 p- 46. 0 p-50 .4 p-55, 0/ 

24 DATA DTDIRD/— 38.8 p -34 • 8 p-30.8 p -26 .8 p -22 .8 p -18. 9 p - 15 tOp-11 . Op 

25 2 — 8. 0 p— 5. 0 p— 3. 0 p— 1 • 0 p 0 • p 0 » »— 2. 0 p-5.5 p— 9. 0 »— 13. 0 p— 18. 0/ 

26 DATA CTNDIR/— 9.5 p -8.5 p— 7 • 0 p-5. 0 p—2.0 p 0 • p 0 • p-3.5 p—7.5 p— 9. 0 p 

27 2 -9. 5p-10. Op-IO. 5p-11. Op-11. 5p-12. Op-12. 5p-13. Op-13.5/ 

28 C 

29 C ***** 27 1/3 OB CENTER FREQUENCIES 

30 C 

31 DATA F/50 • p 63 • p 80 . p 100 . p 125 . p 160 . p 200 . p 250 . p 315 . p 400 . p 500 . p 630 . p 

32 2 800. p 1000. p 1250. p 1600. p 2000. p 2500. p 3150. p 4000. p 5000. p 6300. p 

33 3 8000. pIOOOO. p12500. P16000. p20000./ 

34 C 

35 C ***** NOTES*************** 

36 C 

37 C NCODE = 01 FOR FIRST-STAGE FAN 

38 C =02 FOR SECOND-STAGE FAN 

39 C =03 FOR COMPRESSOR WITHOUT CTN COMPONENT 

40 C =04 FOR COMPRESSOR WITH CTN COMPONENT 

41 C 

42 C IGV = 00 NO INLET GUIDE VANES FOR 1ST STAGE FAN 

43 C =01 INLET GUIDE VANES PRESENT FOR 1ST STAGE FAN 

44 C 

45 C IGV * S ASSUMED PRESENT FOR 2ND STAGE FANS AND COMPRESSORS 

46 C 

47 C RUN WILL TERMINATE IF THE BLADE PASSAGE FREQUENCY ( BPF ) 

48 C WORKS OUT TO BE LESS THAN 44.7 OR GREATER THAN 22390 HZ. 

49 C 

50 C BPF MUST BE EQUAL TO OR GREATER THAN 355 HZ OR 
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51 C 

52 C 

53 C 

54 C 

55 C 

56 

57 C 

58 C 

59 C 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 C 

73 

74 

75 

76 

77 C 

78 C 

79 C 

80 C 

81 C 

82 C 

83 C 

84 C 

85 C 

86 C 

87 C 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 
9° 

100 C 


GARBAGE MAY RESULT (FOR COMBINATION TONE NOISE) 

THE REF FREQUENCY FOR BROADBAND SPECTRA (2.5*BPF) 

IS ALLOWED TO BE LARGER THAN 22390 HZ. 

IF (MTR.GT .MTRD) GO TO 820 

***** COMPUTE PARAMETERS COMMON TO INLET AND DISCHARGE DUCT ***** 

DENOM=20.*ALOG10 ( DELT/DTREF ) +10 . *ALOG10 ( W/WREF ) +63 . 0 
RSSCBB=-5.0*ALOG10 (RSS/3Q0. ) ‘ 

RSSCDT=-10.0*ALOG10 (RSS/300. ) 

DELTA=ABS ( MT/ ( 1 , 0- ( XNS/XNR ) ) ) • 

BPF=N*XNR/60. 

IF(BPF.LT .44.7.0R.BPF.GT .22390. ) GO TO 810 
CALL BNDN ( BPF f 27 r BN ) 

FB=F ( BN ) 

IREF=BN 

FREF=2.5*FB 

NFREF=8N+4 

IF (NFREF.LE.27) FREF=F ( NFREF ) 

IF (NC0DE.GE.3) GO TO 8 
IF (SCI . EQ, 1 ) GO TO 8 
GO TO 95 
8 CONTINUE 


* COMPUTE NOISE RADIATED FROM THE INLET * 


***** COMPUTE BROADBAND COMPONENT OF INLET RADIATED NOISE ***** 


***** CALCULATE THE PEAK 1/3 OB SPL ( LBBPK ) 

IF ( MTR , GT, 0 . 9 ) GO TO 14 
IF { MTRD , GT , 1 . 0 ) GO TO 16 
LBBAR=58.5 
GO TO 20 

16 LBBAR=58. 5+20. *ALOG10 (MTRD) 

GO TO 20 

14 IF (MTRD , GE. 1 , 0 ) GO TO 18 
LBBAR=57. 6-20. *ALOG10 (MTR) 

GO TO 20 

18 LBBA"R=57. 6+20. *AL0G10 ( MTRD ) -20 , *ALOG10 (MTR) 
20 CONTINUE 

LBBPK=LBBAR+DENOM 
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101 C 

102 C 

103 C 

104 C 

105 C 

106 

107 

108 

109 

110 
111 
112 

113 

114 C 

115 C 

116 C 

117 

118 C 

119 C 

120 C 

121 C 

122 

123 

124 C 

125 C 

126 C 

127 

128 C 

129 C 

130 C 

131 

132 

133 

134 

135 

136 C 

137 C 

138 C 

139 C 

140 C 

141 C 

142 

143 

144 

145 

146 

147 

148 

149 

150 


***** COMPUTE SOURCE EMISSION ANGLE (PHP) AND ACTUAL SOUND 
***** PROPAGATION PATH LENGTH (RP) FROM INPUT VALUES OF 
***** DESIRED ANGLE AND DISTANCE RELATIVE TO THE INLET. 
***** * {FORWARD SPEED EFFECT) 

SC=1 

IF (NC0DE.GE.3) SC=3 

CALL XFORM (X*Y*0.0*MA*SC*DX1* DX2 * PH * PHP r R * RP * ERR ) 

IF ( ERR ,EQ. 1 ) GO TO 900 

PHI=PH 

RI=R 

PHPIrPHP 

RPI=RP 

***** COMPUTE CORRECTION FOR CONVECTIVE AMPLIFICATION 

C A I =-40 . *AL0G10 ( 1 . 0-MA*COS ( PHP/57 . 29577951 ) ) 

***** CORRECT THE PEAK LEVEL FOR DIRECTIVITY * DISTANCE* 
***** LOCAL IMPEDANCE* AND CONVECTIVE AMPLIFICATION 

DEL=GIRC (PHP * ANG* BBDIRI *19*1) 

LBB=LBBPK+DEL-20.*ALOG10 (RP/RREF) +CIMPD+CAI 

***** CORRECT FOR ROTOR-STATOR SPACING 

LBB=LBB+RSSCBB 

***** COMPUTE BROADBAND SPECTRUM 
DO 40 J=1 * 24 

DEL=10.*ALOG10 (EXP (-4.264* (AL0G10 (F( J)/FREF) )**2) ) 

SPLIBB ( J) =LBB+DEL 

IF ( SPLI BB ( J ) . LT .0.0) SPLIBB ( J ) =0 . 0 
40 CONTINUE 

***** COMPUTE DISCRETE TONE SPL’S OF INLET RADIATED NOISE 


***** COMPUTE PEAK LEVEL OF THE FUNDAMENTAL TONE ( LTPK ) 

IF ( MTR , GT, 0 • 72 ) GO TO 48 
IF ( MTRD, GT ,1.0) GO TO 42 
LTBAR=60 • 5 
GO TO 50 

42 LTBAR=60. 5+20. *ALOG10 (MTRD) 

GO TO 50 

48 IF ( MTRD, GE .1.0) GO TO 46 

XLTBAR=60 , 5+50 • *ALOG10 ( MTR/0 , 72 ) 
YLTBAR=59,5-80,*ALOG10 (MTR) 


***** 
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151 

152 

153 

154 

155 

156 

157 

158 C 

159 C 

160 C 

161 C 

162 

163 

164 C 

165 C 

166 C 

167 

168 C 

169 C 

170 C 

171 C 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 

182 C 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 C 

196 C 

197 C 

198 C 

199 

200 


LTBAR=AMIN1 ( XLTBAR ? YLTBAR ) 

GO TO 50 

46 XLTBAR=60.5+20.*AL0G10 ( MTRD ) +50 . *ALOG10 (MTR/0.72) 
YLTBAR=59.5+80 «*ALOG10 < MTRD/MTR ) 

LTBAR=AMIN1 (XLTBAR? YLTBAR) 

50 CONTINUE 
LTPK=LTBAR+DENOM 

***** CORRECT THE PEAK LEVEL FOR DIRECTIVITY ? DISTANCE? 

***** LOCAL IMPEDANCE? AND CONVECTIVE AMPLIFICATION 

DEL=GIRC ( PHP ? ANG ? DTD IR I ? 19 ? 1 ) 

LT=LTPK+DEL-20.*ALOG10 ( RP/RREF ) +C IMPD+CAI 

***** CORRECT FOR ROTOR-STATOR SPACING 

LT=LT+RSSCDT 

***** CORRECT FOR ROTOR-STATOR INTERACTION AND COMPUTE 
***** THE DISCRETE TONE HARMONIC LEVELS 

DO 51 J=l?24 
SPLIDT ( J)=0«0 

51 CONTINUE 

DO 59 K=1 ? 6 
FBK=FB*K 

IF(FBK. GT. 11220. ) GO TO 60 
CALL BNDN ( FBK ?24 ? BN ) 

IF(K.EQ.l) GO TO 53 
IF(BN.EQ.JK) GO TO 60 
53 JK=BN 

IF ( NCODE. EQ. 1 . AND. I GV« EQ. 0 ) GO TO 56 
IF(K.GT.l) GO TO 55 
SPLIDT ( JK ) =LT 

I F ( DELTA . LT . 1 . 05 ) SPL I DT ( JK ) =LT-8 . 0 
GO TO 58 

55 SPLIDT ( JK ) =LT-3 • 0-3 • 0*K 
GO TO 58 

56 ' SPLIDT ( JK ) =LT+3 • 0-3. 0*K 

IF (DELTA.LT. 1.05. AND. K.EQ.l) SPLIDT ( JK ) =LT-8. 0 

58 CONTINUE 

59 CONTINUE 

60 CONTINUE 

***** COMPUTE COMBINATION TONE NOISE (INLET ONLY) ***** 

(FIRST-STAGE FANS OR COMPRESSORS IF DESIRED WHERE MTR > 1.0) 

DO 71 I=l?24 
SICTN ( I ) =0 . 0 
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201 71 CONTINUE 

202 IF(NCODE.EQ.2.OR.NCODE.EQ.3.OR.MTR.LE.1.0) GO TO 75 

203 C 

204 C ***** COMPUTE PEAK LEVEL FOR EACH OF 3 COMPONENTS 

205 C 

206 70 IF(MTR.GT. 1.146) GO TO 61 

207 LCB12=785.68*ALOG10 (MTR)+30.0 

208 GO TO 62 

209 61 LCB12=-49.62*ALOG10 <MTR) +79.44 

210 62 CONTINUE 

211 IFtMTR.GT. 1.322) GO TO 63 

212 LCB14=391 . 81*ALOG10 ( MTR ) +30 . 0 

213 GO TO 64 

214 63 LCB14=-50 • Q6*ALOG10 ( MTR ) +83. 57 

215 64 CONTINUE 

216 IF(MTR.GT, 1.610) GO TO 65 

217 LCB18=199.20*ALOG10< MTR) +30.0 

218 GO TO 66 

219 65 LCB18=-49.89*AL0G10<MTR)+81.52 

220 66 CONTINUE 

221 LC12PK=LCB12+DEN0M 

222 LC14PK=LCB14+DEN0M 

223 LC18PK=LCB18+DEN0M 

224 C 

225 C ***** CORRECT THE PEAK LEVELS FOR DIRECTIVITY r DISTANCE? 

226 C ***** LOCAL IMPEDANCE? AND CONVECTIVE AMPLIFICATION 

227 C 

228 DEL=G IRC ( PHP ? ANG ? CTNDIR ? 19 ? 1 ) 

229 LC12=LC12PK+DEL-20 , *AL0G10 (RP/RREF) +CIMPD+CAI 

230 LC14=LC14PK+DEL-20,*ALOG10 ( RP/RREF ) +C IMPD+CAI 

231 LC18=LC18PK+DEL-20 . *AL0G10 ( RP/RREF )+C IMPD+CAI 

232 C 

233 C ***** CORRECT FOR IGV*S? IF PRESENT 

234 C 

235 IF (NCODE.EQ.l.AND. IGV.EQ.O) GO TO 72 

236 LC12=LC12-5.0 

237 LC14=LC14-5.0 

238 LC18=LC18-5.0 

239 72 CONTINUE 

240 C 

241 C ***** COMPUTE SPECTRUM SHAPE FOR EACH OF 3 CTN COMPONENTS 

242 C 

243 IR12=IREF-3 

244 IR14=IREF-6 

245 IR18=IREF-9 

246 DO 80 J=l?24 

247 IF (U.LT, IR12) SPLDEL=30 , 0*ALOG1 0 ( F ( J ) /FB ) +9 . 03 

248 IF( J.EQ. IR12) SPLDEL=0.0 

249 IF( J.GT.IR12) SPLDEL=-30 • 0*AL0G10(F(J)/FB) -9 ,03 

250 SCTN12 ( U ) =LC12+SPLDEL 
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251 

252 C 

253 

254 

255 

256 

257 

258 C 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 C 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 C 

296 

297 C 

298 

299 C 

300 C 


IF(SCTN12( J) .LT.0.0) SCTN12 ( J ) =0 . 0 

IF ( J.LT.IR14) SPLDEL=50.0*ALOG10 (F( J)/FB)+30.1 
IF ( J.EQ.IR14) SPLDEL=0.0 

IF ( J.GT.IR14) SPLDEL=-50.0*AL0G10(F( J)/FB)-30.1 
SCTN14 ( J ) =LC14+SPLDEL 
IF <SCTN14(J) .LT.0.0) SCTN14( J)=0.0 

IF< J.LT.IR18) SPLDEL=50.0*ALOG10(F( J)/FB) +45.15 
IF(J.EQ.IRIS) SPLDEL=0.0 

IF { J.GT.IR18) SPLDEL=-30.0*ALOG10 ( F ( J)/FB)-27.09 

5CTN18 ( J ) =LC18+SPLDEL 

IF (SCTN18 ( J) .LT.0.0) SCTN18 ( J ) =0 . 0 

IF(SCTN12(J).GT.O. 0 .OR .SCTN14 ( J ) . GT . 0 . 0 . OR . SCTN18 ( J ) .GT.O.O) 

2 SICTN(J)=10.*ALOG10{10.**(SCTN12(J)/10. ) +10.** (SCTN14 ( J)/10. )+ 
2 10.**(SCTN18(J)/10. ) ) 

80 CONTINUE 
75 CONTINUE 

IF (NC0DE.LE.2) 60 TO 87 
IF(TC3.EQ.1.0R.TC3.EQ,3) GO TO 81 
IF ( TC3 . EQ .2 ) GO TO 86 
GO TO 89 

87 IF(TCl.GE.l) GO TO 86 
GO TO 89 

81 CALL DCTRT(DSPLr XLC > HDC > MDC r CDC > FDC > 0 . 0 r 0 . 0 » PHP , 1 ) 

DO 84 1=1 r 24 

SPLI BB ( I ) =SPL IBB ( I ) -DSPL ( I ) 

SPLIDT ( I ) =SPLIDT ( I ) -DSPL ( I ) 

SICTN ( I ) =SICTN ( I ) —DSPL ( I ) 

84 CONTINUE 

IF ( TC3, EQ. 1 ) GO TO 82 

86 CALL DCTRT ( DSPL r XLFF r HDFF . MDFF t CDFF rFDFF?NSFF»MI t PHP » TCI ) 

DO 83 1=1 f 24 

SPLIBB ( I ) =SPLIBB ( I ) -DSPL ( I ) 

SPLIDT ( I ) =SPLIDT ( I ) -DSPL ( I ) 

SICTN ( I ) =S I CTN ( I ) -DSPL ( I ) 

83 CONTINUE 

82 DO 85 1=1 r24 

IF ( SPLIBB ( I) .LT.0.0) SPLIBB ( I ) =0 . 0 
IF ( SPLIDT ( I ) .LT.0.0) SPLIDT ( I ) =0 , 0 
IF(SICTN(I) .LT.0.0) SICTN(I)=0.0 

85 CONTINUE 
89 CONTINUE 

IF (NCODE.GE.3.0R.SC2.EQ.O) GO TO 195 
95 CONTINUE 
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301 C 

302 C 

303 C 

304 C 

305 C 

306 C 

307 C 

308 C 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 C 

319 C 

320 C 

321 C 

322 C 

323 C 

324 

325 
• 326 

327 

328 

329 

330 C 

331 C 

332 C 

333 

334 C 

335 C 

336 C 

337 C 

338 

339 

340 C 

341 C 

342 C 

343 

344 C 

345 C 

346 C 

347 

348 

349 

350 C 


* COMPUTE NOISE RADIATED FROM THE FAN DISCHARGE DUCT * 
******************************************* ***************** 

***** CALCULATE THE BROADBAND NOISE COMPONENT ***** 


***** COMPUTE THE PEAK 1/3 OB SPL ( XLBBPK ) 

IF(MTR.GT.l.O) GO TO 110 
IF(MTRD.GT.l.O) GO TO 112 
XLBBAR=60.0 
GO TO 116 

112 XLBBAR=60.0+20.*ALOG10 (MTRD) 

GO TO 116 

110 XLBBAR=60. 0+20. *AL0G10 (MTRD) -20. *ALOG10 (MTR) 

116 CONTINUE 

XLBBPK=XLBBAR+DENOM 

***** COMPUTE SOURCE EMISSION ANGLE (PHP) AND ACTUAL SOUND 
***** PROPAGATION PATH LENGTH (RP) FROM INPUT VALUES OF 
***** DESIRED ANGLE AND DISTANCE RELATIVE TO THE DISCHARGE 
***** DUCT. (FORWARD SPEED EFFECT) 

CALL XFORM(Xf Yr 0 . 0 r MA r 2 r DX1 > DX2 > PH » PHP r R , RP f ERR ) 

IF(ERR.EQ.l) GO TO 900 

PHPD=PHP 

RPD=RP 

PHD=PH 

RD=R 


***** COMPUTE CORRECTION FOR CONVECTIVE AMPLIFICATION 

CAD=-40 .*AL0G10 ( 1 . 0-MA*COS ( PHP/57 . 29577951 ) ) 

***** CORRECT THE PEAK LEVEL FOR DIRECTIVITY t DISTANCEr 
***** LOCAL IMPEDANCE f AND CONVECTIVE AMPLIFICATION 

DEL=GIRC(PHPr ANG r BBDIRD f 19 r 1) 

XLBB=XLBBPK+DEL-20.*ALOG10 (RP/RREF ) +CIMPD+CAD 
***** CORRECT FOR ROTOR-STATOR SPACING 
XLBB=XLBB+RSSCBB 

***** CORRECT FOR IGV’Sf IF PRESENT 

IF(NCODE.EQ.l.AND. IGV.EQ.O) GO TO 134 
XLBB=XLBB+3.0 
134 CONTINUE' 
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351 C ***** COMPUTE THE BROADBAND SPECTRUM 

352 C 

353 DO 140 J=1 * 24 

354 DEL=10.*ALOG10 (EXP (-4.264* (ALOGIO (F( J)/FREF) ) **2 ) ) 

355 SPLDBB ( J ) =XLBB+DEL 

356 IF ( SPLDBB ( J ) . LT .0.0) SPLDBB ( J ) =0 . 0 

357 140 CONTINUE 

358 C 

359 C ***** COMPUTE THE DISCRETE TONE LEVELS OF THE DISCHARGE DUCT ***** 

360 C 

361 C 

362 C ***** COMPUTE THE PEAK LEVEL OF THE FUNDAMENTAL TONE ( XLTPK ) 

363 C 

364 IF(MTR.GT.l.O) GO TO 148 

365 IF(MTRD.GT.l.O) GO TO 142 

366 XLTBAR=63.0 

367 GO TO 150 

368 142 XLTBAR=63. 0+20 • * ALOGIO ( MTRD ) 

369 GO TO 150 

370 148 XLTBAR=63 » 0+20 • * ALOGIO ( MTRD ) -20 • * ALOGIO ( MTR ) 

371 150 CONTINUE 

372 XLTPK=XLTBAR+DENOM 

373 C 

374 c ***** CORRECT THE PEAK LEVEL FOR DIRECTIVITY * DISTANCE* 

375 C ***** LOCAL IMPEDANCE* AND CONVECTIVE AMPLIFICATION 

376 C 

377 DEL=GIRC ( PHP * ANG * DTDIRD * 19*1) 

378 XLT=XLTPK+DEL-20 • * ALOGIO ( RP/RREF ) +C IMPD+CAD 

379 C 

380 C ***** CORRECT FOR ROTOR-STATOR SPACING 

381 C 

382 XLT=XLT+RSSCDT 

383 C 

384 C ***** CORRECT FOR IGV’S* IF PRESENT 

385 C 

386 IF ( NCODE • EQ • 1 . AND. IGV « EQ. 0 ) GO TO 154 

387 XLT=XLT+6.0 

388 154 CONTINUE 

389 C 

390 C ***** COMPUTE LEVELS OF THE HARMONICS OF THE FUNDAMENTAL TONE 

391 C 

392 DO 157 J=l*24 

393 SPLDDT ( J ) =0 • 0 

394 157 CONTINUE 

395 DO 165 K=l *6 

396 FBK=FB*K 

397 IFtFBK. GT. 11220. ) GO TO 170 

398 CALL BNDN ( FBK * 24 * BN ) 

399 IF(K.EQ.l) GO TO 159 

400 IFCBN.EQ.JK) GO TO 170 
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401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 C 

415 

416 

417 . 

418 

419 

420 

421 

422 

423 C 

424 

425 C 

426 

427 C 

428 C 

429 C 

430 . 

431 

432 

433 

434 

435 

436 

437 

438 

439 

440 

441 

442 


159 JK=BN 

IF(NCODE.EQ.l.AND. IGV.EQ.O) 60 TO 156 
IF(K.GT.l) GO TO 155 
SPLDDT ( JK ) =XLT 

IF (DELTA.LT. 1.05) SPLDDT ( JK)=XLT-8.0 
GO TO 158 

155 SPLDDT ( JK ) = XLT-3 . 0-3 . 0*K 
GO TO 158 

156 SPLDDT ( JK)=XLT+3 .0-3.0 *K 

IF (DELTA.LT. 1.05. AND. K.EQ.l) SPLDDT ( JK ) =XLT-8 . 0 
158 CONTINUE 
165 CONTINUE 
170 CONTINUE 

IF ( TC2 . EQ. 0 ) GO TO 195 

CALL DCTRT(DSPL,XLFA,HDFA,MDFA, CDFA,FDFA, 0 . 0 , 0 . 0 , PHP , 4 ) 

DO 173 1=1,24 

SPLDBB ( I ) =SPLDBB ( I ) -DSPL ( I ) 

SPLDDT ( I ) =SPLDDT ( I ) -DSPL ( I ) 

IF (SPLDBB ( I ) .LT. 0 . 0 ) SPLDBB ( I ) =0 . 0 
IF (SPLDDT ( I ) «LT. 0 • 0 ) SPLDDT ( I ) =0 . 0 
173 CONTINUE 

195 CONTINUE 

GO TO 900 

***** ERROR MESSAGES ***** 

810 IF (NC0DE.LE.2 ) WRITE(6r320) 

IF ( NCODE. GE • 3 ) WRITE(6,321) 

320 FORMAT (///,2X, ’FAN BPF IS OUTSIDE ALLOWABLE RANGE - RUN*, 

2 ’ TERMINATED*,/) 

321 FORMAT (//,2X,41HC0MPRESS0R BPF IS OUTSIDE ALLOWABLE RANGE, 
2 17H - RUN TERMINATED,/) 

GO TO 890 
820 WR ITE ( 6 , 322 ) 

322 FORMAT (//,2X, »MTR MUST BE LES THAN OR EQUAL TO MTRD * , 

2 * FOR FAN OR COMPRESSOR NOISE - RUN TERMINATED*,/) 

890 ERR=1 
900 RETURN 
END 
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1 

2 

3 

4 

5 

6 

7 

8 C 

9 

10 

11 

12 

13 

14 

15 C 

16 

17 

18 

19 C 

20 C 

21 C 

22 C 

23 

24 C 

25 C 

26 C 

27 C 

28 

29 

30 C 

31 C 

32 C 

33 

34 

35 

36 C 

37 C 

38 C 

39 C 

40 

41 C 

42 C 

43 C 

44 C 

45 C 

46 

47 

48 C 

49 C 

50 C 


SUBROUTINE TURB ( F t PR , N2 t D t AE t XNB t X t Y , RREF t TC4 , FDT ) 

DIMENSION F ( 24 ) >BBSPEC(30) >FBPF(30) 'OASDIR(IO) rSPLDIR(lO ) t 
2 DIRANG ( 10 ) 

REAL N2 p MAfMDTPf DSPL (24 ) 

INTEGER ERRfBN,TC4 

COMMON /TC/SPL (24 ) t AREF t UREF t WREF ► DTREF , XLTP t HDTP , MDTP t CDTP 
COMMON MAfRHOAr CAf RHOO > CO > C IMPD > DX1 1 DX2 f ERR 

DATA BBSPEC/-60.0»-41.0»-36.0'-29.5'-26.0»-21.5»-14.5'-10.9» 

2 -8* 3 1-6 .6 1 -5*4 1 -4 .4>-3.7f-3*0f-2»5r-l«7f- , 1.0r-0«5r -0.25 r 0 • 00 t 

3 — 0.2f— 0.4f— 0.6f— 1«0»-1.3>— 1.7f-2.2f —2 • 8 1 —3 • 5 r —4.5/ 

DATA FBPF/ .001; • 005 r .01; . 02 f . 03 f «05r . 10 r . 15 f . 20 > . 25 f .30 t • 35 1 

2 .40' .45r .50f .60f .70* ,80 t ,90 t 1.00 » 1.10 » 1.20 » 1.30 » 1.40 r 1.50 r 

3 1.60'1.70'1.80rl.90'2.00/ 

DATA 0 AS DIR/— 13 • 8 t —4 . 6 ' —4 .5'— l»5'0.0i>— 1.3r —2 . 8 t —4.2' —5 • 4 ' —6 .6/ 
DATA SPLDIR/-25.6'-9.5'-8.8»-2.5»0.0'-4.3'-6.5'-11.8'-14.4»-19.5/ 
DATA DIRANG/0.0'90. »92.»111. '120. f 130. »140. »150. rl60.rl80./ 


*** COMPUTE DOMINANT (LAST) STAGE TIP SPEED *** 

UT= ( D/2. )* (2.0*3. 14*N2/6Q. ) 

*** COMPUTE TOTAL OASPL (LOSSLESS) AT 120 DEG t 200 FT *** 

*** SIDELINE CORRECTED FOR LOCAL IMPEDANCE *** 

OASPK=40.*ALOG10(1.0-(1.0/PR)**0.286)-20.*ALOG10 (UT/UREF)+ 

2 10,*AL0G10 ( A E/ AREF ) +C IMPD+1 09 . 0 

INITIALIZE SPECTRUM 

DO 10 1=1 '24 
SPL ( I ) =0 • 0 
10 CONTINUE 

*** COMPUTE DISCRETE 1/3 O.B. SPL (LOSSLESS) AT 120 DEG t *** 
*** 200 FT SIDELINE *** 

SPLPK=OASPK— 5.0 

*** COMPUTE SOURCE EMISSION ANGLE AND ACTUAL DISTANCE *** 

*** I.E.r MEASURED ANGLE AND DISTANCE CORRECTED FOR *** 

*** FORVMRD SPEED *** 

CALL XFORM(X»Y»0.0»MA'4rDXlFDX2fPHf PHPrRfRPp ERR) 

IF ( ERR , EQ . 1 ) GO TO 50 

*** APPLY DIRECTIVITY AND DISTANCE CORRECTIONS TO PEAK *** 
*** TOTAL OASPL AND PEAK TONE SPL *** 
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51 C 

52 

53 

54 C 

55 

56 

57 C 

58 C 

59 C 

60 
61 
62 

63 

64 

65 C 

66 C 

67 C 

68 

69 C 

70 C 

71 C 

72 

73 C 

74 C 

75 C 

76 

77 

78 

79 

80 
81 
82 

83 C 

84 C 

85 C 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 C 

96 

97 

98 

99 
100 


OASDEL=GIRC(PHP* DIRANG* OASDIR * 10 * 1) 
OASPL=OASPK+OASDEL-20 • *AL0G1 0 ( RP/RREF ) 

SPLDEL=GIRC(PHP*DIRANG* SPLDIRrlOrl) 
SPLT=SPLPK+SPLDEL-20.*ALOG10 (RP/RREF) 

*** APPLY CORRECTION FOR CONVECTIVE AMPLIFICATION *** 

CAMP=— 40.*ALOG10 ( 1 . 0-MA*COS ( PHP/57 . 3 ) ) 

OASPL=OASPL+CAMP 
5PLT =SPLT +CAMP 

IF(OASPL.LE.O.O.OR.SPLT.LE.O.O) GO TO 50 
IF(SPLT.GE.OASPL) GO TO 48 

*** COMPUTE BROADBAND 0A5PL (LOSSLESS) *** 

OASBB=10.*ALOG10 ( 10 . ** ( OASPL/IO . ) -10 . ** ( SPLT/10 . ) ) 

*** COMPUTE BLADE PASSAGE FREQUENCY *** 

BPF=XNB*N2/60. 

*** FIND 1/3 O.B. CONTAINING THE BPF *** 

IF ( BPF, GE • 44 . 7 • AND , BPF, LE. 11220. ) GO TO 19 

FPK=BPF 

BN=0 

GO TO 21 

19 CALL BNDN { BPF *24 > BN ) 

FPK=F(BN) 

21 CONTINUE 

*** COMPUTE TOTAL SPECTRUM (LOSSLESS) *** 

DO 30 J=l*24 
FRAT=F( J)/FPK 

DEL=GIRC ( FRATr FBPF r BBSPEC » 30 f 1 ) 

SPL ( J ) =DEL-7 • 3+OASBB 

IF(SPL( J) .LT. 0.0. AND. J.NE.BN) SPL (J) =0.0 
IF(J.NE.BN) GO TO 30 
SPLPBB=SPL( J) 

SPL( J) =10, *ALOG10 (10. *+(SPLPBB/10. )+10.**( SPLT/10. ) ) 

30 CONTINUE 

IF ( TC4 . EQ , 0 ) GO TO 50 

CALL DCTRT(DSPL>XLTPr HDTPfMDTPrCDTPrFDTf O.Or 0.0rPHP f 4) 

DO 32 1=1*24 

SPL ( I ) =SPL ( I ) -DSPL ( I ) 

IF ( SPL ( I ) • LT . 0 , 0 ) SPL ( I ) =0 , 0 


) 
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101 32 CONTINUE 

102 .GO TO 50 

103 C 

104 C 

105 48 WR ITE ( 6 r 100 ) 

106 100 FORMAT (//f5X»45HT0NE SPL NOT LESS THAN TOTAL OASPL IN TURBINE t 

107 2 23H NOISE - RUN TERMINATED t / ) 

108 ERR=1 

109 50 RETURN 

110 END 
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1 SUBROUTINE CORE ( FREQ * W3 * RH03* T3 * T4* DT * X* Y* NFLAG*RREF * TC5 * FDCR ) 

2 DIMENSION FREQ ( 24 ) * DELPWL ( 24 ) * DIRDF ( 24 * 19 ) * 

3 2 ANGDF (19) * D IRSF ( 19 ) * ANGSF ( 19 ) * SPWL (24 ) * SPLA ( 24 ) 

4 REAL MA* MDTP * DSPL (24 ) 

5 INTEGER ERR * TC5 

6 COMMON /TC/SPL (24) * AREF*UREF * WREF *DTREF * XLTP * HDTP * MDTP * CDTP 

7 COMMON MA*RHOA*CA*RHOO*CO*CIMPD*DX1* DX2 * ERR 

8 C 

9 C POWER SPECTRUM SHAPE * REF P156* FAA-RD-77-4 

10 C 

11 DATA DELPWL/-24 • 0 * -20 • 0*-16.0*-13»0*-10.0* -7 • 0 * -4 . 5 * -2*5 r 

12 2 -1.0r0.0r-1.0r-2.5»-4.5»-7.0r-10.0»-13.0»-16.0»-20.0»-24.0» 

13 3 -27.5»-31.5f-36.0»-40.0»-45.0/ 

14 C 

15 C DIRECTIVITY INDICES FOR DUAL-FLOW ENGINE EXHAUST 

16 C GEOMETRY* REF P155* FAA-RD-77-4. EXTRAPOLATED BELOW 40 DEG AND 

17 C ABOVE 140 DEG TO COVER RANGE FROM 0 TO 180 DEG. OCT 77. 

18 C 

19 C 

20 DATA ( (DIRDF(I*J} ,1=1*24) * J=1 * 19 ) /4*-5 . 1 * 3*-9 . 4 r 17*-15 . 9 r 

21 2 4*-4 . 8 * 3*-8 • 6 * 17*-14 , 4 * 4*-4 , 5 * 3*-8 ,0*17*-12.9* 

22 3 4*-4 . 3 * 3*-7 • 2 * 17*-11 « 5 * 4*-4, 0 * 3*-6. 5 * 17*-10 • 0 * 

23 2 4*-3.8 * 3*-5.8* 17*-8.5* 4*-3.2 *-5. 0 *-4, 5 *-5. 0 * 17*-6.5 * 4*-3. 0 * 

24 3 -4.0 *-3.5 * 18*-4,5*4*-2.7*-3. 0 *-2 .5 *-4, 0 * 17*-2.5 * 4*-2 . 0 * 

25 4 -1.5*-1.5*-3.0* 17*-0 .5 *4*-0 .8 * 0. 0*0. 0*-l. 8*17*1. 0*4*0. 8* 

26 5 1.8*1.5*1.0*17*2.5*4*3.0*4.0*3.5*3.5*17*5.0*5*5.0*4.8* 

27 6 5.5*17*4.5*4*7.0*5.0*6.0*6.5*17*3.5*4*7.2*4.7*5.9*6.6*17*2.2* 

28 7 4*7 • 0 * 4 • 0 *5 .3 *6 .3 * 17*0 .7 * 4*6 .4 *3 • 1 *4 .4 *5 .6 * 17*— 1 • 1 * 

29 8 4*5.6 *1.8*3, 2*4, 5* 17*-3 ■ 2/ 

30 DATA ANGDF/0. 0*10, *20. *30, *40. *50. *60. *70. *80, *90. *100. *110. * 

31 2 120. *130. *140. *150. *160. *170. *180./ 

32 C 

33 C DIRECTIVITY INDICES FOR SINGLE-FLOW ENGINE EXHAUST GEOMETRY 

34 C REF GE*S PROPOSED APPENDIX TO ARP 876 FOR COMBUSTOR NOISE 

35 C PREDICTION* JULY 1977. EXTRAPOLATED BELOW 10 DEG AND ABOVE 

36 C 160 DEG TO COVER RANGE FROM 0 TO 180 DEG. OCT 77, 

37 C 

38 DATA DIRSF/-8. 5 *-8.0 *-7. 5 *-7.0 *-6.5 *-6.0 *-5.3 *-4.6 *-4.0* -1.7* 

39 2 0 • 7 * 3 . 0 * 5 • 0 * 3 . 5 * 1 • 2 * -1 . 9 * -5 . 1 * -8 . 8 * -12 • 6/ 

40 DATA ANGSF/0. 0*10. *20, *30. *40. *50. *60. *70. *80. *90. *100. *110. * 

41 2 120. *130. *140. *150. *160. *170. *180./ 

42 C 

43 C 

44 C THIS PROGRAM COMPUTES A CORE NOISE PREDICTION AT A POINT USING 

45 C THE GENERAL ELECTRIC ENGINE PREDICTION METHOD DESCRIBED IN 

46 C FAA-RD-77-4* DATED FEB 1977* AND THEIR PROPOSED APPENDIX TO 

47 C ARP 876 (COMBUSTOR NOISE PREDICTION) FOR THE SAE A-21 JET 

48 C NOISE SUBCOMMITTEE* JULY 1977. 

49 C SPL’S ARE DB RE 0.00002 N/M2 FREE-FIELD (LOSSLESS) 

50 C 
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51 C 

52 C 

53 C 

54 

55 

56 C 

57 C 

58 C 

59 C 

60 
61 
62 

63 C 

64 C 

65 C 

66 C 

67 C 

68 

69 

70 

71 C 

72 C 

73 C 

74 C 

75 C 

76 C 

77 C 

78 C 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 C 

94 C 

95 C 

96 C 

97 C 

98 C 

99 C 
100 


COMPUTE OVERALL SOUND POWER 

OAPWL=10.*ALOG10 ( (W3/WREF)*( ( T4-T3 ) /DTREF ) **2* ( RHO3/RHO0 ) **2 ) 
2 —40 , *ALOG10 (DT/DTREF)+123.5 

COMPUTE POWER SPECTRUM REF. P156 OF FAA-RD-77-4 
PEAK FREQUENCY IS 400 HZ. 

SPWLPK=0APWL-6.8 
DO 10 1=1 f 24 

10 SPWL ( I ) =SPWLPK+DELPWL ( I ) 

CONVERT POWER SPECTRUM TO SOUND PRESSURE LEVEL SPECTRUM 
IMPEDANCE COMPUTED FROM AMBIENT STATIC CONDITIONS 
ATMOSPHERIC ATTENUATION IS NEGLECTED 

DO 20 1=1 f 24 

SPLA ( I ) =SPWL ( I ) +CIMPD-20 .8 
20 CONTINUE 


APPLY DIRECTIVITY FOR SINGLE-FLOW OR DUAL-FLOW ENGINE 
EXHAUST TO OBTAIN SPL*S AT DESIRED ANGLE. 

COMPUTE EQUIVALENT STATIC DIRECTIVITY ANGLE AND DISTANCE 
IF MA > 0.0 

ALSO CORRECT FOR SPHERICAL DISPERSION 

CALL XFORM(X>Y>O.OfMA*5fDXlfDX2fPHf PHP ?R r RP r ERR ) 

IF ( ERR . EQ, 1 ) GO TO 70 
IF (NFLAG.EQ.2 ) GO TO 40 
DIRCOR=GIRC( PHP t ANGSF fDIRSF»19>l) 

DO 30 1=1 f 24 

SPL ( I ) =SPLA ( I ) +DIRCOR-20 . *ALOG10 ( RP/RREF ) -37 . 0 
30 CONTINUE 
GO TO 50 
40 DO 45 1=1 f 24 
F=FREQ( I ) 

DC=DT AB2 ( PHP , F , ANGDF , FREQ r 19 1 24 , 1 r 1 , D I RDF 1 24 t I ERR ) 

SPL ( I ) =SPLA ( I ) +DC-20 . *ALOG10 ( RP/RREF ) -37 . 0 
45 CONTINUE ' 

50 CONTINUE 

NO RELATIVE VELOCITY EFFECT IS APPLIED TO POWER LEVEL 
NO DYNAMIC AMPLIFICATION EFFECT IS APPLIED 

CORRECT 1/3 O.B. SPL'S FOR CONVECTIVE AMPLIFICATION AT 
FORWARD SPEED AS FOLLOWS: 

C=-40 • *ALOG10 ( 1 . 0-MA*COS ( PHP/57 • 3 ) ) 
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101 DO 47 I=lf24 

102 SPL ( I ) =SPL ( I ) +C 

103 IF ( SPL ( I ) • LT# 0 « 0 ) SPL(I)=0.0 

104 47 CONTINUE 

105 C 

106 IF ( TC5 • EQ . 0 ) GO TO 70 

107 CALL DCTRT ( DSPL* XLTP t HDTP *MDTP t CDTP *FDCR*0.0*0«0* PHP r 4 ) 

108 DO 49 1=1 *24 

109 SPL ( I ) =SPL ( I ) -DSPL ( I ) 

110 IF(SPL(I) .LT.0.0) SPL ( I ) =0 • 0 

111 49 CONTINUE 

112 C 

113 70 RETURN 

114 END 
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3 

4 

5 

6 
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8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 C 

38 

39 C 

40 C 

41 C 

42 C 

43 C 

44 C 

45 C 

46 C 

47 

48 

49 

50 


SUBROUTINE JETMX { FREQ r TO r D r VJ1 r TJ1 r NPLU6 r NCOAX r X r Y ) 

DIMENSION A(4rl7) r 

2 ARAT ( 5 ) rHN(4) r FREQ (■ 24 ) r XOR ( 4 ) rY0R(4) rRS(4) rTHS{4) rRSTAR(4) r 

3 THSTAR (4) r P2PRF2 ( 4 ) r FI { 4 ) r F2 (4 ) 

REAL MC ( 5 r 17 ) r LGSTN ( 17 ) r LP1 ( 4 ) rLP2(4) r MA r M r NAR r MSTAR r ME r 
2 MCX 

INTEGER ERR 

COMMON / JMX/SPL1 (24) rSPL2(24) r TREF r RGr GC r H r RW r 
2 D2 r V J2 » T J2 f NAR 

COMMON MAjRHOAfCArRHOOrCOr C IMPDr DX1 r DX2 r ERR 

DATA ( (A(I»J) t J=lfl7) » 1=1 r4)/0. 8239 r-1. 24 r-2.3019f 1.4745 » 

2 45. 759 r-2. 1086 r 3. 3692 r-2. 3864 t 12. 808 »-l. 8167 r-1. 7894 rl. 1641 t 

3 4. 512 f -0. 028729 , 0.789 t -0.51772 t 2. 293E14r 

4 0. 4176 »-l. 7269 r-1. 3226 rl. 3507 r 18. 932 r -3. 2199 rl. 3276 r-0. 1765 r 

5 8. 4045 r -1. 5494 r -0. 16571 rO. 22007 r 5. 2846 rl. 3015 rO. 85116 r 

6 -0.60053rl. 993E13 r 

7 0. 32555 r 0. 2482 r -4. 0206 r 2. 6807 r 202. 88 r 3. 8424 r 10. 520 r -8. 0248 r 

8 10. 491 t -0. 019983 rO. 73551 r-0. 76495 t 4. 5428 rO. 20763 rO. 27691 t 

9 0 .036298 r6.086E13r 

1 0. 74038 r -3, 4538 r -4. 7181 r 3. 2662 r 331. 55 rO. 53294 r 9. 5983 r -6. 829 r 

2 10. 429r 0. 36732 r 0.89132 r-0. 87600 r 2. 6963 rl. 8687 r 3. 0174r 

3 -2 • 1329 r 4 • 794E13/ 

DATA ( (MC(Ir J) r J=lrl7) r 1=1 r 5 ) /2 . 00 r 1 .22 r 0 . 46 r 0 . 74 r 1 . 06 r 1 . 41 r 

2 1.78r2.16r2.55r2.95r3.36r3.80r4.24r4.70r5.00r5.10r5.10r 

3 0.85r-0.80rl.02rl.34rl.68r2.07r2.49r2.93r3.38r3.85r4.31r4.72r 

4 5.04r5.10r5.10r5.10r5.10r 

5 0 • 00 r 0 • 66 r 1 .45 r 2. 16 r 2 .82 r 3.47 r 4 . 07 r 4.62 r 4. 97 r5.04r5.06r5.08r 

6 5.10r5.10r5.10r5.10r5.10r 

7 0.00r0.74rl» 66 r2.46r3.15r3,74r4,23r4» 62 r4«91r5.04r5.06r5.08r 

8 5.10r5.10r5.10r5.10r5.10r 

9 0.95rl.60r2.63r3.60r4.32r4.77r4.99r5.02r5.03r5.04r5.06r5.08r 

1 5.10r5.10r5.10r5.10r5.10/ 

DATA AR AT/1, 0r3.0r6.0rl0.0rl5.0/ 

DATA LGSTN/-1 »15r-1.0r-0.8r-0.6r-0,4r-0.2r0.0r0,2r0.4r0,6r 

2 0.8rl.0rl.2rl.4rl.6rl.8r2.0/ 

DATA HN/1 . 25 r0.221r0« 442 r 0 . 884/ 

VA=MA*CA 


***** STEP 1 ***** 


SET PRIMARY OR SINGLE NOZZLE EXIT MACH NO.r ORr 

IF MA > O.Or COMPUTE RELATIVE MACH NO. FOR A SINGLE NOZZLE 

NSEC=0 

TJ=TJ1 

M=1 . 0/SQRT ( ( 1 . 33*RG*T J*GC ) / ( V J1*V J1 ) - ( 1 , 33-1 . 0 ) /2 . 0 ) 
IF(MA.LE.O.O) GO TO 15 
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51 IF(NCOAX.EQ.l) GO TO 15 

52 MSTAR=M*(1.0-VA/VJ1)**0.75 

53 M=MSTAR 

54 15 CONTINUE 

55 C 

56 C ***** STEP 2 ***** 

57 C 

58 C 

59 C COMPUTE EFFECTIVE NOZZLE DIAMETER t RD PARAMETER FOR PLUG 

60 C NOZZLES AND REDUCED COORDINATES OF RECEIVER LOCATION 

61 C 

62 RD=1.0 

63 DE=D 

64 IF(NPLUG.NE.l) GO TO 18 

65 DE=2.0*SQRT(H*(D-H) ) 

66 RD=2.0*H/DE 

67 18 CONTINUE 

68 C 

69 C ***** STEP 3 ***** 

70 C 

71 XR=X/DE 

72 YR=Y/DE 

73 19 CONTINUE 

74 IF(XR,GT. 30,0. OR.XR.LT, -30.0) GO TO 81 

75 I F ( YR , LT .O.O.OR.YR, GT , 30,0) GO TO 81 

76 C 

77 C ***** STEP 4 ***** 

78 C 

79 C 

80 C COMPUTE POLAR COORDINATES OF RECEIVER RELATIVE TO THE 

81 C SOURCE LOCATION AS A FUNCTION OF FREQUENCY t CORRECTED 

82 C FOR FORWARD SPEED IF NECESSARY. 

83 C 

84 DO 20 1 = 1 r 4 

85 X0R(I)=A<Irl3)*M**A{I»14)*(TJ/TREF)**(A(Irl5)+A{ Ir 16)*M) 

86 YOR ( I ) =0 .5+0 , 132*X0R ( I ) 

87 IF ( YR.GT ,YOR ( I ) ) GO TO 7 

88 CK=-DX1 

89 IF(XR.GE.CK) GO TO 81 

90 THS(I)=3. 14159265 

91 RS ( I ) =XOR ( I ) -XR 

92 THST AR ( I ) =THS ( I ) 

93 RST AR ( I ) =RS ( I ) / ( 1 , 0-MA ) 

94 GO TO 12 

95 7 THS ( I ) =ATAN ( ( YR-YOR < I ) )/(XR-XOR(I) ) ) 

96 IF (THS (I) .LT.0,0) THS ( I ) =THS ( I ) +3 . 14159265 

97 IF(THSd) ,LT. 0.1309) GO TO 81 

98 RS ( I ) =SQRT ( (XR-XOR(I) ) **2+ ( YR-YOR ( I ) )**2) 

99 C 

100 C ***** STEP 5 INTEGRATED INTO DO LOOP OF STEP 4 ***** 
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101 C 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 C 

114 C 

115 C 

116 C 

117 C 

118 C 

119 C 

120 
121 
122 

123 

124 

125 C 

126 C 

127 C 

128 C 

129 C 

130 C 

131 C 

132 C 

133 C 

134 C 

135 C 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 C 


1 THST AR ( I ) =THS ( I ) 

RSTAR ( I ) =RS ( I ) 

IF(MA.LE.O.O) 60 TO 12 

CTHSTR=(1.0/(1.0-MA*MA) )*( ( 1 . O/TAN ( THS { I ) ) ) -MA*SQRT ( 1 . 0-MA*MA+ 
2 ( 1. O/TAN ( THS ( I ) ) )**2) ) 

IF ( CTHSTR . NE . 0 . 0 ) THST AR ( I ) =ATAN ( 1 . O/CTHSTR ) 

IF ( THSTAR ( I ) . LT.0.0 ) THSTAR ( I ) =THS.TAR ( I ) +3 . 14159265 
I F ( CTHSTR . EQ . 0 . 0 ) THSTAR ( I ) =1 . 570796 
RSTAR ( I ) = ( YR-YOR ( I ) } /SIN (THSTAR ( I ) ) 

12 CONTINUE 
20 CONTINUE 

***** STEP 6 ***** 


COMPUTE PREDICTION EQN COEFFICIENTS WHICH ARE NOT A 
FUNCTION OF FREQUENCY 

21 C1=M**2 • 34 

C2=10.65*(TJ/T0 )**0.93 
C3=-15 . 18* ( T J/TO ) **1 . 11*M**0 .89 
C6=17.5*(TJ/(3.0*TREF) ) ** ( 0 .89* (M*M-1 . 0 ) ) 
C7=0.41*(TJ/(3.6*TREF) ) ** ( 0 .566* (M*M-1 . 0 ) ) 

***** STEP 7 ***** 


COMPUTE PREDICTION EQN COEFFICIENTS WHICH ARE A FUNCTION 
OF FREQUENCY AND THE MEAN SQUARE ACOUSTIC PRESSURE 
FOR EACH OF THE THREE OCTAVE BANDS AND THE OVERALL 
INCLUDING A CORRECTION FOR DYNAMIC AMPLIFICATION 
(EFFECTS OF FORWARD SPEED ON DIRECTIVITY ANGLE AND 
PROPAGATION PATH LENGTH ARE ALSO INCLUDED) 

DO 30 1=1 r 4 
R=RST AR ( I ) 

THET=THST AR ( I ) 

C4=A(I>5)*M**A(Ir6)*(TJ/TREF)**(A{If7)+A(Ir8)*M) 

C5=A ( I t 9 ) *M**A ( I r 10 ) * ( T J/TREF ) ** ( A ( I t 11 ) +A ( I r 12 ) *M ) 

ALP2=A ( I > 1 ) *M**A ( I r 2 ) * ( TJ/TREF )**(A(I»3)+A(I>4)*M) 

C9= ( (1.0+ALP2*M*M)/(ALP2*M*M+(1.0-(M*COS(THET) )/ 

2 (1.0+C6*EXP(-C7*R) ) )**2) )**(5,/2. ) 

C10=1.0+(C4*EXP(-C5*THET) ) / ( 1 . 0+C6*EXP ( -C7*R/4 , ) ) 

P2PRF2 ( I ) = ( ( A ( I t 17 ) *C9/C10 ) * ( TJ/TREF ) **1 . 54*M**4 . 0* ( 1 . O+COS ( THET ) 

2 **4 . 0 ) * ( Cl/ ( R*R ) +C2/R**4 . 0+C3/R**6 . 0 ) ) / ( 1 . 0-MA*COS ( 3 . 14159265 

3 -THET) ) 

30 CONTINUE 

IF ( NSEC # EQ . 1 ) GO TO 70 
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151 C 

152 C 

153 C 

154 C 

155 C 

156 C 

157 C 

158 

159 

160 
161 
162 

163 

164 C 

165 C 

166 C 

167 C 

168 C 

169 

170 

171 

172 

173 C 

174 C 

175 C 

176 C 

177 C 

178 C 

179 C 

180 
181 
182 

183 

184 

185 

186 

187 C 

188 C 

189 C 

190 C 

191 C 

192 C 

193 C 

194 

195 

196 

197 

198 

199 

200 C 


***** STEP 8 AND 17 ***** 


COMPUTE SPL IN EACH FREQUENCY BAND FOR THE PRIMARY 
OR SINGLE NOZZLE JET AND CORRECT FOR LOCAL ACOUSTIC 
IMPEDANCE 

RCC=C IMPD 
DO 42 1=1 f 4 

LP1 ( I ) =10 , *ALOG10 ( P2PRF2 ( I ) )+RCC 

IF (NPLUG.EQ.l ) LP1 ( I ) =10 . *ALOG10 ( P2PRF2 ( I ) ) +3. 0*ALOG10 ( 0 . 10+ 
2 2 • Q*H/DE) +RCC 
42 CONTINUE 


COMPUTE CENTER FREQUENCIES FOR EACH FREQUENCY BAND OF 
PRIMARY OR SINGLE NOZZLE JET 

DO 48 1=1 t 4 

FI ( I )=HN( I )*CA/(DE*RD**0.4) 

48 CONTINUE 

IF(NCOAX.NE.l) GO TO 55 
***** STEP 9 ***** 


FOR COAXIAL JETS f COMPUTE CORRECTION TO PRIMARY JET 
SPL * S TO ACCOUNT FOR THE PRESENCE OF THE SECONDARY JET 

DO 54 1=1 f 4 

ABCIS=AL0G10 (FI ( I ) *DE*RD**0 . 4/V J1 ) 

MCX=DT AB2 ( ABC IS f NAR f LGSTN f ARAT f17f5f1f1fMCf5f I ERR ) 
DLP=1 0 . *MCX*ALOG10 ( 1 . 0-V J2/V J1 ) 

LP1 ( I ) =LPl ( I ) +DLP 
54 CONTINUE 
GO TO 55 

***** STEP 11 ***** 


FOR COAXIAL JETSf COMPUTE EFFECTIVE SECONDARY JET 
FLOW PARAMETERS 

57 NSEC=1 
DE=D2 
RD=1 • 0 

TTE= ( T Jl+R W*T J2 ) / ( 1 . 0+RW ) 

V JE= ( V J1+RW*V J2 ) / ( 1 # 0+RW ) 

ME=SQRT ( 1 ./ ( (1.4*RG*TTE*GC)/(VJE*VJE)-(1.4-1.0)/2.0) ) 
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201 C 

202 

203 

204 

205 C 

206 C 

207 C 

208 

209 

210 
211 

212 C 

213 C 

214 C 

215 C 

216 C 

217 C 

218 C 

219 C 

220 
221 
222 

223 

224 

225 C 

226 C 

227 C 

228 C 

229 C 

230 C 

231 C 

232 C 

233 C 

234 

235 

236 

237 C 

238 C 

239 

240 C 

241 C 

242 C 

243 C 

244 

245 

246 

247 C 

248 C 

249 C 

250 C 


IF(MA.GT.O.O) ME=ME*(1.0-VA/VJE>**0.75 

M=ME 

TJ=TTE 

***** STEP 12 ***** 

X2=X+DX2 
XR=X2/DE 
YR=Y/DE 
GO TO 19 

***** STEP 8 FOR THE SECONDARY JET ***** 


FOR COAXIAL JETS* COMPUTE SECONDARY JET SPL*S IN EACH 
FREQUENCY BAND AND EACH BANDS CENTER FREQUENCY 
INCLUDING CORRECTION FOR LOCAL ACOUSTIC IMPEDANCE 

70 CONTINUE 
DO 75 1=1*4 

LP2 ( I ) =10 • *ALOG10 ( P2PRF2 ( I ) )+RCC 
F2 ( I ) =HN ( I ) *CA/DE 
75 CONTINUE 

CALL SUBROUTINE TO EXPAND FREQUENCY RANGE AND COMPUTE 
1/3 O.B. AND SPECTRAL LEVELS FOR 50 TO 10*000 HZ 


***** STEPS 10*14*15* AND 20 ***** 

SECONDARY JET: 

CALL EXTFRQ ( SPL2 * FREQ * LP2 { 2 ) *LP2(3) *LP2(4) * F2 ( 2 ) * F2 ( 3 ) * 
2 F2 ( 4 ) *CA*VJE*DE*XR*YR*ERR*F2MIN*F2MAX*RD) 

IF(ERR.EQ.l) GO TO 250 


GO TO 65 


PRIMARY OR SINGLE NOZZLE JET : 

55 CALL EXTFRQ ( SPL1 * FREQ * LP1 ( 2 ) * LP1 (3) *LP1(4) *F1(2) *F1(3) * 
2 FI (4) *CA*VJ1*DE*XR* YR * ERR * F1MIN * F1MAX * RD ) 

IF(ERR.EQ.l) GO TO 250 
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251 C 

252 C 

253 IF ( NCOAX « EQ . 1 ) GO TO 57 

254 GO TO 250 

255 C 

256 C 

257 65 CONTINUE 

258 C 

259 GO TO 250 

260 81 WRITE (6 r 300 ) 

261 300 FORMAT ( // » 2X » 44H0BSERVER LOCATION OUTSIDE ALLOWABLE ENVELOPE t 

262 2 38H FOR JET MIXING NOISE - RUN TERMINATED r / ) 

263 ERR=1 

264 250 RETURN 

265 END 
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1 

2 

3 

4 

5 

6 
7 

a 

9 

10 

n 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 C 

27 C 

28 

29 C 

30 C 

31 C 

32 

33 

34 

35 

36 

37 

38 C 

39 C 

40 C 

41 C 

42 C 

43 

44 

45 

46 

47 

48 

49 

50 


SUBROUTINE EXTFRQ ( SPL p FQ r LPA r LPB p LPC r FA r FB p FC r CO r VU t D r XR p 
2 YR t ERR r FM IN p FMAX p R ) 

DIMENSION F(20) rFQ(24) rXD<5) r SPL { 24 ) 

REAL LSA r LSX rLSCrLSYr LP (20 ) rLF(20) rLSTN(31) rLR(5r31) r 
2 LPApLPB,LPC»LPXfLPY 
INTEGER ERR 

DATA LSTN/~1.40f-1.30f-1.22r--1.15r-1.10r-1.04r-1.00f-0.92p-~0.80r 

2 -0.70r-0.60p-0.52f-0.40p-0.30r-0.22f-0.15r-0.10r--0.04p0.00p 

3 0.08r0.20r0.30r0.40r0.48r0.60f 0.70 r 0.78 r 0.84 r 0.90 r 0.95 r 1.00/ 
DATA XD/-30. r5. rlO. r20. r30./ 

DATA ( (LR(IrU) pJ=1p31) p 1=1 p 5 ) /23 . 8 p 21 . 9 p 20 . 5 p 19 .2 p 18 .2 p 17 . 3 p 16 . 6 i 

2 15.3p13.4p11.9p10.4p9.2p7.5p6.2p5.2p4.5p3.8p3.4p3.0p2.6p2.2p 

3 2.2p2.6p2.9p3.6p4.2p4.9p5.5p6«0p6.6p7.1p 

4 23.8 p21.9p 20.5 p 19.2 p 18.2 p 17.3 p 16.6 p 15.3 p 13.4 p 11.9 p 10.4 p 9.2 p 

5 7.5p6.2p5.2p4.5p3.8p3«4p3.0p2.6p2.2p2.2p2.6p2.9p3.6p4.2p4,9p 

6 5 • 5 p 6 . 0 p 6 • 6 p 7. 1 p 

7 19.9p18.1p16.8p15.6p14,6p13.8p13.0p11.9p10.0p8.6p7.4p6.5p5.3p 

8 4.3p3.8p3«2p3.0p2.8p2.6p2.7p3.2p3.8p4,7p5.4p6.8p7.9p9.0p 

9 10.0p10.9p11.7p12.4p 

1 16.5p14.9p13.3p12.3p11.4p10.6p9.9p8.9p7.3p6.2p5.3p4.8p4.1p3,8p 

2 3.8p4.0p4.2p4.6p5.0p5.6p6.9p8.0p9.3p10.4p12.4p14.0p15.5p16.6p 

3 17.7p18.7p19.5p 

4 13«6p11.4p9.9p8.6p7.6p6»8p6.2p5.2p4.1p3.6p3.6p3«8p4.7p5.8p6.8p 

5 7.6p8.5p9.2p10.0p11.3p13.6p15.5p17.4p19.2p21.9p24.1p26.0p27.5p 

6 28.9p30.1p31.4/ 


N=3 

SET CENTER FREQS FOR THE THREE PREDICTED OCTAVE-BANDS 

F ( 1 ) =FA 
F (2)=FB 
F ( 3 ) =FC 
LP ( 1 ) =LPA 
LP ( 2 ) =LPB 
LP ( 3 ) =LPC 

EXTEND THE LOW-FREQ RANGE TO AT LEAST 50 HZp IF NECESSARY 


IF ( FA , LE.50 . 0 ) GO TO 40 
SA=0 ,221*C0/VU 
LSA=ALOG10 (SA) 

IF (LSA.LT , -1,40. OR, LSA, GT ,1.0 ) GO TO 89 
DLPA=DTAB2(LSApXRpLSTNpXDp31p5p1p1pLRp5p I ERR) 
7 FX=F(l)/2.0 
SX= ( FX/FA ) *SA 
LSX=ALOG10 (SX) 
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51 IF ( LSX. LT .—1.40) WRITE(6,87> 

52 DLPX=DTAB2 ( LSX , XR , LSTN , XD , 31 , 5 , 1 , 1 , LR , 5 , 1 ERR ) 

53 LPX=LPA-DLPX+DLPA 

54 IF(LPX.LT.O.O) LPX=0.0 

55 N=N+1 

56 DO 10 NK=2,N 

57 I=N+2-NK 

58 F(I)=F(I-1) 

59 10 LP(I)=LP(I-1) 

60 F(1)=FX 

61 LP(1)=LPX 

62 IF(F(1) .GT.50.0) GO TO 7 

63 C 

64 C EXTEND HI-FREQ RANGE TO AT LEAST 10000 HZ, IF NECESSARY 

65 C 

66 40 IF(F(N) .GE. 10000.0) GO TO 50 

67 SC=0.884*C0/VJ 

68 LSC=AL0G10(SC) 

69 IF ( LSC • GT« 1 , 0 ) GO TO 89 

70 DLPC=DTAB2(LSC,XR,LSTN, XD , 31 , 5 , 1 , 1 , LR , 5 , IERR) 

71 42 FY=2.0*F<N) 

72 SY=(FY/FC)*SC 

73 LSY=ALOG10 (SY) 

74 DLPY=DTAB2(LSY,XR,LSTN,XD,31,5,1,1,LR,5, IERR) 

75 LPY=LPC-DLPY+DLPC 

76 IF(LPY.LT.O.O) LPY=0.0 

77 N=N+1 

78 F(N)=FY 

79 LP(N)=LPY 

80 IF(F(N) .LT. 10000.0) GO TO 42 

81 C 

82 C FROM THE NOW GENERATED OCT-BAND SPL SPECTRUM COVERING THE 

83 C FREQ RANGE FROM 50 TO 10000 HZ r INTERPOLATE THE OCTAVE- 

84 C BAND LEVELS AT EACH 1/3 OCT-BAND CENTER FREQ FROM 50 TO 

85 C 10000 HZ AND CORRECT TO 1/3 OCT-BAND LEVELS BY SUBTRACTING 

86 C 4.85 DB 

87 C 

88 C 

89 50 DO 58 I=lfN 

90 58 LF(I)=ALOG10(F(I) ) 

91 DO 60 1=1,24 

92 FK=ALOG10 ( FQ ( I ) ) 

93 OBSPL=GIRC(FK,LF,LP,N,l) 

94 SPL ( I ) =0BSPL-4.85 

95 60 IF ( SPL (I) .LT.0.0) SPL(I)=0.0 

96 C 

97 87 FORMAT (//, 2 X,48H WARNING - IN JET MIXING NOISE, ATTEMPT TO EXPAND, 

98 2 24H LOW FREQ RANGE TO 50 HZ , / , 12X , 18HEXCEEDED LIMITS OF, 

99 3 24H THOMSON SPECTRUM SHAPES,/) 

100 GO TO 100 
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101 

102 

103 

104 

105 

106 

107 

108 


89 WRITE(6,80) 

ERR=1 

80 FORMAT ( // , 2X , 46HF0R JET MIXING NOISE, ONE OR MORE OF THE THREE, 

2 50H PREDICTED OCTAVE BANDS FALLS OUTSIDE THE STROUHAL, 

3 / , 2X , 40H NO. RANGE OF THE REDUCED SPECTRA SHAPES, 

4 17H - RUN TERMINATED,/) 

100 RETURN 

END 
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1 

2 C 

3 C 

4 C 

5 

6 

7 

8 C 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 C 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 C 

32 C 

33 C 

34 

35 

36 

37 

38 

39 C 

40 C 

41 C 

42 C 

43 C 

44 

45 

46 C 

47 C 

48 C 

49 

50 C 


SUBROUTINE SHOCK ( FREQ p VJ t TJt DE p HJ p SC p RG p GC p X p Y p SPL ) 

***** SHOCK BROADBAND NOISE FOR THE LFC PROGRAM ************ 


DIMENSION FREQ (24) p SX (20 ) p HX ( 20 ) p CX ( 20 ) p SPL (24) 

REAL KO p K1 p M J pLOpLIpMCpMAp LOGS I 
INTEGER SpSCpERR 

COMMON MA pRHOApCApRHOOpCOpCIMPDpDXIp DX2 p ERR 

DATA SX/-0. 6990 p -0. 5229 p -0. 3979 p -0. 1549 pO.OpO. 1761 pO. 3010 t 

2 0. 4771 p 0. 5441 p 0. 6020 p 0. 6532 pO. 6990 pO. 7782 pO. 8451 pO. 9031 p 

3 1 • 000 p 1 « 301 p 1 .602 p 1 » 832 p 1 « 845/ 

DATA HX/116.0p121.6p125.5p132.5p137.7p142.7p145.7p148.5p 

2 149. 1 p 149. 2 p 149. 1 p 148. 8 p 147. 9 p 146. 7 p 145. 7 p 143. 7 p 137. 4 p 

3 130 • 5 p 125 .4 p125.2/ 

DATA CX/0.70p0.71p0.71p0.72p0.73p0.74p0.74p0.71p0.69p0.67p 
2 0.64p0.62p0 .58 p 0 , 54 p0.50p0. 45 p0.28p0,12p0.02p0. 02/ 


DO 10 1=1 p 24 
SPL ( I ) =0 • 0 
10 CONTINUE 
C=0.70 
K0=1.10 
Kl=l,31 
NS=8 

BC=0.2316 

GAMA=1.33 

IF (SC.EQ.8) GAMA=1.4 

MJ=1 . 0/SQRT ( (GAMA*RG*TJ*GC)/ (VJ*VJ)- (GAMA-1.0 )/2.0) 
IF(MJ.LE.l.O) GO TO 70 

***** STEP 1 ***** 

BET A=SQRT ( M J*M J-l . 0 ) 

L0=K0*H J*BET A 
L1=K1*HJ*BETA 
MC=C*VJ/CA 
VC=C*VJ 

***** STEP 2 ***** 

***** COMPUTE EMISSION ANGLE AND ACTUAL SOUND PATH LENGTH ***** 

CALL XFORM(Xp Y p 0 . 0 p MA » SC p DX1 p DX2 p PH p PHP p R p RP p ERR ) 

IF ( ERR , EQ • 1 ) GO TO 70 

***** STEP 3 ***** 

W1=6,2832*BC*HU*BETA/CA 
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51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 


IF ( BETA. 6T • 1 • 0 ) 60 TO 101 

ANS1=(40.0*ALOG10 (BETA) ) - (20 . 0*ALOG10 (RP/DE ) ) 

GO TO 105 

101 ANS1=(20,0*AL0G10 (BETA) )-( 20 . 0*AI_OG10 ( RP/DE ) ) 

C 

105 CONTINUE 
C 

C ***** STEP 4 ***** 

C 

DF=1 . 0-MC*COS ( ( 180 . -PHP ) /57 .29577951 ) 

W2=L1*DF/VC 

C 

C 

C ***** STEP 5 ***** 

C 

DO 666 1=1 p 24 

SI=2.0*3.14159*FREQ( I )*L0/CA 

IF(SI.LT.0.2,OR.SI.GT.70.0) WRITE (6 p 1500 ) SCpI 
1500 FORMAT (// p 2X p 35H WARNING - FOR SHOCK BROADBAND NOISEp 

2 10H p SOURCE =pI2p 24H p THE STROUHAL PARAMETER p / p 12 X p 

3 16H OF 1/3 O.B. NO,pI3p 20H IS OUTSIDE RANGE OFp 

4 15H MASTER SPECTRAp/) 

LOGSI=ALOG10 (SI ) 

HI=GIRC ( LOGS I pSXpHXp20p1) 

C I=GIRC ( LOGS I pSXpCXp20p1) 

C 

c 

c 

c 

WC=2.0*3.14159*FREQ( I ) 

1 1 END=NS-1 
SUMI=0.0 

DO 210 NzlpIIEND 
CI2=CI**(N*N) 

I SEND=NS-N 
SUMS=0 • 0 

DO 220 S=1pISEND 
IS=S-1 

QNS=W2*N*(1.0-(0.06*( IS+( (N+1.0J/2.0) ) ) ) 

QCOS=COS ( QNS*WC ) 

QSIN=SIN( (QNS*WC*BC)/2.0) 

W0RK3= ( QCOS*QS IN ) /QNS 
SUMS=SUMS+W0RK3 
220 CONTINUE 

W0RK4=CI2*SUMS 
SUMI=SUMI+W0RK4 
210 CONTINUE 
C 

W0RK5= (4. 0*SUMI ) / ( NS*BC*WC ) 

ANS2=1 • 0+W0RK5 
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101 

102 C 

103 

104 C 

105 

106 C 

107 

108 C 

109 C 

110 C 

111 C 

112 

113 

114 

115 

116 

117 

118 C 

119 

120 


ANS2=10 • 0*ALOG10 ( ABS ( ANS2 ) ) 
ANS3=10 • *ALOG10 ( W1*FREQ ( I ) ) 
• SPL(I)=HI+ANS1+ANS3+ANS2 
666 CONTINUE 


*** STEP 5(G) ***** 

CAMP=-40 • *ALOG10 ( 1 • 0-MA*COS ( PHP/57 .29577951 ) ) 
DAMP=-10 • *ALOG10 ( 1 . 0-MA*COS ( PHP/57 .29577951 ) ) 
DO 30 I=lr24 

SPL ( I ) =SPL ( I ) +C IMPD+CAMP+DAMP 
IF(SPL(I) .LT.0.0) SPL ( I ) =0 • 0 
30 CONTINUE 

70 RETURN 
. END 
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1 SUBROUTINE SCRCH ( V J » H J t DJt NPR t A t AREF t X t Y t SC r SPL1 t SPL2> NB1 t NB2 ) 

2 REAL NPR . MA. PR { 9 ) t SPLRF1 (9 ) . SPLRF2 ( 9 ) 

3 INTEGER SC t ERR 

4 COMMON MArRHOA»CArRHOOrCOrCIMPDf DXl»DX2f ERR 

5 C 

6 DATA PR/2. 0 > 2. 5 r 3. 0 f 3 . 5 > 4 . 0 ^ 4 . 5f 5. 0r6. 0>7.0/ 

7 DATA SPLRF1/110. t 127. »136. f141. t 143. 5» 144. 5f144.5f143. t 140.5/ 

8 DATA SPLRF2/110 • t 124 • >130, rl33, rl 34. 5.135, fl34,5r 132 • t 128 . / 

9 C 

10 IF ( NPR , LT .2 , 0 ) GO TO 30 

11 XS=3.5*HU*SQRT<NPR-1.893) 

12 C 

13 CALL XFORM(X>YfXS»MAfSCfDXl> DX2 fPHfPHP^RfRP, ERR ) 

14 IF (ERR.EQ.l ) GO TO 30 

15 C 

16 C STEP 3 

17 C 

18 SPLlRF=GIRC(NPR,PR,SPLRFlr9» 1) 

19 SPL2RF=G IRC ( NPR r PR , SPLRF2 1 9 1 1 ) 

20 C 

21 SPL1R=SPL1RF+10,*AL0G10<A/AREF) 

22 SPL2R=SPL2RF+10.*AL0G10{A/AREF) 

23 C 

24 C STEP 4 

25 C 

26 RREF=4 , 0*D J 

27 DSPL1=20.*ALOG10 ( ( COS (PHP/ ( 2 , *57 . 29577951 ) ) )**2+ 

28 2 0.5*(SIN(PHP/(2. *57, 29577951) ) ) **2 ) -20 . *ALOG10 (RP/RREF ) +2. 5 

29 C 

30 DSPL2=20.*ALOG10 (SIN (PHP/57. 29577951) ) -20 , *ALOG10 (RP/RREF ) 

31 C 

32 C STEP 5 

33 C 

34 SPL1=SPL1R+DSPL1 

35 SPL2=SPL2R+DSPL2 

36 C 

37 C STEP 6 

38 C 

39 DSPL=C IMPD-40 , *ALOG10 <1,0-MA*COS( PHP/57 « 29577951 ) ) 

40 SPL1=SPL1+DSPL 

41 SPL2=SPL2+DSPL 

42 IF(SPLl.LT.O.O) SPL1=0.0 

43 IF ( SPL2 . LT. 0 . 0 ) SPL2=0.0 

44 C 

45 C STEP 7 

46 C 

47 MJA=VJ/CA 

48 C 

49 F1=(CA*(MA+0.625*(MJA-MA) ) * ( 1 . 0-MA ) * ( 1 . 0+0 . 625*M JA ) )/ 

50 2 (1.25*SQRT(NPR-1,893)*(1.0+0,625*(MJA-MA) )*MJA*HJ) 
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51 C 

52 F2=2.0*F1 

53 IF(F1.GE.44.7.AND.F1.LE.11220. ) GO TO 12 

54 WR ITE ( 6 * 100 ) 

55 SPL1=0.0 

56 NB1=0 

57 GO TO 14 

58 12 CALL BNDN(F1*24*NB) 

59 NB1=NB 

60 14 IF(F2.GE.44.7.AND.F2.LE.11220. ) GO TO 16 

61 WRITE (6*102) 

62 SPL2=0.0 

63 NB2=0 

64 GO TO 18 

65 16 CALL BNDN(F2*24*NB) 

66 NB2=NB 

67 18 CONTINUE 

68 100 FORMAT (/*2X*40HF0R SHOCK SCREECH TONES * THE FUNDAMENTAL* 

69 2 45H FREQUENCY IS OUTSIDE RANGE OF 50 TO 10000 HZ*/) 

70 102 FORMAT ( // * 2X * 44HF0R SHOCK SCREECH TONES* THE SECOND HARMONIC* 

71 2 45H FREQUENCY IS OUTSIDE RANGE OF 50 TO 10000 HZ*/) 

72 C 

73 30 RETURN 

74 END 
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1 SUBROUTINE DCTRT ( DSI > L > HD > MD > CD > FD >NS > MI > PHP > IC ) 

2 REAL FQ ( 24 ) > DSI (24 ) rLf MD> MI > NS > M 

3 DATA FQ/50 • *63. *80* > 100 • > 125 • > 160 • > 200 • > 250 • >315* > 400 • > 500 • > 

4 2 630. >800. >1000. >1250. >1600. >2000. >2500. >3150, >4000. >5000. > 

5 3 6300. >8000. >10000./ 

6 DO 5 1=1 >24 

7 DSI ( I ) =0 • 0 

8 5 CONTINUE 

9 IF ( IC.EQ.O ) GO TO 35 

10 IF ( IC.EQ.2) GO TO 25 

11 M=MD 

12 IF ( IC.NE.4) M=-MD 

13 DDBPK=10 . * ( L/HD ) **0 • 7 

14 XLAMD=CD/FD 

15 Fl=l , 0/ ( HD/XL AMD ) **0 . 6 

16 F2 = l . 0- ( M/2.0 ) * (2 . 0- ( HD/XLAMD ) ) 

17 C 

18 IF ( IC.EQ.4) GO TO 10 

19 F3=(140.-PHP)/80. 

20 IF ( PHP . LE, 60 • 0 ) F3= ( ( 4 . 0-NS ) /4 . ) *( PHP/100 . ) +0 . 15*NS+0 . 4 

21 GO TO 15 

22 C 

23 10 F3=PHP/130 • 

24 IF (PHP.GT ,130. ) F3=(205.-PHP)/75. 

25 15 CONTINUE 

26 C 

27 DDBFD=DDBPK*F1*F2*F3 

28 C 

29 DO 20 1=1 >24 

30 DSI ( I ) =DDBFD 

31 IF(FQ(I) .GT.FD.OR.FQ(I) .LT.FD) DSI(I)= 

32 2 DDBFD*EXP<-( (ABS(ALOG10(FQ(I)/FD) ) )**1.3/0.35) > 

33 20 CONTINUE 

34 IF ( IC.NE.3) GO TO 35 

35 25 DSHMI=0 . 0 

36 IF ( MI , LE. 0 • 5 ) GO TO 35 

37 DSHMI=216.*(MI-0.5)**2.5 

38 DO 30 1=1 > 24 

39 DSI ( I ) =DSI ( I ) +DSHMI 

40 30 CONTINUE 

41 35 CONTINUE 

42 RETURN 

43 END 
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1 SUBROUTINE TEN ( FREQ t L r DELTA r C r MUA t VREF t X t Y t SPL ) 

2 DIMENSION FREQ (24) t SPL (24) 

3 REAL Lr MA»MUA 

4 INTEGER ERR 

5 COMMON MArRHOAr CArRHOO r CO r CIMPDr DX1 * DX2 > ERR 

6 C 

7 VA=MA*C A 

8 C 

9 C ***** STEP 1 ***** 

10 C 

11 C ***** COMPUTE ESTIMATED TURB. BL THICHNESS AT T.E. ***** 

12 C 

13 IF(DELTA.LE.O.O) DELTA=0 .376*0/ ( VA*C*RHOA/MUA ) **0 . 2 

14 C 

15 C ***** STEP 2 ***** 

16 C 

17 C ***** CALCULATE FORWARD SPEED TRANSFORMATIONS ***** 

18 C 

19 CALL XFORMtXrYf O.OrMAr llrDXlf DX2fPHfPHPrRfRPf ERR) 

20 IF(ERR.EQ.l) GO TO 910 

21 C 

22 C ***** STEP 3 ***** 

23 C 

24 BET AP=2» 0*57 .29577951* AT AN ( L/ (2. 0*RP ) ) 

25 C 

26 C ***** STEP 4 ***** 

27 C 

28 C 

29 C ***** COMPUTE THE OVERALL SOUND PRESSURE LEVEL ***** 

30 C 

31 0ASPL1=50 • *AL0G10 ( VA/VREF ) +10 • *ALOG10 ( DELT A*BET AP/RP ) 

32 2 +10. *ALOG10( (COS( ( PHP/2 . 0 ) /57 . 29577951) )**2)+80, 7 

33 C 

34 C ***** STEP 5 ***** 

35 C 

36 C ***** CORRECT OASPL FOR LOCAL ACOUSTIC IMPEDANCE ***** 

37 C ***** AND CONVECTIVE AMPLIFICATION ***** 

38 C 

39 CAMP=-40 »*ALOG10 ( 1 . 0-MA*COS (PHP/57 .29577951 ) ) 

40 0ASPL1=0ASPL1+C IMPD+CAMP 

41 C 

42 C ***** STEP 6 ***** 

43 C 

44 c ***** COMPUTE THE PEAK FREQUENCY ***** 

45 C 

46 FPK=0.1*VA/DELTA 

47 C 

48 C ***** STEP 7 ***** 

49 C 

50 C ***** FIND 1/3 O.B. CONTAINING THE PEAK FREQ t FPK ***** 
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51 C 

52 DO 6 1=1,24 

53 SPL ( I ) =0 • 0 

54 6 CONTINUE 

55 IF(FPK.LT.44.7.0R.FPK.GT.11220. ) GO TO 905 

56 CALL BNDN ( FPK , 24 , NB ) 

57 IREF=NB 

58 C 

59 C ***** COMPUTE THE 1/3 O.B. SPECTRUM ***** 

60 C 

61 DO 30 1=1,24 

62 SPL ( I ) =OASPL1 + 10 « *AL0G1Q ( 0 . 613* ( ( FREQ ( I ) /FREQ ( IREF ) )**4)/ 

63 2 ( ABS ( ( FREQ ( I ) /FREQ ( IREF ) ) **1 , 5+0 .5 ) ) **4 ) 

64 IF(SPL{I) .LT.0.0) SPL(I)=0.0 

65 30 CONTINUE 

66 C 

67 C ***** STEP 8 SHIELDING - - NOT APPLIED ***** 

68 GO TO 910 

69 905 WRITE (6,200) 

70 200 FORMAT ( // , 2X , 43HPEAK FREQ OF TRAILING-EDGE NOISE IS OUTSIDE, 

71 2 45H RANGE OF 50 TO 10000 HZ - NO T.E.N. COMPUTED,/) 

72 910 RETURN 

73 END 
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1 SUBROUTINE TBLN ( X t Y , 2 t A t L t DELS t NA J t MUA t IC ) 

2 C ***** TURBULENT BOUNDARY LAYER NOISE PREDICTION ***** 

3 C 

4 C 

5 DIMENSION A(20) » X (20 ) * Y (20 ) » Z (20 ) t DELS (20) r SPL (24 > 20 ) 

6 REAL MAf MUAf L (20 ) 

7 INTEGER ERR 

8 COMMON /TBL/UREF,AREF,RREF,SPLT(24) 

9 COMMON MArRHOAf CArRHOOf COrCIMPDf DXlf DX2f ERR 

10 c 

11 c ***** BEGIN ELEMENTAL AREA LOOP ***** 

12 C 

13 U=MA*CA 

14 DO 6 J=1 t NAU 

15 DO 4 1 = 1 r 24 

16 SPL ( I > J ) =0 • 0 

17 4 CONTINUE 

18 6 CONTINUE 

19 DO 200 J=1 t NA J 

20 C 

21 C ***** STEP 1 ***** 

22 C 

23 C ***** FORWARD SPEED TRANSFORMATION ***** 

24 C 

25 IF(Y(U) .LE.0.0) GO TO 200 

26 XX=X(J) 

27 YY=SQRT(Y( J)*Y( J)+Z( J)*Z( J) ) 

28 CALL XFORMtXXr YYrO.O^MAf 12f DXl»DX2»PHfPHP,R»RP f ERR) 

29 IF (ERR • EQ. 1 ) GO TO 910 

30 C 

31 C ***** STEP 1(D) ***** 

32 C 

33 GMP=57 • 29577951* ACOS ( Y ( J ) /RP ) 

34 C 

35 C ***** STEP 2 ***** 

36 C 

37 OASPL=0 • 0 

38 IF ( GMP . LT • 90 « 0 ) OASPL=60 . *ALOG10 ( U/UREF ) +10 . *AL0G1 0 ( A ( J ) /AREF ) 

39 2 -20 »*ALOG10 ( RP/RREF ) +20 • *ALOG10 ( COS ( GMP/57 .29577951) ) 

40 3 +91.1 

41 C 

42 C ***** STEP 3 ***** 

43 C 

44 c ***** CORRECT FOR CONVECTIVE AND DYNAMIC AMPLIFICATION ***** 

45 C 

46 CAMP=-30.*ALOG10 (1. 0-0. 18*MA*C0S ( PHP/57 . 29577951) ) 

47 DAMP=-10 » *ALOG10 ( 1 • 0-MA*C05 ( PHP/57 • 29577951 ) ) 

48 IF ( GMP , LT, 90 , 0 ) OASPL=OASPL+CAMP+DAMP 

49 C 

50 C ***** STEP 4 ***** 
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51 C 

52 C ***** CORRECT FOR LOCAL ACOUSTIC IMPEDANCE ***** 

53 C 

54 IF ( GMP .LT.90.0) OASPL=OASPL+C IMPD 

55 C 

56 IF(OASPL.LE.O.O) GO TO 200 

57 C 

58 C ***** STEP 5 ***** 

59 C 

60 IF(IC.EQ.l) GO TO 10 

61 DELTA=0.376*L( J ) / ( RHOA*U*L ( J)/MUA)**0.2 

62 DELTAS=DELTA/8,0 

63 GO TO 15 

64 10 DELT AS=DELS ( J ) 

65 15 FPK=0.01102*U/DELTAS 

66 C 

67 C ***** STEP 6 ***** 

68 C 

69 C 

70 C ***** FIND THE 1/3 OB CONTAINING THE PEAK FREQ # FPK ***** 

71 C 

72 IF ( FPK ■ LT ■ 44 • 7 • OR • FPK .GT . 11220 . ) GO TO 40 

73 CALL BNDN(FPK#24#NB) 

74 NU=NB 

75 C 

76 SPLPK=0ASPL-7.0 

77 DO 38 1=1 #24 

78 IF(I.EQ.NJ) SPLd# J)=SPLPK 

79 IFCI.EQ. (NJ-1) ) SPL( I # J5=SPLPK-1.0 

80 IF d • EQ« ( N J-2 ) ) SPL( I# J5=SPLPK-2.0 

81 IF ( I «LT • (NJ— 2 ) ) SPLd# J } =SPLPK-2 . 0- ( N J-2-I ) *2 . 7 

82 IF(I.GT.NJ) SPLd# J ) =SPLPK+ < NJ-I ) *2 .2 

83 IF(SPL(I»J) .LT.0.0) SPLd#U)=0.0 

84 38 CONTINUE 

85 C 

86 GO TO 200 

87 40 WRITE { 6 » 300 ) J 

88 300 FORMAT (/ #2X#46HF0R TURB. B.L. NOISE# PEAK FREQ OF ELEMENT NO.# 

89 2 13 # 17H IS OUTSIDE RANGE#/#2X#29H OF 50 TO 10000 HZ - NO NOISE# 

90 3 39H CONTRIBUTION COMPUTED FOR THIS ELEMENT#/) 

91 C 

92 200 CONTINUE 

93 C 

gq. c ***** STEPS 8 AND 9 ***** 

95 C 

96 C ***** COMPUTE TOTAL 1/3 OB SPECTRUM AND OASPL ***** 

97 C 

98 DO 48 1=1 # 24 

9 Q SPLT ( I ) =0 « 0 

100 48 CONTINUE 
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101 DO 50 1=1 f 24 

102 SUM=0.0 

103 DO 58 J=1 p NA J 

104 IF(SPL( If J) .GT.0.0) SUM=SUM+10 . ** ( SPL ( I t J ) /10 . ) 

105 58 CONTINUE 

106 IF ( SUM. GT .0.0) SPLT ( I )=10«*ALOG10 (SUM) 

107 50 CONTINUE 

108 C 

109 C 

110 910 CONTINUE 

111 RETURN 

112 END 
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1 SUBROUTINE XFORM ( XP , YA r XS »M t SC t DX1 t 0X2 r PH ? PHP t R t RP t ERR ) 

2 C 

3 C AIRPLANE MACH NO. (M) MUST BE LESS THAN 1.0 

4 C 

5 INTEGER SC » ERR 

6 REAL M 

7 Y=YA 

8 IF(SC.EQ.1.0R.SC.EQ.3) GO TO 10 

9 IFCSC.EQ.4.0R.SC.EQ.5) GO TO 20 

10 IF(SC.EQ.7.0R.SC.EQ.9) GO TO 30 

11 IF(SC.EQ.8.0R.SC.EQ.10) GO TO 40 

12 IF (SC,EQ.11,0R,SC.EQ,12) GO TO 20 

13 C SC=2 

14 X=XP-DX2 

15 GO TO 50 

16 10 X=XP-DX1 

17 GO TO 50 

18 20 X=XP 

19 GO TO 50 

20 30 X=XP+XS 

21 GO TO 50 

22 40 X=XP-DX2+XS 

23 50 CONTINUE 

24 C 

25 C GIVEN X AND Y t COMPUTE R AND PH 

26 C G.IVEN R, PHr AND Mr COMPUTE RP AND PHP 

27 C 

28 IF (X.EQ.0.0. AND. Y , EQ , 0 . 0 ) GO TO 120 

29 IF(Y.GT.O.O) GO TO 90 

30 IF(X.GT.O.O) GO TO 82 

31 PH=180.0 

32 PHP=180.0 

33 R=-X 

34 RP=R/C1.0+M) 

35 GO TO 150 

36 82 PH=0.0 

37 PHP=0.0 

38 R=X 

39 RP=R/(1.0-M) 

40 GO TO 150 

41 90 IF(X.GT.O.O.OR.X.LT.O.O) GO TO 92 

42 PH=90 • 0 

43 R=Y 

44 GO TO 94 

45 92 R=SQRT(X*X+Y*Y) 

46 PH=57 • 29577951* AT AN ( Y/X ) 

47 IF(PH.LT.O.O) PH=PH+180.0 

48 94 CONTINUE 

49 IFCM.GT.0.0) GO TO 98 

50 PHP=PH 
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51 RP=R 

52 GO TO 150 

53 98 COTPH=1.0/TAN(PH/57. 29577951) 

54 COTPHP=(1.0/(1.0-M*M) ) * ( COTPH+M*SQRT ( 1 . 0-M*M+COTPH*COTPH ) ) 

55 PHP=90.0 

56 IF(COTPHP.GT. 0.0. OR. COTPHP.LT. 0.0) PHP=57 . 29577951* 

57 2 AT AN ( 1 . O/COTPHP ) 

58 IF(PHP.LT.O.O) PHP=PHP+180 . 0 

59 RP=R*SIN (PH/57.29577951) /SIN (PHP/57. 29577951) 

60 GO TO 150 

61 120 WRITE(6f 122) SC 

62 122 FORMAT (//r5Xp ’FOR SOURCE CODE =*»I3r» THE OBSERVER LOCATION * t 

63 2 * IS SAME AS SOURCE LOCATION - RUN TERMINATED *»/ ) 

64 ERR=1 

65 150 RETURN 

66 END 
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1 SUBROUTINE BNDN ( FX t NTOB t BN ) 

2 INTEGER BN 

3 REAL F ( 27 ) 

4 DATA F/50. >63. >80. >100. #125. >160. >200, >250.^315. >400. >500., 

5 2 630. >800. , 1000. , 1250. >1600. >2000. >2500. ,3150. ,4000. >5000. , 

6 3 6300. ,8000. , 10000. >12500. >16000, *20000./ 

7 C 

8 KK=NT0B-1 

9 BNsNTOB 


10 


DO 10 J=1,KK 

11 


FC=SQRT(F(U)*F( J+l ) ) 

12 


IF(FX.GT.FC) GO TO 10 

13 


BN=J 

14 


GO TO 12 

15 

10 

CONTINUE 

16 

12 

CONTINUE 

17 


RETURN 

18 


END 
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1 


FUNCTION GIRC(ARG*X*Y*NX*IO) 

2 

C 

AITKIN INTERPOLATION 

3 

C 

ARG=INDEPENDENT ARGUMENT 

4, 

C 

X=INDEPENDENT TABLE 

5 

C 

Y=DEPENDENT TABLE 

6 

C 

NX=NUMBER VALVES X TABLE 

7 

C 

10=1 FIRST ORDER INTERPOLATION 

8 

C 

10=2 SECOND ORDER INTERPOLATION 

9 


DIMENSION X ( 1 ) ,Y(1) *XX(4) *YY(4) *EE(3) *FF(2) 

10 


IEND=NX— IO+l 

11 


IF ( X ( 1 ) -X ( 2 ) ) 10 * 20 * 20 

12 

C 

ASCENDING ORDER 

13 

10 

DO 15 1=1 * IEND 

14 


IF ( X ( I ) -ARG ) 15*30*30 

15 

15 

CONTINUE 

16 

16 

IB=NX-IO 

17 


I EX=0 

18 


GO TO 45 

19 

20 

DO 25 1=1 » IEND 

20 


IF(X(I)-ARG)30*30*25 

21 

25 

CONTINUE 

22 


GO TO 16 

23 

30 

IF ( IO-l ) 35*35 *40 

24 

35 

IF(I-1)36, 36*37 

25 

36 

IB=1 

26 


GO TO 45 

27 

37 

IB=I-IO 

28 


GO TO 45 

29 

40 

IF( 1-2)41*41*42 

30 

41 

IEX=0 

31 


GO TO 36 

32 

42 

I EX=1 

33 


GO TO 37 

34 

45 

DO 50 1=1*4 

35 


XX ( I ) =X ( IB ) —ARG 

36 


YY ( I ) =Y ( IB ) 

37 

50 

IB=IB+1 

38 


DO 60 1=1*3 

39 

60 

EE ( I ) =XX ( 1+1 ) -XX ( I ) 

40 


DO 70 1=1*2 

41 

70 

FF ( I ) =EE { I ) +EE ( 1+1 ) 

42 


DO 80 1=1*3 

43 


EE ( I ) = ( YY ( I ) *XX ( 1+1 ) -YY ( 1+1 ) *XX ( I ) ) /EE ( I ) 

44 


IF ( 10— 1 ) 100*100*80 

45 

80 

CONTINUE 

46 


DO 90 1=1*2 

47 

90 

EE ( I ) = ( XX ( 1+2 ) *EE ( I ) -XX ( I ) *EE ( 1+1 ) ) /FF ( I ) 

48 


IF(IEX)100*100*95 

49 

95 

EE ( 1 ) = ( EE ( 1 ) +EE ( 2 ) )/2. 

50 

100 

GIRC=EE(1) 

51 


RETURN 

52 


END 
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1 FUNCTION DTAB2 ( XX , ZI , X , Z , NX , NZ , KX , KZ , Y, Mr IERR ) 

2 C DTAB DOUBLE INTERPOLATION 

3 C XX=VALUE TO BE INTERPOLATED X DIRECTION 

4 C ZI=VALUE TO BE INTERPOLATED Z DIRECTION 

5 C X= X TABLE 

6 C Z= Z TABLE 

7 C NX=NUMBER X VALUES IN TABLE 

8 C NZ=NUMBER Z VALUES IN TABLE 

9 C KX=1 IF LINEAR INTERPOLATION X DIRECTION 


10 

C 

=2 IF CURVILINEAR 

11 

C KZ 

=1 IF LINEAR INTERPOLATION Z DIRECTION 

12 

C 

=2 IF CURVILINEAR 

13 

C Y= 

Y ARRAY 

14 

C M= 

ROW DIMENSION Y ARRAY 

15 

C IERR=1 SUCCESSFUL RETURN 

16 

C 

=2 UNSUCCESSFUL RETURN 

17 


DIMENSION X(l> ,Z(1) ,Y<1> ,X2(4> ,Z2(4> ,Y2(16> ,INX(2) ,IKX<2> , 

18 


1ISV (2) , IXINT (2) 

19 


ZZ=ZI 

20 


INX(1)=NX 

21 


INX(2)=NZ 

22 


IKX(1)=KX 

23 


IKX(2)=KZ 

24 


DO 200 1=1,2 

25 


IF(I NX (I) -2)1000,140, 10 

26 

10 

IF ( 1-1 ) 11 , 11 , 12 

27 

11 

IF ( X ( 1 ) -X (2) )20,1000,30 

28 

c 

ASCENDING ORDER 

29 

20 

DO 21 J=1,NX 

30 


IF(X( J)-XX)21,40,50 

31 

21 

CONTINUE 

32 

C 

VALUE BEYOND END OF TABLE 

33 

22 

ISV(I)=INX(I)-IKX(I) 

34 

23 

IXINT ( I ) =IKX ( I ) 

35 


GO TO 150 

36 

12 

IF ( Z ( 1 ) -Z ( 2 ) ) 25, 1000, 35 

37 

25 

DO 26 J=1,NZ 

38 


IF (Z(J)-ZZ)26,40,50 

39 

26 

CONTINUE 

40 


GO TO 22 

41 

C 

DESCENDING ORDER 

42 

30 

DO 31 J=1 , NX 

43 


IF ( X ( J) — XX >50,40,31 

44 

31 

CONTINUE 

45 


GO TO 22 

46 

35 

DO 36 J=1,NZ 

47 


IF { Z { J) -ZZ ) 50,40,36 

48 

36 

CONTINUE 

49 


GO TO 22 

50 

C 

VALUE EQUAL VALUE IN TABLE 
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DTAB2 


51 

40 

IXINT { I )=0 

52 


ISV(I)=J 

53 


GO TO 150 

54 

50 

IF ( J-l )60r60>70 

55 

60 

ISV ( I ) =1 

56 


GO TO 23 

57 

70 

IF ( J-2 ) 60 r 60 » 80 

58 

80 

ISV ( I ) =J-IKX ( I ) 

59 


IF ( IKX ( I ) -2 ) 23 f 90 f 90 

60 

90 

IF ( ISV ( I ) +2+IKX ( I ) -INX { I ) ) 100 t 100 

61 

100 

IXINT ( I ) =3 

62 


GO TO 150 

63 

140 

IXINT ( I ) =1 

64 


ISV(I)=1 

65 

150 

J1=ISV(I) 

66 


IF ( I-l ) 155 r 155 ? 160 

67 

155 

DO 156 J=1 r 4 

68 


X2( J)=X( Jl) 

69 

156 

J1=J1+1 

70 


GO TO 200 

71 

160 

DO 165 J=1 r 4 

72 


Z2(J)=Z(J1) 

73 

165 

J1=J1+1 

74 

200 

CONTINUE 

75 


IC=1 

76 


K= ( ISV ( 1 ) -1 ) *M+ISV ( 2 ) -M 

77 


DO 220 1=1 >4 

78 


K=K+M 

79 


DO 210 J=1 f 4 

80 


K1=K+ J-l 

81 


Y2(IC)=Y(K1) 

82 

210 

IC=IC+1 

83 

220 

CONTINUE 

84 

225 

IF ( IXINT (2) ) 301 r 301 » 230 

85 

230 

IZO=IXINT (2) 

86 


IZ02=IZ0 

87 


DO 280 K=1 r IZO 

88 


DO 270 1=1 t IZ02 

89 


DZ1=Z2 ( I ) -ZZ 

90 


KI=K+I 

91 


DZ2=Z2 ( K I ) -ZZ 

92 


DZ=DZ2-DZ1 

93 


IXO=IXINT(l) *4+1 

94 


DO 260 J=I t IXOr 4 

95 

260 

Y2 ( J ) = ( Y2 ( J ) *DZ2-Y2 ( J+l ) *DZ1 ) /DZ 

96 

270 

CONTINUE 

97 


IZ02=IZ02-1 

98 

280 

CONTINUE 

99 

301 

IF ( IXINT (1) ) 400 f 400 t 300 

100 

300 

IXINT (2 ) =IXINT ( 1 ) 
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DTAB2 


101 


IXINT (1 ) =0 

102 

302 

IC=1 

103 


DO 310 1=1*4 

104 


Z2 ( I ) =X2 ( I ) 

105 


Y2(I)=Y2(IC) 

106 

310 

IC=IC+4 

107 


zz=xx 

108 


GO TO 225 

109 

400 

DTAB2 =Y2(1) 

110 


IERR=1 

111 


RETURN 

112 

1000 

IERR=2 

113 


DTAB2 =0. 

114 


RETURN 

115 


END 
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SAMPLE CASE 


This sample case is one of the program runs made to generate comparisons 
between the predictions of this program and the noise levels measured in- 
flight on the NASA/Dryden JetStar. The input data used are shown on the 
next two pages, followed by the program output. 
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JETSTAR CRUISE NOISE CHECK CASE 
TP’S 3Af3B»3C MIC LOCATIONS 1 AND 3 
010000010101010000000000 
0000000000 
0.076475 1116.4 

010100 
001 
10.20 
0101 
1.25 

024054 20. 


80. 

315. 

3.281 

080 1 

.455 

0.915 

1116.4 

1.0 


450. 

230.9 

35. 

1400. 

1.079 

00 

1.079 

1000. 


53.34 

003 

32,174 


36000. 

0.70 

1.9 

0.022798 

968.5 

390.3 

16.1 

13911. 


28. 

1.118 

1.044 

3.52 

13911. 


16.1 

1896. 

1292. 

002 

001 

002 

001 

0.0932 

756. 

10.20 

002 

15.38 


10.20 

003 

002 

001 

17.88 


13.80 

002 

10.06 


13.80 

12.56 


36000. 

0.75 

2.0 

0.022798 

968.5 

390.3 

17.2 

14461. 


31. 

1.149 

1.078 

3.68 

14461. 


17.2 

2041. 

1382. 

0.1029 

790. 


1662. 


1789. 
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002 

001 

002 

001 

10.20 

15.38 


002 

10.20 

17.88 


003 

002 

001 

13.80 

10.06 


002 

13.80 

12.56 


36000 . 

0.80 

2.1 

0.022798 

968.5 

390.3 

18.3 

14990 . 


33 , 

1.180 

1.110 

3.81 

14990 . 


18.3 

0.1132 

821 . 

2186 . 

1470 . 

002 

001 

0 Q 2 

001 

10.20 

15.38 


002 

10.20 

17,88 


003 

002 

001 

13.80 

10.06 


002 

13.80 

12.56 



1913 . 



Sample Case Output 
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WARNING 


WARNING 


IN JET MIXING NOISE* ATTEMPT TO EXPAND LOW FREO RANGE TO 50 HZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 


IN JET MIXING NOISE* ATTEMPT TO EXPAND LOW FREO RANGE TO 50 HZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 





L 


JETSTAR CRUISE NOISE CHECK CASE 
TP’S 3Ar3B#3C MIC LOCATIONS 1 AND 3 

FLIGHT CONDITIONS ALT = 36000. MA = .700 POWER SET. CODE = 1.90 

OBSERVER LOCATION NO. = i 
TOTAL ENGINE AND AIRFRAME NOISE 


1/3 OCT. BANO SPL’S AND OASPL’S DB RE 0.00002 N/M2r 


FREO 

(HZ) 

FAN 

FWD 

FAN 

AFT 

FAN 

TOTAL 

CMX 

TURB 

CORE 

JET 

PRI 

MIXING 

SEC 

NOISE 

TOTAL 

SHOCK 

PRI 

BD 

SEC 

SHOCK 

PRI 

SCREECH 

SEC 

AIRFRAME 
TEN TBLN 

TOTAL 

50. 

44.1 

.0 

44.1 

.0 

47.0 

71.4 

74.8 

.0 

74.8 

94.1 

.0 

.0 

• 0 

• 0 

• 0 

94,2 

63 • 

47.1 

.0 

47.1 

.0 

50.3 

75.4 

76.7 

.0 

76.7 

98.0 

.0 

.0 

• 0 

• 0 

• 0 

q 0.1 

60. 

50.2 

.0 

50.2 

.0 

54.7 

79.4 

78.6 

.0 

78.6 

101.7 

.0 

.0 

• 0 

• 0 

• 0 

loi .n 

100* 

53.1 

.0 

53.1 

.0 

58.3 

82.4 

80.5 

.0 

80.5 

104.8 

.0 

.0 

• 0 

• 0 

• 0 

104 .0 

125* 

56.0 

.0 

56.0 

.0 

59.7 

85.4 

82.2 

.0 

82.2 

107.5 

.0 

.0 

• 0 

• 0 

• 0 

107 .6 

160. 

59.2 

.0 

59.2 

.0 

61.6 

88.4 

84.0 

.0 

84.0 

110.4 

.0 

.0 

• 0 

• 0 

• 0 

110.4 

200* 

62.2 

.0 

62.2 

.0 

63.4 

90.9 

85.7 

.0 

85.7 

112.0 

.0 

.0 

• 0 

• 0 

• 0 

112.0 

250* 

65.1 

.0 

65.1 

.0 

65.2 

92.9 

86.9 

.0 

86.9 

112.8 

.0 

.0 

• 0 

• 0 

• 0 

112 . n 

3i5i 

68.2 

.0 

68.2 

.0 

67.5 

94.4 

88.1 

.0 

88.1 

113.7 

.0 

.0 

• 0 

• 0 

• 0 

1 13.0 

400. 

71.4 

.0 

71.4 

.0 

70.0 

95.4 

89.3 

.0 

89.3 

116.9 

.0 

.0 

• 0 

• 0 

• 0 

1 17.0 

500. 

74.5 

.0 

74.5 

.0 

71.9 

94.4 

89.9 

.0 

89.9 

124.3 

.0 

.0 

• 0 

• 0 

• 0 

124.3 

630 . 

77.8 

.0 

77.8 

.0 

73.8 

92.9 

90.5 

.0 

90.5 

131.2 

.0 

.0 

• 0 

• 0 

• 0 

131.2 

000 • 

80.8 

.0 

80.8 

.0 

75.9 

90.9 

91.1 

.0 

91.1 

132.6 

.0 

.0 

• 0 

• 0 

• 0 

1 32 .f» 

1000. 

83.8 

.0 

B3.8 

.0 

78.0 

88.4 

92.2 

.0 

92.2 

130.5 

.0 

.0 

• 0 

• 0 

• 0 

130.5 

1250 . 

87.0 

.0 

87.0 

.0 

79.9 

85.4 

93.3 

.0 

93.3 

130.1 

.0 

.0 

• 0 

• 0 

• 0 

130.1 

1600 • 

89.9 

.0 

89.9 

.0 

82.5 

82.4 

94.5 

.0 

94.5 

129.3 

.0 

.0 

• 0 

• 0 

• 0 

12°. 3 

2000. 

92.8 

.0 

92.8 

.0 

85.0 

79.4 

94.9 

.0 

94.9 

128.6 

.0 

.0 

• 0 

• 0 

• 0 

120.0 

2500 . 

95.7 

.0 

95.7 

.0 

86.9 

75.4 

95,3 

.0 

95.3 

127.3 

.0 

.0 

• 0 

• 0 

• 0 

127.3 

3 1 50 . 

94.0 

.0 

94.0 

.0 

89.1 

71.4 

95.7 

.0 

95.7 

126.2 

.0 

.0 

• 0 

• 0 

• 0 

120.3 

4000. 

93.5 

.0 

93.5 

.0 

91.2 

67.9 

95.3 

.0 

95.3 

125.1 

.0 

.0 

• 0 

• 0 

• 0 

125.1 

123,0 

5000. 

100.4 

.0 

100.4 

.0 

92.8 

63.9 

94.9 

.0 

94.9 

123.8 

.0 

.0 

• 0 

• 0 

• 0 

6300 . 

94.8 

.0 

94.8 

.0 

94.3 

59.4 

94.5 

.0 

94.5 

122.5 

.0 

.0 

• 0 

• 0 

• 0 

122.5 

0000. 

95.6 

.0 

95.6 

.0 

95.7 

55.4 

93.8 

.0 

93.8 

121.2 

.0 

.0 

,0 

• 0 

• 0 

121.2 

10000. 

100.3 

.0 

100.3 

.0 

96.7 

50.4 

93.0 

.0 

93.0 

120.0 

.0 

.0 

• 0 

• 0 

.0 

120.1 

OASPL 

106.1 

.0 

106.1 

.0 

102.1 

102.6 

105.7 

.0 

105.7 

139.6 

.0 

.0 

.0 

• 0 

• 0 

1 3 U .7 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A»3B»3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITIONS ALT = 36000. MA = 
OBSERVER LOCATION NO. = 1 

ENGINE NO. 1 


.700 


POWER SET. CODE 


1.90 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 


*** INLET 

RADIATED NOISE *** 

** EXHAUST DUCT ** 
* RADIATED NOISE * 


FREQ 

(HZ) 

BROAD 

-BAND 

DSCRT 

TONES 

CMBNTN 

TONES 

TOTAL 

BROAD 

-BAND 

DSCRT 

TONES TOTAL 

TOTAL 

FAN 

NOISE 

50. 

.0 

.0 

41.7 

41.7 

.0 

.0 

.0 

41.7 

63. 

.0 

.0 

44.7 

44.7 

.0 

.0 

.0 

44.7 

BO. 

5.9 

.0 

47.8 

47.8 

.0 

.0 

.0 

47.8 

100. 

12.6 

.0 

50.7 

50.7 

.0 

.0 

.0 

50.7 

125. 

19. a 

.0 

53.6 

53.6 

.0 

.0 

.0 

53.6 

160. 

27.4 

.0 

56.8 

56.8 

.0 

.0 

.0 

56.8 

200. 

34.1 

.0 

59.8 

59.8 

.0 

.0 

.0 

59.8 

250. 

40.3 

.0 

62.7 

62.7 

.0 

.0 

.0 

62.7 

315. 

46.5 

.0 

65.8 

65.8 

.0 

.0 

.0 

65.8 

400. 

52.4 

.0 

68.9 

69.0 

.0 

.0 

.0 

69.0 

500. 

57.6 

.0 

72.0 

72.1 

.0 

.0 

.0 

72.1 

630. 

62.6 

.0 

75.1 

75.4 

.0 

.0 

.0 

75.4 

BOO. 

67.4 

.0 

78.0 

78.4 

.0 

.0 

.0 

78.4 

1000. 

71.5 

.0 

81.0 

81.4 

.0 

.0 

.0 

81.4 

1250. 

75.3 

.0 

84.0 

84.5 

.0 

.0 

.0 

84.5 

1600. 

79.0 

.0 

86.9 

87.5 

.0 

.0 

.0 

87.5 

2000. 

82.1 

.0 

89.7 

90.4 

.0 

.0 

.0 

90.4 

2500. 

B4.7 

.0 

92.6 

93.3 

.0 

.0 

.0 

93.3 

3150. 

87.1 

.0 

89.6 

91.6 

.0 

.0 

.0 

91.6 

4000. 

89.2 

.0 

86.5 

91.1 

.0 

.0 

.0 

91.1 

5000. 

90.8 

96.9 

83.6 

98.0 

.0 

.0 

.0 

98.0 

6300. 

92.1 

.0 

80.6 

92.4 

.0 

.0 

.0 

92.4 

BOOO. 

93.1 

.0 

77.5 

93.2 

.0 

.0 

.0 

93.2 

10000. 

93.6 

95.9 

74.6 

97.9 

.0 

.0 

.0 

97.9 

OASPL 

90.6 

99.5 

97.4 

103.7 

.0 

.0 

.0 

103.7 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A>3B»3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION! ALT = 36000. MA = .700 POWER SET. CODE = l.°0 


OBSERVER LOCATION NO. = 1 


ENGINE NO. 2 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 


*** INLET RAOIATED NOISE *** 

** EXHAUST DUCT ** 
* RADIATED NOISE * 

TOTAL 

FAN 

NOISE 

FREO 

(ECZI 

BROAD 

-BAND 

DSCRT 

TONES 

CMBNTN 

TONES 

TOTAL 

BROAD 

-BAND 

DSCRT 

TONES TOTAL 

50. 

.0 

.0 

40.3 

40.3 

.0 

.0 

.0 

40.3 

63 . 

.0 

.0 

43.4 

43.4 

.0 

.0 

.0 

43.4 

60. 

5.0 

.0 

46.5 

46.5 

.0 

.0 

.0 

46.5 

100. 

11.4 

.0 

49.4 

49.4 

.0 

.0 

.0 

49.4 

125. 

16.5 

.0 

52.3 

52.3 

.0 

.0 

.0 

52.3 

160. 

26.1 

.0 

55.5 

55.5 

.0 

.0 

.0 

55.5 

200. 

32.8 

.0 

58.5 

58.5 

.0 

.0 

.0 

58.5 

250. 

39.0 

.0 

61.4 

61.4 

.0 

.0 

.0 

61.4 

315. 

45.2 

.0 

64.5 

64.5 

.0 

.0 

.0 

64.5 

400. 

51.1 

.0 

67.6 

67.7 

.0 

.0 

.0 

67.7 

500. 

56.3 

.0 

70.7 

70.8 

.0 

.0 

.0 

70. B 

630. 

61.3 

.0 

73.8 

74.1 

.0 

.0 

.0 

74.1 

000. 

66.1 

.0 

76.7 

77.1 

.0 

.0 

.0 

77.1 

1000. 

70.2 

.0 

79.7 

80.1 

.0 

.0 

.0 

80.1 

1250. 

74.0 

.0 

82.7 

83.2 

.0 

.0 

.0 

83.2 

1600. 

77.7 

.0 

85.6 

86.2 

.0 

.0 

.0 

B6.2 

2000. 

60.7 

.0 

86.4 

89.1 

.0 

.0 

.0 

B9.1 

2500. 

83.4 

.0 

91.3 

92.0 

.0 

.0 

.0 

92.0 

3150. 

65.6 

.0 

80.3 

90.3 

.0 

.0 

.0 

90.3 

4000, 

87.9 

.0 

85.2 

89.8 

.0 

.0 

.0 

89.8 

5000. 

09.5 

95.6 

82.3 

96.7 

.0 

.0 

.0 

96.7 

6300. 

90.8 

.0 

79.3 

91.1 

.0 

.0 

.0 

91.1 

8000. 

91.8 

.0 

76.2 

91.9 

.0 

.0 

.0 

91 .9 

10000. 

92.3 

94.6 

73.3 

96.6 

.0 

.0 

.0 

96.6 

OASPL 

96.3 

96.1 

96.1 

102.4 

.0 

.0 

.0 

102.4 

WARNING 

- IN 

JET MIXING NOISE » ATTEMPT TO EXPAND LOW 

FREO 

RANGE TO 


EXCEFDED LIMITS OF 

THOMSON 

SPECTRUM SHAPES 




WARNING - IN JET MIXING NOISE t ATTEMPT TO EXPAND LOW FRT 0 PANG! TO 50 117 
EXCEEDED LIMITS OE THOMSON SPECTRUM SHAPES 
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JETSTAR CRUISE NOISE CHECK CASE 
TP»S 3A.3B.3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA = .700 POWER SET. CODE = 1.00 

OBSERVER LOCATION NO. = 3 

TOTAL ENGINE AND AIRFRAME NOISE 


1/3 OCT. BAND SPL'S AND OASPL’S OB RE 0.00002 N/M2» 


FREQ 

FAN 

FAN 

FAN 




JET 

MIXING 

NOISE 

SHOCK 

BB 

SHOCK 

SCREECH 

AIRFRAME 


(HZ) 

FWD 

AFT 

TOTAL 

CMX 

TURB 

CORE 

PR I 

SEC 

TOTAL 

PR I 

SEC 

PR I 

SEC 

TEN 

TBLH 

TOTAL 

50. 

44.9 

.0 

44.9 

.0 

46.7 

71.6 

73.9 

.0 

73.9 

95.1 

.0 

.0 

.0 

.0 

.0 

"5.2 

63. 

47.9 

.0 

47.9 

.0 

50.1 

75.6 

75.8 

.0 

75.8 

99.0 

.0 

.0 

.0 

.0 

.0 

oo.i 

80. 

51.0 

.0 

51.0 

.0 

54.4 

79.6 

77.8 

.0 

77.8 

102.7 

.0 

.0 

.0 

.0 

.0 

102.7 

100. 

53.9 

.0 

53.9 

.0 

58.1 

82.6 

79.6 

.0 

79,6 

105.7 

.0 

.0 

.0 

.0 

.0 

105.8 

125. 

56.8 

.0 

56.8 

.0 

59.4 

85.6 

81.3 

.0 

81.3 

108,4 

.0 

.0 

.0 

.0 

.0 

108.4 

160. 

60.1 

.0 

60.1 

.0 

61.3 

88.6 

83.1 

.0 

83.1 

111.2 

.0 

.0 

.0 

.0 

.0 

111.2 

200. 

63.0 

.0 

63.0 

.0 

63.2 

91.1 

84.8 

.0 

84.8 

112.7 

.0 

.0 

.0 

.0 

.0 

112.7 

250. 

66.0 

.0 

66.0 

.0 

64.9 

93.1 

86.0 

.0 

86.0 

113.5 

.0 

.0 

.0 

.0 

.0 

113.5 

315. 

69.1 

.0 

69.1 

.0 

67.2 

94.6 

87.2 

.0 

87.2 

114.5 

.0 

.0 

.0 

.0 

.0 

114.6 

400. 

72.5 

.0 

72.5 

.0 

69.7 

95.6 

88.4 

.0 

88.4 

118.3 

.0 

.0 

.0 

.0 

.0 

118.4 

500. 

75.7 

.0 

75.7 

.0 

71.6 

94.6 

89.0 

.0 

89.0 

126.1 

.0 

.0 

.0 

.0 

.0 

126.1 

630. 

79.1 

.0 

79.1 

.0 

73.6 

93.1 

89.6 

.0 

89.6 

132.5 

.0 

.0 

.0 

.0 

.0 

132.5 

800. 

82.3 

.0 

82.3 

.0 

75.6 

91.1 

90.3 

.0 

90.3 

133.5 

.0 

.0 

.0 

.0 

.0 

133.5 

1000. 

85.5 

.0 

85.5 

.0 

77.7 

88.6 

91.4 

.0 

91.4 

131.3 

.0 

.0 

.0 

.0 

.0 

131.3 

1250. 

88.8 

.0 

88.8 

.0 

79.6 

85.6 

92.5 

.0 

92.5 

131.3 

.0 

.0 

.0 

.0 

.0 

131.3 

1600. 

91.9 

.0 

91.9 

.0 

82.3 

82.6 

93.6 

.0 

93.6 

130.5 

.0 

.0 

.0 

.0 

.0 

130.5 

2000. 

94.9 

.0 

94.9 

.0 

84.7 

79.6 

94.0 

.0 

94.0 

129.5 

.0 

.0 

.0 

.0 

.0 

12°. 5 

2500. 

97.7 

.0 

97.7 

.0 

86.7 

75.6 

94.4 

.0 

94.4 

128.4 

.0 

.0 

.0 

.0 

.0 

128.4 

3150. 

97.4 

.0 

97.4 

.0 

88.8 

71.6 

94.8 

.0 

94.8 

127.3 

.0 

.0 

.0 

.0 

.0 

127.3 

4000. 

98.2 

.0 

98.2 

.0 

90.9 

68.1 

94.5 

.0 

94.5 

126.1 

.0 

.0 

.0 

.0 

.0 

126.1 

5000. 

105.9 

.0 

105.Q 

.0 

92.5 

64.1 

94.1 

.0 

94.1 

124.8 

.0 

.0 

.0 

.0 

.0 

124.9 

6300. 

100.6 

.0 

100.6 

.0 

94.1 

59.6 

93.7 

.0 

93.7 

123.6 

.0 

.0 

.0 

.0 

.0 

123.6 

8000. 

101.5 

.0 

101.5 

.0 

95.4 

55.6 

92.9 

.0 

92.9 

122.3 

.0 

.0 

.0 

.0 

.0 

122.3 

10000. 

106.0 

.0 

106.0 

.0 

96.5 

50.6 

92.2 

.0 

92.2 

121.0 

.0 

.0 

.0 

.0 

.0 

121.2 

OASPL 

111.1 

.0 

111.1 

.0 

101.9 

102.8 

104.8 

.0 

104.8 

140.7 

.0 

.0 

.0 

.0 

.0 

140.7 
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JETSTAR CRUI5E NOISE CHECK CASE 
TP'S 3A.3B.3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000 


MA = .700 POWER SET. CODE = 


OBSERVER LOCATION NO. = 3 


ENGINE NO. 1 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 

** EXHAUST DUCT ** 

*** INLET RADIATED NOISE *** * RADIATED NOISE * 


FREQ 

BROAD 

DSCRT 

(HZ) 

-BAND 

TONES 

50. 

.0 

.0 

63. 

4.3 

.0 

80. 

11.2 

.0 

100. 

18.6 

.0 

125. 

25.9 

.0 

160. 

33.6 

.0 

200. 

40.3 

.0 

250. 

46.5 

.0 

315. 

52.7 

.0 

400. 

58.6 

.0 

500. 

63.8 

.0 

630. 

68.8 

.0 

800. 

73.6 

.0 

1000. 

77.7 

.0 

1250. 

81.5 

.0 

1600. 

85.2 

.0 

2000. 

88.3 

.0 

2500. 

90.9 

.0 

3150. 

93.3 

.0 

4000. 

95.4 

.0 

5000. 

97.0 

102.7 

6300. 

98.3 

.0 

8000. 

99.3 

.0 

10000. 

99.8 

101.7 

OASPL 

105.8 

105.2 


CMBNTN BROAD 

TONES TOTAL -BAND 

42.4 42.4 .0 

45.5 45.5 .0 

48.6 48.6 .0 

51.5 51.5 .0 

54.4 54.4 .0 

57.6 57.6 .0 

60.5 60.6 .0 

63.5 63.6 .0 

66.5 66.7 .0 

69.7 70.1 .0 

72.7 73.3 .0 

75.9 76.7 .0 

78.8 80.0 .0 

81.8 83.2 .0 

84.8 86.4 .0 

87.7 89.6 .0 

90.5 92.5 .0 

93.4 95.4 .0 

90.4 95.1 .0 

87.3 96,1 .0 

84.4 103.8 .0 

81.4 98.4 .0 

78.3 99.3 .0 

75.3 103.9 .0 

98.2 108.9 .0 


DSCRT 


TOTAL 

FAN 

TONES 

TOTAL 

NOISE 

.0 

.0 

42.4 

.0 

.0 

45.5 

.0 

.0 

48.6 

.0 

.0 

51.5 

.0 

.0 

54.4 

.0 

.0 

57.6 

.0 

.0 

60.6 

.0 

.0 

63.6 

.0 

.0 

66.7 

.0 

.0 

70.1 

.0 

.0 

73.3 

.0 

.0 

76.7 

.0 

.0 

80.0 

.0 

.0 

83.2 

.0 

.0 

86.4 

.0 

.0 

89.6 

.0 

.0 

92.5 

.0 

.0 

95.4 

.0 

.0 

95.1 

.0 

.0 

96.1 

.0 

.0 

103.8 

.0 

.0 

98.4 

.0 

.0 

9'i.3 

.0 

.0 

103.9 

.0 

.0 

108.9 


1.90 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A.3B.3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION! ALT = 36000. MA = .700 POWER SET. CODE = 1.90 

OBSERVER LOCATION NO. = 3 

ENGINE NO. 2 

FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 

** EXHAUST DUCT ** 

*** INLET RADIATED NOISE *** * RADIATED NOISE * 


TOTAL 


FREQ 

BROAD 

OSCRT 

CMBNTN 


BROAD 

DSCRT 


FAN 

(HZ1 

-BAND 

TONES 

TONES 

TOTAL 

-BAND 

TONES 

TOTAL 

NOISE 

50. 

.0 

.0 

41.2 

41.2 

.0 

.0 

.0 

41.2 

63. 

3.1 

.0 

44.2 

44.2 

.0 

.0 

.0 

44.2 

no. 

9.4 

.0 

47.3 

47.3 

.0 

.0 

.0 

47.3 

100. 

16.7 

.0 

50.2 

50.2 

.0 

.0 

.0 

50.2 

125. 

24.0 

.0 

53.1 

53.1 

.0 

.0 

.0 

53.1 

160. 

31.7 

.0 

56.4 

56.4 

.0 

.0 

.0 

56.4 

200. 

38.3 

.0 

59.3 

59.3 

.0 

.0 

.0 

59.3 

250. 

44.6 

.0 

62.2 

62.3 

.0 

.0 

.0 

62.3 

315. 

50.7 

.0 

65.3 

65.4 

.0 

.0 

.0 

65.4 

400. 

56.6 

.0 

68.5 

68.7 

.0 

.0 

.0 

68.7 

500. 

61.8 

.0 

71.5 

71.9 

.0 

.0 

.0 

71.9 

630. 

66.8 

.0 

74.7 

75.3 

.0 

.0 

.0 

75.3 

800. 

71.6 

.0 

77.6 

78.5 

.0 

.0 

.0 

78.5 

1000. 

75.7 

.0 

80.5 

81,7 

.0 

.0 

.0 

81.7 

1250. 

79.5 

.0 

83.5 

85.0 

.0 

.0 

.0 

85.0 

1600. 

83.3 

.0 

86.4 

88.1 

.0 

.0 

.0 

68.1 

2000. 

86.3 

.0 

89.2 

91.0 

.0 

.0 

.0 

91.0 

2500. 

89.0 

.0 

92.1 

93.9 

.0 

.0 

.0 

93.9 

3150. 

91.4 

.0 

89.1 

93.4 

.0 

.0 

.0 

93.4 

4000. 

93.5 

.0 

86.0 

94.2 

.0 

.0 

.0 

94.2 

5000. 

95.1 

100.7 

83.1 

101.8 

.0 

.0 

.0 

101. 8 

6300. 

96.4 

.0 

80.1 

96.5 

.0 

.0 

.0 

96.5 

8000. 

97.3 

.0 

77.0 

97.4 

.0 

.0 

.0 

97.4 

10000. 

97.8 

99.7 

74.1 

101.9 

.0 

.0 

.0 

101.9 

OASPL 

103.8 

103.2 

97.0 

107.0 

.0 

.0 

.0 

107.0 


WARNING - IN JET MIXING NOISE. ATTEMPT TO EXPAND LOW FREO RANGE TO 50 HZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 


WARNING - IN JET MIXING NOISE. ATTEMPT TO EXPAND LOW FRtO RANGE TO 50 117 
FXCEECrO LIMITS OF THOMSON SPECTRUM SIIAPFS 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A»3B»3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA = .750 POWER SET. CODE = £>.00 

OBSERVER LOCATION NO. = 1 

TOTAL ENGINE AND AIRFRAME NOISE 


1/3 OCT. BAND SPL'S AND OASPL'S DB RE 0.00002 N/M2r 


FREQ 

(HZ) 

FAN 

FWD 

FAN 

AFT 

FAN 

TOTAL 

CMX 

TURB 

CORE 

JET 
PR I 

MIXING 

SEC 

NOISE 

TOTAL 

SHOCK 
PR I 

BB 

SEC 

SHOCK 
PR I 

SCREECH 

SEC 

AIRFRAME 
TEN TBLM 

TOTAL 

50. 

51.8 

.0 

51.8 

.0 

48.4 

75.3 

75.4 

.0 

75.4 

101.6 

.0 

.0 

.0 

.0 

.0 

101 .6 

63 . 

54.8 

.0 

54.8 

.0 

51.6 

79.3 

77.3 

.0 

77.3 

105.4 

.0 

.0 

.0 

.0 

.0 

105,4 

80. 

57.9 

.0 

57.9 

.0 

55.8 

83.3 

79.2 

.0 

79.2 

108.9 

.0 

.0 

.0 

• 0 

• 0 

100,° 

100. 

60.8 

.0 

60.8 

.0 

60.0 

86.3 

81.1 

.0 

81.1 

111.8 

.0 

.0 

.0 

• 0 

• 0 

ii j ,n 

125 . 

63.7 

.0 

63.7 

.0 

61,3 

89.3 

82.8 

.0 

82. B 

114.5 

.0 

.0 

.0 

• 0 

• 0 

114.5 
1 1 (y , q 

160 . 

66.9 

.0 

66.9 

.0 

63.1 

92.3 

84.7 

.0 

84.7 

116.9 

.0 

.0 

.0 

.0 

• 0 

200 . 

69.8 

.0 

69.8 

.0 

65.0 

94.8 

86.3 

.0 

86.3 

117.9 

.0 

.0 

.0 

• 0 

• 0 

lin.o 

250. 

72.8 

.0 

72.8 

.0 

66.7 

96.8 

87.5 

.0 

87.5 

118.6 

.0 

.0 

• 0 

• 0 

• 0 

110.7 

315. 

75.8 

.0 

75.8 

.0 

68.9 

98.3 

88.8 

.0 

88.8 

120.3 

.0 

.0 

.0 

.0 

• 0 

120,3 

400. 

78.9 

.0 

78.9 

.0 

71.6 

99.3 

90.0 

.0 

90.0 

125.8 

.0 

.0 

.0 

• 0 

• 0 

125,0 

500. 

81.8 

.0 

81.8 

.0 

73.4 

98.3 

90.5 

.0 

90.5 

133.6 

.0 

.0 

.0 

• 0 

• 0 

133.0 

630 • 

84.9 

.0 

84.9 

.0 

75.4 

96.8 

90.9 

.0 

90.9 

137.8 

.0 

.0 

.0 

.0 

• 0 

137,0 

800. 

88.1 

.0 

80.1 

.0 

77.4 

94.8 

91.4 

.0 

91.4 

136.9 

.0 

.0 

.0 

• 0 

• 0 

136. Q 

1000. 

91.0 

.0 

91.0 

.0 

79.6 

92.3 

92.6 

.0 

92.6 

135.1 

.0 

.0 

.0 

• 0 

• 0 

135.1 
135 • 0 

1250 . 

93.9 

.0 

93.9 

.0 

81.4 

89.3 

93.7 

.0 

93.7 

135.0 

.0 

.0 

• 0 

.0 

• 0 

1600 . 

97.2 

.0 

97.2 

.0 

83.9 

86.3 

95.0 

.0 

95.0 

134.4 

.0 

.0 

.0 

• 0 

• 0 
• 0 

134.4 

132.0 

132.0 

2000. 

100.0 

.0 

100.0 

.0 

86.6 

83.3 

95.4 

.0 

95.4 

132.8 

.0 

.0 

.0 

• 0 

2500 . 

102.9 

.0 

102.9 

.0 

88.4 

79.3 

95.9 

.0 

95.9 

132.0 

.0 

.0 

.0 

• 0 

• 0 

3150. 

105.9 

.0 

105.9 

.0 

90.6 

75.3 

96,4 

.0 

96.4 

130.8 

.0 

.0 

.0 

• 0 

• 0 

130.0 

4000. 

103. 1 

.0 

103.1 

.0 

92.7 

71.8 

96.1 

.0 

96.1 

129.5 

.0 

.0 

.0 

• 0 

• 0 

12”. 5 

5000. 

101.0 

.0 

101.0 

.0 

94.4 

67.8 

95.8 

.0 

95.8 

128.2 

.0 

.0 

.0 

• 0 

• 0 

120.2 

6300 . 

104.1 

.0 

104.1 

.0 

95.9 

63.3 

95.4 

.0 

95.4 

126.9 

.0 

.0 

.0 

, 0 

# 0 

126.° 

8000. 

9 q .5 

.0 

99.5 

.0 

97.3 

59.3 

94.7 

.0 

94.7 

125.7 

.0 

.0 

.0 

.0 

• 0 

125.7 

loono • 

99.8 

.0 

99.8 

.0 

98.4 

54.3 

94.0 

.0 

94.0 

124.5 

.0 

.0 

.0 

.0 

.0 

124.5 

OASPL 

111.9 

.0 

111.9 

.0 

103.0 

106.4 

106.3 

.0 

106.3 

144.8 

.0 

.0 

.0 

.0 

.0 

144.0 
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JETSTAR CRUISE NOISE CHECK CASE 
TP»S 3A.3B.3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION! ALT = 36000. MA = .750 POWER SET. CODE = 2.00 


OBSERVER LOCATION NO. = 1 


ENGINE NO. 1 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL*S AND OASPL'S - DB RE 0.00002 N/M2 


** EXHAUST DUCT ** 

**« INLET RADIATED NOISE *** * RADIATED NOISE * 

TOTAL 


FREO 

BROAD 

DSCRT 

CMBNTN 


BROAD 

DSCRT 


FAN 

(HZ) 

-BAND 

TONES 

TONES 

TOTAL 

-BAND 

TONES 

TOTAL 

NOISE 

50. 

.0 

.0 

49.4 

49.4 

.0 

.0 

.0 

49.4 

63 . 

.0 

.0 

52.4 

52.4 

.0 

.0 

.0 

52.4 

80. 

.0 

.0 

55.5 

55.5 

.0 

.0 

.0 

55.5 

100. 

8.3 

.0 

58.4 

58.4 

.0 

.0 

.0 

58.4 

125. 

15.4 

.0 

61.3 

61.3 

.0 

.0 

.0 

61.3 

160. 

23.5 

.0 

64.5 

64.5 

.0 

.0 

.0 

64.5 

200. 

30.5 

.0 

67.4 

67.4 

.0 

.0 

.0 

67.4 

250. 

37.1 

.0 

70.3 

70.3 

.0 

.0 

.0 

70.3 

315. 

43.7 

.0 

73.4 

73.4 

.0 

.0 

.0 

73.4 

400. 

50.0 

.0 

76.5 

76.5 

.0 

.0 

.0 

76.5 

500. 

55.6 

.0 

79.4 

79.4 

.0 

.0 

.0 

79.4 

630. 

61.0 

.0 

B2.5 

82.5 

.0 

.0 

.0 

82.5 

800. 

66.2 

.0 

85.6 

85.7 

.0 

.0 

.0 

85.7 

1000. 

70.7 

.0 

88.5 

88.6 

.0 

.0 

.0 

88.6 

1250. 

74.8 

.0 

91.4 

91.5 

.0 

.0 

.0 

91.5 

1600. 

79.0 

.0 

94.7 

94.8 

.0 

.0 

.0 

94.8 

2000. 

82.4 

.0 

97.5 

97.6 

.0 

.0 

.0 

97.6 

2500. 

85.5 

.0 

100.3 

100.5 

.0 

.0 

.0 

100.5 

3150. 

88.3 

.0 

103.3 

103.5 

.0 

.0 

.0 

103.5 

4000. 

90.8 

.0 

100.2 

100.7 

.0 

.0 

.0 

100,7 

5000. 

92.8 

.0 

97.3 

98.6 

.0 

.0 

.0 

98.6 

6300. 

94.5 

99.6 

94.3 

101.7 

.0 

.0 

.0 

101.7 

8000. 

95.9 

.0 

91.2 

97.1 

.0 

.0 

.0 

97.1 

10000. 

96.8 

.0 

88.3 

97.4 

.0 

.0 

.0 

97.4 

OASPL 

102.0 

9«.6 

108.1 

109.5 

.0 

.0 

.0 

109.5 
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JETSTAR CRUISE NOISE CHECK CASE 
TP’S 3A*3B*3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITIONS ALT = 36000. MA = .750 POWER SET. CODE 


OBSERVER LOCATION NO. = 1 


2.00 


ENGINE NO. 2 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL’S - DB RE 0.00002 N/M2 







** EXHAUST DUCT 

** 


INLET 

RADIATED NOISE *** 

* RADIATED NOISE * 


FREQ 

BROAD 

DSCRT 

CMBNTN 


BROAD 

DSCRT 


TOTAL 

FAN 

(HZ I 

-BAND 

TONES 

TONES 

TOTAL 

-BAND 

TONES TOTAL 

NOISE 

50. 

.0 

.0 

48.1 

48.1 

.0 

.0 

.0 

48.1 

63 . 

.0 

.0 

51.1 

51.1 

.0 

.0 

.0 

51.1 

60. 

.0 

.0 

54.2 

54.2 

.0 

.0 

.0 

54.2 

100. 

7.2 

.0 

57.1 

57.1 

.0 

.0 

.0 

57.1 

125. 

14.2 

.0 

60.0 

60.0 

.0 

.0 

.0 

60.0 

160. 

22.2 

.0 

63.2 

63.2 

.0 

.0 

.0 

63.2 

200. 

29.2 

.0 

66.1 

66.1 

.0 

.0 

.0 

66.1 

250. 

35.8 

.0 

69.0 

69.0 

.0 

.0 

.0 

69.0 

315. 

42.4 

.0 

72.1 

72.1 

.0 

.0 

.0 

72.1 

400. 

48.7 

.0 

75.2 

75.2 

.0 

.0 

.0 

75.2 

500. 

54.3 

.0 

78.1 

78.1 

.0 

.0 

.0 

78.1 

630. 

59.7 

.0 

81.2 

81.2 

.0 

.0 

.0 

81.2 

800. 

64.9 

.0 

84.3 

84.4 

.0 

.0 

.0 

84.4 

1000. 

69.4 

.0 

87.2 

87.3 

.0 

.0 

.0 

87.3 

1250. 

73.5 

.0 

90.1 

90.2 

.0 

.0 

.0 

90.2 

1600. 

77.7 

.0 

93.4 

93.5 

.0 

.0 

.0 

93.5 

2000. 

81.1 

.0 

96.1 

96.3 

.0 

.0 

.0 

96.3 

2500. 

84.2 

.0 

99.0 

99.2 

.0 

.0 

.0 

99.2 

3150. 

87.0 

.0 

102.0 

102.2 

.0 

.0 

.0 

102.2 

4000. 

89.5 

.0 

98.9 

99.4 

.0 

.0 

.0 

9°. 4 

5000. 

°1 .5 

.0 

96.0 

97.3 

.0 

.0 

.0 

97.3 

6300. 

93.2 

98.3 

93.0 

100.4 

.0 

.0 

.0 

100.4 

6000. 

94.6 

.0 

89.9 

95.8 

.0 

.0 

.0 

95.8 

10000. 

95.5 

.0 

87.0 

96.0 

.0 

.0 

.0 

96.0 

OASPL 

100.7 

98.3 

106.8 

108.2 

.0 

.0 

.0 

108.2 


WARNING - IN JET MIXING NOISE* ATTEMPT TO EXPAND LOW FREQ RANGE TO 50 IIZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 


WARNING 


IN JET MIXING NOISE* ATI E MPT TO EXPAND LOW rur.O RANGE TO 50 M7 
EXCEEDED LIMITS Or THOMSON SI'ErTIMIM SHAPES 

1 1 it 



JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A,3B,3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA = .750 POWER SET. CODE = 2.01) 

OBSERVER LOCATION NO. = 3 

TOTAL ENGINE AND AIRFRAME NOISE 


1/3 OCT. BAND SPL'S AND OASPL'S DB RE 0.00002 N/M2» 


FREQ 

FAN 

FAN 

FAN 




JET 

MIXING 

NOISE 

SHOCK 

BB 

SHOCK 

SCREECH 

AIRFRAME 


(HZ) 

FWD 

AFT 

TOTAL 

CMX 

TURB 

CORE 

PR I 

SEC 

TOTAL 

PR I 

SEC 

PR I 

SEC 

TEN 

TBLN 

TOTAL 

50. 

52.7 

.0 

52.7 

.0 

48.0 

75.3 

74.4 

.0 

74.4 

102.4 

.0 

.0 

.0 

.0 

.0 

102.4 

63 . 

55.7 

.0 

55.7 

.0 

51.2 

79.3 

76.3 

.0 

76,3 

106.2 

.0 

.0 

.0 

.0 

.0 

106.2 

60. 

58.9 

.0 

58.9 

.0 

55.4 

83.3 

78.2 

.0 

78.2 

109.6 

.0 

.0 

.0 

.0 

.0 

109.6 

100. 

61.8 

.0 

61.8 

.0 

59.6 

86.3 

80.1 

.0 

80.1 

112.5 

.0 

.0 

.0 

.0 

.0 

112.5 

125. 

64.7 

.0 

64.7 

.0 

60.9 

89.3 

81.8 

.0 

81.8 

115.2 

.0 

.0 

.0 

.0 

.0 

115.2 

160. 

67.9 

.0 

67.9 

.0 

62.7 

92.3 

83.7 

.0 

83.7 

117.5 

.0 

.0 

.0 

.0 

.0 

117.5 

200. 

70.8 

.0 

70.8 

.0 

64.7 

94.8 

85.3 

.0 

85.3 

118.5 

.0 

.0 

.0 

.0 

.0 

118.5 

250. 

73.7 

.0 

73.7 

.0 

66.3 

96.8 

86.5 

.0 

86.5 

119.2 

.0 

.0 

.0 

.0 

.0 

119.2 

315. 

76.7 

.0 

76.7 

.0 

68.5 

98.3 

87.8 

.0 

87.8 

121.1 

.0 

.0 

.0 

.0 

.0 

121.1 

400. 

79.9 

.0 

79.9 

.0 

71.2 

99.3 

89.0 

.0 

89.0 

127.1 

.0 

.0 

.0 

.0 

.0 

127.1 

500. 

82.8 

.0 

82.8 

.0 

73.0 

98.3 

89.5 

.0 

89.5 

134.9 

.0 

.0 

.0 

.0 

.0 

134." 

630. 

85.9 

.0 

85.9 

.0 

75.0 

96.8 

89.9 

.0 

89.9 

138.6 

.0 

.0 

.0 

.0 

.0 

138.6 

600. 

89.1 

.0 

89.1 

.0 

77.0 

94.8 

90.4 

.0 

90.4 

137.5 

.0 

.0 

.0 

.0 

.0 

137. r» 

1000. 

92.0 

.0 

92.0 

.0 

79.2 

92.3 

91.6 

.0 

91.6 

135.9 

.0 

.0 

.0 

.0 

.0 

135." 

1250. 

95.0 

.0 

95.0 

.0 

81.0 

89.3 

92.7 

.0 

92.7 

135.8 

.0 

.0 

.0 

.0 

.0 

135.8 

1600. 

98.4 

.0 

98.4 

.0 

83.5 

86.3 

94,0 

.0 

94.0 

135.2 

.0 

.0 

.0 

.0 

.0 

135.2 

2000. 

101.2 

.0 

101.2 

.0 

86.2 

83.3 

94.4 

.0 

94.4 

133.6 

.0 

.0 

.0 

.0 

.0 

133.6 

2500. 

104.1 

.0 

104.1 

.0 

88.1 

79.3 

94.9 

.0 

94.9 

132.8 

.0 

.0 

.0 

.0 

.0 

132.8 

3150. 

107.1 

.0 

107.1 

.0 

90.2 

75.3 

95.4 

.0 

95.4 

131.6 

.0 

.0 

.0 

.0 

.0 

131.6 

4000. 

104.8 

.0 

104.8 

.0 

92.4 

71.8 

95.1 

.0 

95.1 

130.3 

.0 

.0 

.0 

.0 

.0 

130.3 

5000. 

103.6 

.0 

103.6 

.0 

94.0 

67.8 

94.8 

.0 

94.8 

129.0 

.0 

.0 

.0 

.0 

.0 

12". 0 

6300. 

106.7 

.0 

108.7 

.0 

95.6 

63.3 

94.4 

.0 

94.4 

127.7 

.0 

.0 

.0 

.0 

.0 

127.8 

8000. 

104.0 

.0 

104.0 

.0 

96.9 

59.3 

93.7 

.0 

93.7 

126.5 

.0 

.0 

.0 

.0 

.0 

126.5 

10000. 

104.6 

.0 

104.6 

.0 

98.1 

54.3 

93.0 

.0 

93.0 

125.3 

.0 

.0 

.0 

.0 

.0 

125.3 

OASPL 

114.5 

.0 

114.5 

.0 

103.4 

106.5 

105.3 

.0 

105.3 

145.7 

.0 

.0 

.0 

.0 

.0 

145.7 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A»3B>3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA = .750 POWER SET. CODE = 2.00 


OBSERVER LOCATION NO. = 3 


ENGINE NO. 1 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 


*** INLET RADIATED NOISE *** 


FREO 

BROAD 

DSCRT 

CMBNTN 


(HZ ) 

-BAND 

TONES 

TONES 

TOTAL 

50. 

.0 

.0 

50.4 

50.4 

63. 

.0 

.0 

53.4 

53.4 

80. 

6.2 

.0 

56.5 

56.5 

100. 

13.3 

.0 

59.4 

59.4 

125. 

20.9 

.0 

62.3 

62.3 

160. 

29.0 

.0 

65.5 

65.5 

200. 

36.0 

.0 

68.5 

68.5 

250. 

42.6 

.0 

71.4 

71.4 

315. 

49.2 

.0 

74.4 

74.4 

400. 

55.5 

.0 

77.5 

77.5 

500. 

61.1 

.0 

80.4 

80.5 

630. 

66.5 

.0 

83.5 

83.6 

800. 

71.7 

.0 

86.6 

86.8 

1000. 

76.2 

.0 

89.5 

89.7 

1250. 

80.3 

.0 

92.4 

92.7 

1600. 

84.5 

.0 

95.7 

96.0 

2000. 

87.9 

.0 

98.5 

98.8 

2500. 

91.0 

.0 

101.4 

101.7 

3150. 

93.8 

.0 

104.4 

104.7 

4000. 

96.3 

.0 

101.3 

102.5 

5000. 

98.3 

.0 

98.3 

101.3 

6300. 

100.0 

105.0 

95.3 

106.6 

8000. 

101.4 

.0 

92.2 

101.9 

10000. 

102.3 

.0 

89.3 

102.5 

OASPL 

107.5 

105.0 

10°. 1 

112.3 


** EXHAUST DUCT 

** 


* RADIATED NOISE 

* 

TOTAL 

BROAD 

DSCRT 


FAN 

-BAND 

TONES TOTAL 

NOISE 

.0 

.0 

.0 

50.4 

.0 

.0 

.0 

53.4 

.0 

.0 

.0 

56.5 

.0 

.0 

.0 

59.4 

.0 

.0 

.0 

62.3 

.0 

.0 

.0 

65.5 

.0 

.0 

.0 

68.5 

.0 

.0 

.0 

71.4 

.0 

.0 

.0 

74.4 

.0 

.0 

.0 

77.5 

.0 

.0 

.0 

80.5 

.0 

.0 

.0 

83.6 

.0 

.0 

.0 

86.8 

.0 

.0 

.0 

69.7 

.0 

.0 

.0 

92.7 

.0 

.0 

.0 

96.0 

.0 

.0 

.0 

98.8 

.0 

.0 

.0 

101.7 

.0 

.0 

.0 

104.7 

.0 

.0 

.0 

102.5 

.0 

.0 

.0 

101.3 

.0 

.0 

.0 

106.6 

.0 

.0 

.0 

101.9 

.0 

.0 

.0 

102.5 

.0 

.0 

.0 

112.3 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A,3B*3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA = .750 POWER SET. CODE = 2.00 

OBSERVER LOCATION MO. = 3 

ENGINE NO. 2 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 


*** INLET 

RADIATED NOISE *** 

** EXHAUST DUCT 
* RADIATED NOISE 

** 
: * 


FREQ 

(HZ) 

BROAD 

-BAND 

DSCRT 

TONES 

CMBNTN 

TONES 

TOTAL 

BROAD 

-BAND 

DSCRT 

TONES TOTAL 

TOTAL 

FAN 

NOISE 

50. 

.0 

.0 

48.9 

4B.9 

.0 

.0 

.0 

48.9 

63. 

.0 

.0 

51.9 

51.9 

.0 

.0 

.0 

51.9 

BO. 

4.B 

.0 

55.1 

55.1 

.0 

.0 

.0 

55.1 

100. 

11.5 

.0 

5B.0 

5B.0 

.0 

.0 

.0 

58.0 

125. 

IB. 9 

.0 

60.9 

60.9 

.0 

.0 

.0 

60.9 

160. 

27.0 

.0 

64.1 

64.1 

.0 

.0 

.0 

64.1 

200. 

34.0 

.0 

67.0 

67.0 

.0 

.0 

.0 

67.0 

250. 

40.7 

.0 

69.9 

69.9 

.0 

.0 

.0 

69.9 

315. 

47.2 

.0 

72.9 

72.9 

.0 

.0 

.0 

72.9 

400. 

53.6 

.0 

76.1 

76.1 

.0 

.0 

.0 

76.1 

500. 

59.2 

.0 

79.0 

79.0 

.0 

.0 

.0 

79.0 

630. 

64.6 

.0 

82.0 

82.1 

.0 

.0 

.0 

82.1 

800. 

69. B 

.0 

05. 2 

85.3 

.0 

.0 

.0 

85.3 

1000. 

74.3 

.0 

Bfl.l 

88.2 

.0 

.0 

.0 

BB.2 

1250. 

7B.4 

.0 

91.0 

91.2 

.0 

.0 

.0 

91.2 

1600. 

02.6 

.0 

94.3 

94.5 

.0 

.0 

.0 

94.5 

2000. 

06. 0 

.0 

97.0 

97.4 

.0 

.0 

.0 

97.4 

2500. 

B9.1 

.0 

99.9 

100.3 

.0 

.0 

.0 

100.3 

3150. 

91.9 

.0 

102.9 

103.2 

.0 

.0 

.0 

103.2 

4000. 

94.4 

.0 

9°,0 

100.9 

.0 

.0 

.0 

100.9 

5000. 

96.4 

.0 

96.9 

99.7 

.0 

.0 

.0 

99.7 

6300. 

9B.1 

103.1 

93.9 

104.7 

.0 

.0 

.0 

104.7 

0000. 

9«.4 

.0 

90. B 

100.0 

.0 

.0 

.0 

100.0 

10000. 

100.3 

.0 

B7.9 

100.6 

.0 

.0 

.0 

100.6 

OASPL 

105.6 

103.1 

107.7 

110.6 

.0 

.0 

.0 

110.6 


WARNING - IN JET MIXING NOISE* ATTEMPT TO EXPAND LOW FREO RANGE TO 50 HZ 
EXCEFDED LIMITS OF THOMSON SPECTRUM SHAPES 


WARNING - FOR SHOCK BROADBAND NOISE* SOURCE = 7 *THE STROUIIAL PARAMETER 
OF 1/3 O.B. NO. 24 IS OUTS IDF R At IGF OF MASTER SITCTRA 
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WARNING - IN JET MIXING NOISE# ATTEMPT TO EXPAND LOW FREQ RANGE TO 50 HZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 


WARNING 


FOR SHOCK BROADBAND NOISE. SOURCE = 7 #THE STROUIIAL PARAMETER 
OF 1/3 O.B. NO. 24 IS OUTSIDE RANGE OF MASTER SPECTRA 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3Ar 3B.3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITIONS ALT = 36000. MA s .800 POWER SET. CODE = 2.10 
OBSERVER LOCATION NO. = 1 

TOTAL ENGINE AND AIRFRAME NOISE 


1/3 OCT. BAND SPL'S AND OASPL'S DB RE 0.00002 N/M2» 


FREQ 

FAN 

FAN 

FAN 




JET 

MIXING 

NOISE 

SHOCK 

BB 

SHOCK 

SCREECH 

AIRFRAME 


(HZ) 

FWD 

AFT 

TOTAL 

CMX 

TURB 

CORE 

PR I 

SEC 

TOTAL 

PR I 

SEC 

PR I 

SEC 

TEN 

TBLN 

TOTAL 

50. 

53.3 

.0 

53.3 

.0 

50.1 

79.5 

75.7 

.0 

75.7 

10B.9 

.0 

.0 

.0 

.0 

.0 

10". 0 

63. 

56.3 

.0 

56.3 

.0 

53.2 

83.5 

77,6 

.0 

77.6 

112.5 

.0 

.0 

.0 

.0 

.0 

112.5 

80. 

59.5 

.0 

59.5 

.0 

57.3 

87.5 

79.6 

.0 

79.6 

115.9 

.0 

.0 

.0 

.0 

.0 

115." 

100. 

62.4 

.0 

62.4 

.0 

62.0 

90.5 

81.4 

.0 

81.4 

118.7 

.0 

.0 

.0 

.0 

.0 

110.7 

125. 

65.3 

.0 

65.3 

.0 

63.3 

93.5 

83.2 

.0 

83.2 

121.3 

.0 

.0 

.0 

.0 

.0 

121.3 

160. 

68.5 

.0 

68.5 

.0 

65.0 

96.5 

85.1 

.0 

85.1 

123.2 

.0 

.0 

.0 

.0 

.0 

123.2 

200. 

71.4 

.0 

71.4 

.0 

67.0 

99.0 

86.8 

.0 

86.8 

124.0 

.0 

.0 

.0 

.0 

.0 

124.0 

250. 

74.3 

.0 

74.3 

.0 

68.6 

101.0 

88.0 

.0 

88.0 

124.7 

.0 

.0 

.0 

.0 

.0 

124. H 

315. 

77.4 

.0 

77.4 

.0 

70.8 

102.5 

89.2 

.0 

89,2 

127.4 

.0 

.0 

.0 

.0 

.0 

127.4 

400. 

80.5 

.0 

80.5 

.0 

73.5 

103.5 

90.5 

.0 

90.5 

134.9 

.0 

.0 

.0 

.0 

.0 

134." 

500. 

83.5 

.0 

83.5 

.0 

75.3 

102.5 

90.8 

.0 

90.8 

141.6 

.0 

.0 

.0 

.0 

.0 

141.6 

630. 

86.6 

.0 

86.6 

.0 

77.3 

101.0 

91.2 

.0 

91.2 

143.4 

.0 

.0 

.0 

.0 

.0 

143.4 

BOO. 

89.8 

.0 

89.8 

.0 

79.3 

99.0 

91.6 

.0 

91.6 

141.2 

.0 

.0 

.0 

.0 

.0 

141.2 

1000. 

92.7 

.0 

92.7 

.0 

81.5 

96.5 

92.7 

.0 

92.7 

140.6 

.0 

.0 

.0 

.0 

.0 

140.6 

1250. 

95.8 

.0 

95.8 

.0 

83.3 

93.5 

93.9 

.0 

93.9 

139.9 

.0 

.0 

.0 

.0 

.0 

13"." 

1600. 

99.2 

.0 

99.2 

.0 

85.7 

90.5 

95.2 

.0 

95.2 

139.2 

.0 

.0 

.0 

.0 

.0 

13". 2 

2000. 

101.7 

.0 

101.7 

.0 

80.5 

87.5 

95.8 

.0 

95.8 

137. B 

.0 

.0 

.0 

.0 

.0 

137.0 

2500. 

104.5 

.0 

104.5 

.0 

90.3 

83.5 

96.3 

.0 

96.3 

136.8 

.0 

.0 

.0 

.0 

.0 

136.8 

3150. 

107.5 

.0 

107.5 

.0 

92.5 

79.5 

96.9 

.0 

96.9 

135.7 

.0 

.0 

.0 

.0 

.0 

135.7 

4000. 

104.9 

.0 

104.9 

.0 

94.7 

76.0 

96.6 

.0 

96.6 

134.3 

.0 

.0 

.0 

.0 

.0 

134.3 

5000. 

103.2 

.0 

103.2 

.0 

96.4 

72.0 

96.4 

.0 

96.4 

133.0 

.0 

.0 

.0 

.0 

.0 

133.0 

6300. 

106.7 

.0 

106.7 

.0 

97.9 

67.5 

*96.1 

.0 

96.1 

131.8 

.0 

.0 

.0 

.0 

.0 

131.0 

8000. 

102.6 

.0 

102.6 

.0 

99.3 

63.5 

95.4 

.0 

95.4 

130.5 

.0 

.0 

.0 

.0 

.0 

130.5 

10000. 

103.0 

.0 

103.0 

.0 

100.5 

58.5 

94.8 

.0 

94.8 

129.5 

.0 

.0 

.0 

.0 

.0 

12". 5 

OASPL 

114.0 

.0 

114.0 

.0 

105.8 

110.6 

106.8 

.0 

106.8 

150.3 

.0 

.0 

.0 

.0 

.0 

150.3 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A>3B»3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA 

OBSERVER LOCATION NO. = 1 

ENGINE NO. 1 


.BOO POWER SFT . CODF = 2.10 


FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL'S - DB RE 0.00002 N/M2 


*** INLET RADIATED NOISE *** 


FREQ 

BROAD 

DSCRT 

CMBNTN 


(HZ) 

-BAND 

TONES 

TONES 

TOTAL 

50. 

.0 

.0 

50.9 

50.9 

63. 

.0 

.0 

53.9 

53.9 

80. 

4.7 

.0 

57.0 

57.0 

100. 

11.3 

.0 

60.0 

60.0 

125. 

18.8 

.0 

62.9 

62.9 

160. 

26.9 

.0 

66.1 

66.1 

200. 

33.9 

.0 

69.0 

69.0 

250. 

40.5 

.0 

71.9 

71.9 

315. 

47.1 

.0 

75.0 

75.0 

400. 

53.4 

.0 

78.1 

70.1 

500. 

59.0 

.0 

81.1 

81.1 

630. 

64.4 

.0 

84.1 

84.2 

800. 

69.6 

.0 

87.4 

87.4 

1000. 

74.1 

.0 

90.2 

90.3 

1250. 

78.2 

.0 

93.2 

93.4 

1600. 

02. 4 

.0 

96.6 

96.0 

2000. 

85.8 

.0 

99.1 

90.3 

2500. 

88.9 

.0 

101.9 

102.1 

3150. 

91.7 

.0 

104.9 

105.1 

4000. 

94.2 

.0 

101.8 

102.5 

5000. 

96.2 

.0 

90.9 

100.8 

6300. 

97.9 

102.2 

95.9 

104.3 

8000. 

9°. 3 

.0 

92.8 

100.1 

10000. 

100.2 

.0 

89.9 

100.6 

OASPL 

105.4 

102.2 

109.7 

111.6 


** EXHAUST DUCT ** 
* RADIATED NOISE * 


BROAD 

DSCRT 


TOTAL 

FAN 

-BAND 

TONES 

TOTAL 

NOISE 

.0 

.0 

.0 

50.9 

.0 

.0 

.0 

53.9 

.0 

.0 

.0 

57.0 

.0 

.0 

.0 

60.0 

.0 

.0 

.0 

62.9 

.0 

.0 

.0 

66.1 

.0 

.0 

.0 

69.0 

.0 

.0 

.0 

71.9 

.0 

.0 

.0 

75.0 

.0 

.0 

.0 

78.1 

.0 . 

.0 

.0 

81.1 

.0 

.0 

.0 

84.2 

.0 

.0 

.0 

87.4 

.0 

.0 

.0 

90.3 

.0 

.0 

.0 

93.4 

.0 

.0 

.0 

96.8 

.0 

.0 

.0 

90.3 

.0 

.0 

.0 

102.1 

.0 

.0 

.0 

105.1 

.0 

.0 

.0 

102.5 

.0 

.0 

.0 

100.8 

.0 

.0 

.0 

104.3 

.0 

.0 

.0 

100.1 

.0 

.0 

.0 

100.6 

.0 

.0 

.0 

111.6 
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JETSTAR CRUISE NOISE CHECK CASE 
TP'S 3A*3B*3C MIC LOCATIONS I AND 3 


FLIGHT CONDITION* ALT = 36000. MA = .BOO POWER SET. CODE = 2.10 
OBSERVER LOCATION NO. = 1 

ENGINE NO. 2 

FAN NOISE COMPONENTS 

1/3 OCT. BANO SPL'S AND OASPL'S - OB RE 0.00002 N/M2 


** EXHAUST DUCT ** 
*** INLET RADIATED NOISE *♦* * RADIATED NOISE * 


FREQ 

(HZ) 

BROAD 

-BAND 

DSCRT 

TONES 

CMBNTN 

TONES 

TOTAL 

BROAD 

-BAND 

DSCRT 

TONES 

TOTAL 

TOTAL 

FAN 

NOISE 

50. 

.0 

.0 

49.6 

49.6 

.0 

.0 

.0 

49.6 

63. 

.0 

.0 

52.6 

52.6 

.0 

.0 

.0 

52.6 

BO. 

3.9 

.0 

55.7 

55.7 

' .0 

.0 

.0 

55.7 

100. 

10.1 

.0 

58.6 

58.6 

.0 

.0 

.0 

58.6 

125. 

17.5 

.0 

61.6 

61.6 

.0 

.0 

.0 

61.6 

160. 

25.6 

.0 

64.8 

64.8 

.0 

.0 

.0 

64. B 

200. 

32.6 

.0 

67.7 

67.7 

.0 

.0 

.0 

67.7 

250. 

39.2 

.0 

70.6 

70.6 

.0 

.0 

.0 

70.6 

315. 

45.8 

.0 

73.6 

73.7 

.0 

.0 

.0 

73.7 

400. 

52.1 

.0 

76.8 

76.8 

.0 

.0 

.0 

76.8 

500. 

57.7 

.0 

79.7 

79.8 

.0 

.0 

.0 

79.8 

630. 

63.1 

.0 

82.8 

82.9 

.0 

.0 

.0 

82.9 

BOO. 

6B.3 

.0 

86.1 

86.1 

.0 

.0 

.0 

86.1 

1000. 

72. B 

.0 

80.9 

89.0 

.0 

.0 

.0 

89.0 

1250. 

76.9 

.0 

91.9 

92.1 

.0 

.0 

.0 

92.1 

1600. 

Bl.l 

.0 

95.3 

95.5 

.0 

.0 

.0 

95.5 

2000. 

B4.5 

.0 

97.8 

98.0 

.0 

.0 

.0 

98.0 

2500. 

87.6 

.0 

100.6 

100.8 

.0 

.0 

.0 

100.8 

3150. 

R0.4 

.0 

103.6 

103.8 

.0 

.0 

.0 

103.8 

4000. 

92.9 

.0 

100.5 

101.2 

.0 

.0 

.0 

101.2 

5000. 

94.9 

.0 

97.6 

99.5 

.0 

.0 

.0 

99.5 

6300. 

96.6 

100.9 

94.6 

103.0 

.0 

.0 

.0 

103.0 

BOOO. 

9B.0 

.0 

91.5 

98,0 

.0 

.0 

.0 

90.8 

10000. 

9B.9 

.0 

OB. 5 

99.3 

.0 

.0 

.0 

9«,3 

OASPL 

104.1 

100.9 

108.4 

110.3 

.0 

.0 

.0 

1)0.3 


WARNING - IN JET MIXING NOISE* ATTEMPT TO EXPAND LOW FREQ RANGE TO 50 HZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 


WARNING - FOR SHOCK BROADBAND MOIST* SOURCE = 7 * Tl IF STROUIIAl I’AP AMI Tf I* 
OF 1/3 O.B. NO. 24 IS ODTSIDI HANOI Ol HASTE R «.l*l C TRA 
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WARNING 


WARNING 


IN JET MIXING NOISE* ATTEMPT TO EXPAND LOW FREQ RANGE TO 50 MZ 
EXCEEDED LIMITS OF THOMSON SPECTRUM SHAPES 


FOR SHOCK BROADBAND NOISE* SOURCE = 7 * THE STROUHAL PARAMETER 
OF 1/3 O.B. NO. 24 IS OUTSIDE RANGE OF MASTER SPECTRA 
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JETSTAR CRUISE NOISE CHECK CASE 
TP’S 3A»3B»3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION: ALT = 36000. MA = .800 POWER SET. CODE = 2.10 

OBSERVER LOCATION NO. = 3 

TOTAL ENGINE AND AIRFRAME NOISE 


1/3 OCT. BAND SPL’S AND OASPL’S DB RE 0.00002 N/M2, 


FPEO 

FAN 

FAN 

FAN 




vJET 

MIXING 

NOISE 

SHOCK 

BB 

SHOCK 

SCREECH 

AIRFRAME 


(HZ) 

FWD 

AFT 

TOTAL 

CMX 

TURB 

CORE 

PR I 

SEC 

TOTAL 

PR I 

SEC 

PR I 

SEC 

TEN 

TBLN 

TOTAL 

'50. 

51.8 

.0 

51.8 

.0 

49.6 

79.2 

74.6 

.0 

74.6 

109.4 . 

.0 

.0 

.0 

.0 

.0 

109.4 

63. 

57.9 

.0 

57.9 

.0 

52.7 

83.2 

76.5 

.0 

76.5 

113.0 

.0 

.0 

.0 

.0 

.0 

113.0 

00. 

61.0 

.0 

61.0 

.0 

56.7 

87.2 

78.5 

.0 

78.5 

116.3 

.0 

.0 

.0 

.0 

.0 

116.3 

^10(1. 

63.9 

.0 

63.9 

.0 

61.4 

90.2 

80.3 

.0 

80.3 

119.1 

.0 

.0 

.0 

.0 

.0 

119.1 

125. 

66.0 

.0 

66.8 

.0 

62.7 

93.2 

82.1 

.0 

82.1 

121.7 

.0 

.0 

.0 

.0 

.0 

121.7 

160. 

70.0 

.0 

70.0 

.0 

64.5 

96.2 

84.0 

.0 

84.0 

123.5 

.0 

.0 

.0 

.0 

.0 

123.5 

-20(1 . 

72.9 

.0 

72.9 

.0 

66.5 

98.7 

05.7 

.0 

85.7 

124.3 

.0 

.0 

.0 

.0 

.0 

124.3 

250. 

75.9 

.0 

75.9 

.0 

68.1 

100.7 

86.9 

.0 

86.9 

125.1 

.0 

.0 

.0 

.0 

.0 

125.1 

315. 

70.9 

.0 

78.9 

.0 

70.2 

102.2 

88.1 

.0 

88.1 

128.0 

.0 

.0 

.0 

.0 

.0 

128.0 

100. 

82.1 

.0 

82.1 

.0 

73.0 

103.2 

09.4 

.0 

B9.4 

135.7 

.0 

.0 

.0 

.0 

.0 

135.7 

500. 

85.0 

.0 

85.0 

.0 

74.7 

102.2 

89.7 

.0 

89.7 

142.3 

.0 

.0 

.0 

.0 

.0 

142.3 

630. 

00.2 

.0 

80.2 

.0 

76.8 

100.7 

90.0 

.0 

90.0 

143.8 

.0 

.0 

.0 

.0 

.0 

143.8 

‘'BOO. 

91.5 

.0 

91.5 

.0 

70.7 

98.7 

90.4 

.0 

90.4 

141.6 

.0 

.0 

.0 

.0 

.0 

141.6 

1000. 

91.1 

.0 

91.4 

.0 

81.0 

96.2 

91.6 

.0 

91.6 

141.2 

.0 

.0 

.0 

.0 

.0 

141.2 

1250. 

97.5 

.0 

97.5 

.0 

82.7 

93.2 

92.8 

.0 

92. B 

140.4 

.0 

.0 

.0 

.0 

.0 

140.4 

'1600 • 

100.9 

.0 

100,9 

.0 

85.2 

90.2 

94.1 

.0 

94.1 

139.6 

.0 

.0 

.0 

.0 

.0 

139.6 

2000. 

103.1 

.0 

103.4 

.0 

88.0 

87.2 

94.6 

.0 

94.6 

138.3 

.0 

.0 

.0 

.0 

.0 

138.3 

2500. 

106.3 

.0 

106.3 

.0 

89.8 

83.2 

95.2 

.0 

95.2 

137.3 

.0 

.0 

.0 

.0 

.0 

137.3 

‘3150. 

109.3 

.0 

109.3 

.0 

91.9 

79.2 

95.7 

.0 

95.7 

136.2 

.0 

.0 

.0 

.0 

.0 

136.2 

1000. 

107.2 

.0 

107.2 

.0 

94.2 

75.7 

95.5 

.0 

95.5 

134.8 

.0 

.0 

.0 

.0 

.0 

134.8 

‘5000 • 

106.3 

.0 

106.3 

.0 

95.9 

71.7 

95.2 

.0 

95.2 

133.5 

.0 

.0 

.0 

.0 

.0 

133.5 

6300. 

111.1 

.0 

111.4 

.0 

97.4 

67.2 

95.0 

.0 

95.0 

132.2 

.0 

.0 

.0 

.0 

.0 

132.3 

0000. 

107.2 

.0 

107.2 

.0 

98.0 

63.2 

94.3 

.0 

94.3 

131.0 

.0 

.0 

.0 

.0 

.0 

131.0 

lonou. 

107.9 

.0 

107.9 

.0 

99.9 

58.2 

93.6 

.0 

93.6 

130.0 

.0 

.0 

.0 

.0 

.0 

130.0 

OASPL 

117.2 

.0 

117.2 

.0 

105.2 

110.4 

105.7 

.0 

105.7 

150.8 

.0 

.0 

.0 

.0 

.0 

150.8 


123 



JETSTAR CRUISE NOISE CHECK CASE 
TP’S 3A»3B»3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITION! ALT = 36000. MA = .BOO POWFR SET. CODE = 2.10 
OBSERVER LOCATION NO. = 3 

ENGINE NO. 1 

FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL’S - DB RE 0.00002 N/M2 


** EXHAUST DUCT ** 
*** INLET RADIATED NOISE *** * RADIATED NOISE * 


FREQ 

(HZ) 

BROAD 

-BAND 

DSCRT 

TONES 

CMBNTN 

TONES 

TOTAL 

BROAD 

-BAND 

DSCRT 

TONES 

TOTAL 

TOTAL 

FAN 

NOISE 

50. 

.0 

.0 

52.5 

52.5 

.0 

.0 

.0 

52.5 

63. 

.0 

.0 

55.5 

55.5 

.0 

.0 

.0 

55.5 

BO. 

8.9 

.0 

58.6 

58.6 

.0 

.0 

.0 

58.6 

100. 

16,4 

.0 

61.5 

61.5 

.0 

.0 

.0 

61.5 

125. 

24.1 

.0 

64.4 

64.4 

.0 

.0 

.0 

64.4 

160. 

32.2 

.0 

67.7 

67.7 

.0 

.0 

.0 

67.7 

200. 

39.2 

.0 

70.6 

70.6 

.0 

.0 

.0 

70.6 

250. 

45.9 

.0 

73.5 

73.5 

.0 

.0 

.0 

73.5 

315. 

52.4 

.0 

76.5 

76.5 

.0 

.0 

.0 

76.5 

400. 

58.8 

.0 

79.7 

79.7 

.0 

.0 

.0 

79.7 

500. 

64.3 

.0 

B2.6 

82.7 

.0 

.0 

.0 

82.7 

630. 

69.8 

.0 

85.7 

85.8 

.0 

.0 

.0 

85.8 

800. 

75.0 

.0 

88.9 

89.1 

.0 

.0 

.0 

89.1 

1000. 

79,4 

.0 

91.8 

92.0 

.0 

.0 

.0 

92.0 

1250. 

83.6 

.0 

94.8 

95.1 

.0 

.0 

.0 

95.1 

1600. 

87.8 

.0 

98.2 

98.6 

.0 

.0 

.0 

98.6 

2000. 

91.2 

.0 

100.6 

101.1 

.0 

.0 

.0 

101.1 

2500. 

94.3 

.0 

103.5 

104.0 

.0 

.0 

.0 

104.0 

3150. 

97.1 

.0 

106.5 

107.0 

.0 

.0 

.0 

107.0 

4000. 

99.6 

.0 

103.4 

104.9 

.0 

.0 

.0 

104.9 

5000. 

101.6 

.0 

100.5 

104.1 

.0 

.0 

.0 

104.1 

6300. 

103.3 

107.6 

97.4 

109.2 

.0 

.0 

.0 

109.2 

0000. 

104.6 

.0 

94.3 

105.0 

.0 

.0 

.0 

105.0 

10000. 

105.5 

.0 

91.4 

105.7 

.0 

.0 

.0 

105.7 

OASPL 

110.8 

107.6 

111.3 

114.9 

.0 

.0 

.0 

114.9 



JETSTAR CRUISE NOISE CHECK CASE 
TP’S 3A.3B.3C MIC LOCATIONS 1 AND 3 


FLIGHT CONDITIONS ALT = 36000. MA = .000 POWER SET. CODE = ?.10 

OBSERVER LOCATION NO. = 3 

ENGINE NO. 2 

FAN NOISE COMPONENTS 

1/3 OCT. BAND SPL'S AND OASPL’S - D8 RE 0.00002 N/M2 

** EXHAUST DUCT ** 

*** INLET RADIATED NOISE *** * RADIATED NOISE * 


TOTAL 


FREQ 

BROAD 

OSCRT 

CMBNTN 


BROAD 

DSCRT 


FAN 

(HZ) 

-BAND 

TONES 

T0NE5 

TOTAL 

-BAND 

TONES 

TOTAL 

NOISE 

50. 

.0 

.0 

51.0 

51.0 

.0 

.0 

.0 

51.0 

63. 

.0 

.0 

54.1 

54.1 

.0 

.0 

.0 

54.1 

BO. 

7.2 

.0 

57.2 

57.2 

.0 

.0 

.0 

57.2 

100. 

14.6 

.0 

60.1 

60.1 

.0 

.0 

.0 

60.1 

125. 

22.2 

.0 

63.0 

63.0 

.0 

.0 

.0 

63.0 

160. 

30.3 

.0 

66.2 

66.2 

.0 

.0 

.0 

66.2 

200. 

37.3 

.0 

69.1 

69.1 

.0 

.0 

.0 

69.1 

250. 

44.0 

.0 

72.1 

72.1 

.0 

.0 

.0 

72.1 

315. 

50.5 

.0 

75.1 

75.1 

.0 

.0 

.0 

75.1 

000. 

56.9 

.0 

70.2 

78.3 

.0 

.0 

.0 

78.3 

500. 

62.4 

.0 

81.2 

81.2 

.0 

.0 

.0 

81.2 

630. 

67.8 

.0 

84.3 

84.4 

.0 

.0 

.0 

84.4 

BOO. 

73.0 

.0 

87.5 

87.6 

.0 

.0 

.0 

87.6 

1000. 

77.5 

.0 

90.4 

90.6 

.0 

.0 

.0 

90.6 

1250. 

81.7 

.0 

93.4 

93.6 

.0 

.0 

.0 

93.6 

1600. 

85.9 

.0 

96.7 

97.1 

.0 

.0 

.0 

97.1 

2000. 

09.3 

.0 

90.2 

99.6 

.0 

.0 

.0 

99.6 

2500. 

92.4 

.0 

102.0 

102.5 

.0 

.0 

.0 

102.5 

3150. 

95.2 

.0 

105.0 

105.5 

.0 

.0 

.0 

105.5 

4000. 

97.7 

.0 

101.9 

103.3 

.0 

.0 

.0 

103.3 

5000. 

99.7 

.0 

9°.0 

102.4 

.0 

.0 

.0 

102.4 

6300 . 

101.4 

105.7 

96.0 

107.4 

.0 

.0 

.0 

107.4 

BOOO. 

102.7 

.0 

92.9 

103.1 

.0 

.0 

.0 

103.1 

10000. 

103.6 

.0 

90.0 

103.0 

.0 

.0 

.0 

103.0 

OASPL 

108.9 

105.7 

109.0 

113.2 

.0 

.0 

.0 

113.2 
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COMPARISONS WITH MEASURED DATA 


One of the objectives of the present effort was to compare acoustic 
measurements made in flight with predictions of the Cruise Noise Pre- 
diction Program. As a result of these comparisons, it was planned to 
modify and improve the prediction program where necessary and where 
possible with the available data. This effort, as anticipated, was 
significantly constrained by time and budget. 

The measured data provided for these comparisons consist of 50 Hz band- 
width sound pressure levels recorded at 5 microphone locations on the 
NASA/Dryden JetStar. The data covered the frequency range of the pre- 
dictions (50 to 10,000 Hz) and beyond to about 20,000 Hz. Three of the 
microphones were mounted approximately 3 feet apart (spanwise) on the 
upper surface leading edge section just outboard of the wing slipper 
tank. For reference herein, these microphones are numbered 1 to 3 
beginning with the outboard microphone. One of the remaining micro- 
phones was located on the slipper tank and the other inside the inlet 
of engine number 1. The data from the slipper tank microphone were 
reported to be erroneous and were not used for comparisons. Data were 
supplied for three test flights and about 18 test conditions per flight. 
The test points (altitude, Mach number, power setting, C , etc) were 
essentially the same for each flight. They covered the range from 
ground, static operation to cruise at 40,000 feet. Spectral data were 
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supplied for only one of the three flights so the data of this flight 
were selected for comparison purposes. The overall sound pressure 
levels (reported for all flights) indicated no significant difference 
in measured noise levels among the three flights. 

Cruising flight at 36,000 feet was, somewhat arbitrarily, selected as 
the baseline flight condition. Figure 5 shows the measured and predicted 
acoustic signature at microphone number 1 for the baseline case at an 
engine pressure ratio of 1.9. The predicted levels shown are 1/3 octave- 
band sound pressure levels corrected for bandwidth to obtain 50 Hz band- 
width levels. As shown, the, spectrum shapes are similar but the pre- 
dicted levels run about 15 db higher than measured. The predicted com- 
ponent noise levels for this case are shown in Figure 6. Jet shock 
broadband noise is seen to dominate the predicted noise spectrum.- 
Figure 7 provides a comparison for the baseline case at a higher cruise 
Mach number and power setting. For this flight condition, the predicted 
and measured spectra have moved farther apart. Also, note that the 
measured compressor tone which showed reasonable agreement with pre- 
dicted tones in Figure 5 has disappeared in Figure 7. Comparisons 
made at other flight conditions and microphone locations on the wing 
show similar results. As a point of reference, Figure 8 presents a 
comparison for static, ground operation. Under these conditions, the 
comparison shows fairly good agreement. At the present, the reasons 
for the differences seen here in the cruise environment have not been 
identified. It was hoped that Mach number and power setting effects 
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could be derived from the measured data. Unfortunately, it appears 
that the flight cases which might provide an isolated effect of either 
parameter have turbulent flow occurring at the microphone locations. 

The flight data included a notation as to the occurrence of laminar or 
turbulent flow at each microphone location for each test point. The 
measured noise spectra seem to confirm these observations. This, in 
effect, has "washed-out" the noise signature impinging from outside the 
boundary layer, i.e. the acoustic record appears to be dominated by the 
fluid pressure fluctuations of the turbulent flow. This problem also 
occurs in attempting to compare cases with and without shocks (or 
supercritical nozzle pressure ratios). 

In an alternate approach, perturbations to the prediction procedures 
have been made in an attempt to produce agreement with measured data. 

The results of this investigation are inconclusive. For example, 
removing the shock noise or the convective amplification effect from the 
predictions removes most of the differences seen in the noise comparisons 
above. However, as noted, these individual effects are not verifiable 
from the measured data. Consequently, it is suggested that more flight 
noise data be obtained with an improved (but as yet undefined) noise 
data acquisition system. Hopefully this will occur eventually. 

Several areas have been identified as possible contributors to the 
apparent lack cf agreement seen in these comparisons. These areas have 
been explored to only a limited extent due to budget constraints. They 
include (1) the shock noise prediction, (2) the convective amplification 
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correction, (3) the altitude correction, and (4) shielding effects. 
These, of course, are not the only possibilities, but appear to be the 
most obvious or promising in terms of producing explanations for the 
results seen here. Another area of interest is the compressor noise 
and its tones as noted previously. A brief discussion of these aspects 
of the prediction procedure follow. 

Jet shock-associated broadband noise clearly dominates the predicted 
noise spectra of the JetStar at cruise conditions. Unfortunately, as 
noted previously, it appears that the noise spectra from test points 
without shocks, which might allow source separation of shock noise, 
are not available. Otherwise, there is no clear reason to suspect the 
basic shock prediction procedure. High forward speed, as mentioned 
before, is not expected to significantly affect the shock noise source 
strength or frequency content, but this remains to be shown. Another 
possibility is shielding of shock noise by the engine nacelle. If the 
source were located nearer the nozzle exit than predicted, as may be 
possible with the JetStar nozzle, significant shielding could occur. 
This is especially true since the source emission angles for radiation 
to the wing microphone positions at cruise are very small. 

Predictions were generated with shock noise excluded to show how the 
comparisons might look if shock noise was being severely over predicted 
at the measurement locations. The results are shown in Figure 9. In 
the absence of shock noise, the predictions run, on the average, only 
a few dB below the measured data. However, the trends with increasing 
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Mach number and power setting are still reversed. As discussed below, 
aerodynamic shielding may be a factor in this difference in trends. 

Convective amplification is a phenomenon associated with motion of an 
acoustic source relative to the surrounding medium. Its effect on air- 
craft noise radiation patterns has been demonstrated but only for Mach 
numbers up to about 0.25. In the absence of data for higher Mach 
numbers, the correction for this effect is assumed to apply at transonic 
Mach numbers as well. This results in rather large corrections to the 
predicted noise levels at cruise conditions, e.g. an increase of 28 dB 
in front of a source and a 10 dB reduction behind a source at 0.80 
Mach number. Obviously, this extrapolation to cruise Mach numbers 
needs experimental verification or redefinition. Until these data 
become available, this correction must be viewed with some skepticism. 

As noted above, the present data do not appear to provide valid aircraft 
noise signatures with aircraft Mach number as the only variable. As 
was done in the case of shock noise, several computer runs were made 
excluding the correction for convective amplification. The results 
are seen in Figure 9, and continue to show an increase in level with 
Mach number and power setting. The levels agree with the measured data 
at 0.70 Mach number but overall average about 4 to 5 dB higher than 
measured . 

The effect of cruise altitude has been accounted for in the present 
procedure by applying a sound level correction based on the ratio of 
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the local characteristic acoustic impedance (pc) to that at sea level, 

standard day (pc). The correction applied is 10 log„„ ( pc/ pc ). 

oo 10 o o 

However, it now appears that this correction procedure requires further 

investigation and possible modification, especially in light of the 

present comparisons. The correction which may be more appropriate here 

is in fact given in Section 3.4.2 of Reference 1. This correction, 

2 2 

20 log„„ ( pc / p c ), will result in about a 7 dB reduction in the 
10 r o o 

program predictions at 36,000 ft. However, it may turn out that the 
proper correction will be a function of the nature of the noise source 
(monopole, dipole, or quadrapole). Alternatively, the problem may be 
resolved by further investigation into the applicability of the basic 
prediction procedures to the high altitude environment. It is antici- 
pated that further investigation will show the later correction above 
to be more applicable to cruise noise predictions and, consequently, a 
correction to the program would be in order. 

Shielding effects are believed to play a key role in the differences 
between measured and predicted noise levels seen in the present study. 
The most influential shielding effects, however, are not structural but 
appear to be the result of the boundary layer and induced flow field 
over the wing. Previous JetStar aerodynamic studies and data in the 
literature indicate that the supersonic flow region over the upper wing 
surface may extend to at least a height equivalent to about 20% of 
wing chord. This effect may be expected to result in a significant 
shielding effect with the shielding effect increasing as the Mach 
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number goes up. Considering the geometry of the engine/microphone loca- 
tions for the present case, this may be at least a partial explanation 
for the observed decrease in the measured noise levels with increasing 
Mach number and power setting and the observed differences in measured 
and predicted noise levels. Unfortunately, methods for predicting this 
effect are not presently available. Furthermore, as stated previously, 
the present procedure was developed for free-field conditions. The 
effects of refraction, for example, in the wing flow field and boundary 
layer are not included. Even in the absence of a supersonic flow region 
on the wing, for our engine/wing/microphone arrangement, refraction in 
the boundary layer can be expected to produce some shielding. Consid- 
ering the complexity of the mechanisms involved, which are not predict- 
able with any confidence with available procedures, and their potential 
importance, the structural shielding procedure of Reference 2 does not 
seem applicable. Consequently, a shielding routine has not been devel- 
oped for the prediction program. 

In conclusion, the results of the comparisons made to date clearly show 
that further investigations are required to explain the differences 
seen in the measured and predicted noise levels. It is believed, at 
this point, that the experimental isolation of the effects discussed 
in this section will go a long way toward resolving these differences, 
and that the prediction program presented here will form an adequate 
basis for aircraft preliminary design work. 
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Figure 5. Comparison of measured and predicted noise levels for Lockheed JetStar at 
" ‘ cruise conditions. _l 
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Figure 8. Comparison of measured and predicted noise levels for Lockheed JetStar during static, 
ground operation. 
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